T • H • E 


OHD 

skie 

UNIVERSITY 


Vibration T ransmission 
Through Rolling Element Bearings 
In Geared Rotor System 


Rajendra Singh 

Department of Mechanical Engineering 
Teik Chin Lim 

Graduate Research Associate 


National Aeronautics and Space Administration 

Lewis Research Center 
Cleveland, Ohio 44135 

Grant No. NAG 3-773 
Final Report - Part I 
RF Project 765863/719176 


December 1989 


ACKNOWLEDGEMENTS 


We are very grateful for the financial support provided throughout this Ph.D. 
research by the NASA Lewis Research Center, Department of Mechanical Engineering, 
and Graduate School of the Ohio State University. Appreciation is also extended to J.S. 
Lin, D.R. Houser, J.J. Zakrajsek and F.B. Oswald for providing some of the 
experimental results used in Chapters III-V. 


1 


TABLE OF CONTENTS 


ACKNOWLEDGEMENTS i 

LIST OF TABLES vi 

LIST OF FIGURES viii 

LIST OF SYMBOLS xvii 

ABSTRACT xxvii 

CHAPTER PAGE 

I INTRODUCTION 1 

1.1 Introduction 1 

1.2 Literature Review 6 

1.3 Scope and Objectives 7 

II BEARING STIFFNESS FORMULATION 12 

2.1 Introduction 12 

2.2 Literature Review 14 

2.3 Assumptions and Objectives 16 

2.4 Bearing Load-Displacement Relations 18 

2.5 Development of Bearing Stiffness Matrix [K]^ 24 

2.6 Numerical Estimation of [K]jj m 31 

2.7 Validation of Proposed Model 33 

2.8 Parametric Studies 36 

2.9 Conclusions 46 

ii 



HI BEARING SYSTEM STUDIES 51 

3.1 Introduction 51 

3.2 Literature Review 53 

3.3 Assumptions and Objectives 54 

3.4 System Governing Equations 56 

3.4.1 Method A: Lumped Parameter Model 56 

3.4.2 Method B: Dynamic Finite Element Formulation 57 

3.4.3 Other Methods 58 

3.5 Bearing System Stability 60 

3.6 System Reponse 63 

3.7 Example Case I: Rigid Shaft and Plate System 71 

3.7.1 Vibration Models 71 

3.7.2 Bearing Transmissibility 75 

3.8 Example Case II: Rigid Shaft and Plate Supported on Flexible Mounts ... 80 

3.8.1 Vibration Models 80 

3.8.2 Frequency Reponse 84 

3.9 Example Case HI: Experimental Study 88 

3.9.1 Physical Setup 88 

3.9.2 Bearing Analysis 88 

3.9.3 System Study 91 

3.10 Concluding Remarks 96 

IV GEARED ROTOR SYSTEM STUDIES 97 

4.1 Introduction 97 

4.2 Literature Review 101 

4.2.1 Casing and Mount Dynamics 101 

4.2.2 Gear Dynamics 102 

4.2.3 Rotor Dyamics 102 

4.3 Assumptions and Objectives 103 

4.4 Theory 105 


in 


4.4.1 Method A: Lumped Parameter Model 105 

4 .4 . 1 . 1 Equations of Motion 105 

4.4. 1.2 System Kinetic Energy 107 

4.4. 1.3 Shaft Stiffness Matrix 107 

4. 4. 1.4 Gear Mesh Stiffness Matrix 110 

4.4. 1.5 Flexible Mount Stiffness Matrix Ill 

4.4. 1.6 Bearing Stiffness Matrix 112 

4.4.2 Method B: Dynamic Finite Element Formulation 113 

4.4.3 Other Methods 113 

4.5 Example Case I: Single-stage Rotor System with Rigid Casing 

and Flexible Mounts 114 

4.5.1 Vibration Models 114 

4.5.2 Eigensolution 118 

4.5.3 Transmissibility 126 

4.6 Example Case II: Geared Rotor System with Rigid Casing 

and Flexible Mounts 131 

4.6.1 Bearing Analysis 131 

4.6.2 Vibration Models 133 

4.6.3 Eigensolution 136 

4.6.4 Transmissibility Spectra 140 

4.7 Example Case HI: Geared Rotor System with Rigidly 

Mounted Flexible Casing 144 

4.7.1 Physical Setup 144 

4.7.2 Vibration Models 146 

4.7.3 Casing Response 147 

4.8 Concluding Remarks 153 

V STATISTICAL ENERGY ANALYSIS 155 

5.1 Introduction 155 

5.2 Modal Analysis of Gear Casing and Mounts 156 

5.2.1 Finite Element Model 156 

5.2.2 Experiments and Modal Validation 158 


tv 


5.2.3 Parametric Studies 160 

5.3 Justification for Using SEA Approach 165 

5.3.1 Modal Densities 165 

5.3.2 Literature Review 166 

5 .4 Example Case I: Coupling Loss Factor of Plate-Cantilevered 

Beam System 167 

5.5 Example Case II: Coupling Loss Factor of Circular 

Shaft-Bearing-Plate System 169 

5.6 Example Case III: A Circular Shaft-Bearing-Plate System 175 

5.6.1 Theory 175 

5.6.2 Validation and Parametric Studies 182 

5.7 Example Case IV: A Geared Rotor System 188 

5.7.1 Assumptions : 188 

5.7.2 Coupling Loss Factor 190 

5.7.3 Vibroacoustic Response 192 

5.7.4 Experimental Validation and Parametric Studies 194 

5.8 Concluding Remarks 200 

VI CONCLUSIONS AND RECOMMENDATIONS ! 201 

6.1 Summary 201 

6.2 Contributions ; 205 

6.2 Future Research 206 

LIST OF REFERENCES 208 


v 



LIST OF TABLES 


Table Tide Page 

2. 1 Design parameters for typical ball and roller bearings used for 

parametric studies 38 

2 . 2 Comparison between the proposed and current ball bearing 

stiffness coefficients 47 

2.3 Comparison between the proposed and current roller bearing 

stiffness coefficients 48 

3 . 1 Bearing system stability criteria for the proposed ball bearing model 64 

3.2 Bearing system stability criteria for the proposed roller bearing model 65 

3 . 3 Design parameters for typical ball and roller bearings used for 

system studies 69 

3.4 Bearing stiffness coefficients and undamped natural frequencies (Hz) 

of example case I 75 

3.5 Results of example case II 83 

3.6 Design and estimated parameters for two identical rolling element 

bearings used in example case III 91 


vi 


3 .7 Predicted and measured cross point mobility level Ly as defined 

by equation (3.22) : . . . 95 

4. 1 Design parameters for deep groove ball bearings used for example 

cases 115 

4 . 2 System parameters of example case I: single-stage rotor system shown 

inFigure4.2 119 

4.3 Comparison of natural frequencies in Hz for example case 1 119 

4.4 Effect of bearing preload F zbm on coj (Hz) for example case I 124 

4.5 System parameters of example case II: geared rotor system shown 

in Figure 4.1 132 

4.6 Eigensolution of example case II 137 

4.7 Example case II: effect of mean bearing displacement 5^ on tOj (Hz) 140 

4.8 Comparison of L A (dB) for example case in 151 

5 . 1 Number of modes for the stiffened and unstiffened gearboxes 

mounted rigidly 163 


5.2 Typical bearing stiffness coefficients of Figure 5.9 for example case II ... . 175 


vu 


LIST OF FIGURES 


Figure Title Page 

1.1 Schematic of a typical shaft-bearing-plate-mount system 3 

1 . 2 Schematic of typical rolling element bearings, (a) Bearing forces and 
moments, (b) Deep groove ball bearing, (c) Angular contact ball bearing. 

(d) Taper roller bearing, (e) Cylindrical roller bearing 4 

1.3 Generic rotating mechanical systems, (a) A geared rotor system. 

(b) A single-stage rotor system 5 

1 .4 Vibration energy transmission in geared rotor system, (a) Structure-borne 

noise paths, (b) Statistical energy analysis model 10 

2.1 Schematic representation of the vibration transmission problem. Here 

the flexible shaft is subjected to mean forces F wsm and torques T wsm 
where w = x, y or z, is the direction and subscript m and s implies mean 
and shaft respectively. Also, 9 is the angular displacement and u is the 
translational displacement 13 

2.2 Rolling element bearing kinematics and coordinate system. Here the 
following nomenclature is used: djj 0 is the outer raceway diameter, dj, m 
is the bearing pitch diameter, dy is the inner raceway diameter, \| f is the 
angular position of rolling element, 5^ is the mean translational 
displacement, P pm is the mean angular displacement, F w (, m is the mean 


vrn 


bearing force, and M pbm is the mean bearing moment where w=x,y,z, 

and p=x ,y , are the directions ; . . . 17 

2.3 Elastic deformation of rolling element for non-constant contact angle ctj 
given by the change in the distance between the inner a, and outer a*, 
raceway groove radius curvature centers due to the mean bearing loads 

or displacements 20 

2.4 Elastic deformation of rolling element for constant contact angle ccj = a 0 
given by the change in the relative position of the inner and outer 

raceways due to the mean bearing loads or displacements 21 

2.5 Decomposition of the effective radial (8)^ and axial (8) z j deformations of 

the j-th rolling element in terms of the mean bearing displacements [q} bm . 
Here G is the bearing outer ring geometrical center 22 

2.6 Comparison between the proposed theory and Gargiulo's formulas [14] for 
axial k bzz and radial k brr stiffness coefficients of ball and roller bearings ... 34 

2.7 Comparison between the proposed theory and the experimental results of 

Kraus et al.[24] for k brr as a function of the mean axial preload 35 

2.8 Comparison between the proposed theory k brr , estimated k brr , and White's 

analytical results [20]. a. Ball bearing, b. Roller bearing 37 

2.9 Dominant stiffness coefficients of ball bearing set A for 0° <; a 0 <; 90° and 

given a constant mean radial bearing displacement 8^ = 0.025 mm 39 

2.10 Dominant stiffness coefficients of roller bearing set B for 0° < a D <, 90° and 

given a constant mean radial bearing displacement 8^ = 0.025 mm 40 

2.11 Dominant stiffness coefficients of ball bearing set A for 0° < a 0 < 90° and 

given a constant mean axial bearing displacement 8^ = 0.025 mm 42 


IX 



2.12 Dominant stiffness coefficients of roller bearing set B for 0° < a 0 < 90° and 

given a constant mean axial bearing displacement 8^, = 0.025 mm 43 

2.13 Dominant stiffness coefficients of ball bearing set A for 0° < a 0 < 90° and 

given a constant misalignment Pyn, = 0.015 rad 44 

2.14 Dominant stiffness coefficients of roller bearing set B for 0° < Oq < 90° and 

given a constant misalignment Pyn, = 0.015 rad 45 

3 . 1 Schematic representation of the vibration transmission problem. Here 

the flexible shaft is subjected to alternating forces Fj^ft) and torques TjsaW 
where j = x, y or z, is the direction and subscript a implies alternating. Also, 

0 is the angular displacement and u is the translational displacement 52 

3 . 2 Discretization method for implementation of [K] bm in finite element model 

of a system similar to Figure 3.1. Here {q(t)} sa and { q(t) } ca are the alternating 
shaft and plate displacement vectors at bearing location respectively. 

Subscript a implies alternating component 59 

3 . 3 Plot of stability criteria functions for ball and roller bearings subjected to 

mean bearing radial deflection 8^. (a) SI = Oj(2,4) and S2 = <E>j(l,3). 
(b)S3=k bll 0 1 (3,5) + 0 2 (35,13,15,33,44) 66 

3.4 Plot of stability criteria functions for ball and roller bearings subjected 
to mean bearing axial deflection 8^ Here, SI = Oj(2,4) 

and S4 = 0^1,5) 67 

3 . 5 Plot of stability criteria functions for ball and roller bearings subjected 

to mean angular misalignment p ym . (a) SI = 0^2,4) and S2 = <£>2(1,3). 

(b) S3 =k bll <D 1 (3,5)+0 2 (35,13,15,33,44) 68 

3.6 Example case I: freely suspended rigid shaft, ball bearing and rigid plate 
system subjected to alternating radial force F^ft) applied at the shaft. 


x 


(a) Physical system, (b) Proposed multi-degree of freedom vibration model 
with DOF = 12. (c) Simple model by Kraus et al. [24] with DOF = 1 72 

3.7 Bearing transmissibility spectra ^f wba ,F xsa ( <0 o) f° r example case I. 

Here, Rl: force transmissibility with f wba = F xba and R2: moment 
transmissibility with f wba = M yba , as predicted by our model with DOF=2 
and the simple model by Kraus et al. [24] with DOF =1 78 

3 . 8 Effect of mean axial preload on the bearing transmissibility spectra 

Rf p (co 0 ) defined by equation (3.17) for example case I. (a) Force 

woa’ xsa 

transmissibility with f wba = F xba . (b) Moment transmissibility with 

fwba = Myba 79 

3.9 Example case II: rigid shaft, ball bearing and rigid plate system supported 
by flexible mounts and subjected to alternating radial force F^t) applied 
at the shaft, (a) Physical system, (b) Proposed multi-degree of freedom 
vibration model with DOF =12. (c) White's vibration model [20] with 

DOF = 2 81 

3.10 Accelerance spectra Aq w j a ,F xsa ( w o) for example case II as predicted by our 

formulation and White's model [20]. (a) Driving point accelerance with 
q w j a = u Ka . (b) Cross point accelerance with q w j a = u xca . (c) Accelerance 
with q w j a = u^jj for shaft and q w j a = u Ka for casing, (d) Accelerance with 
q w j a = 9ysa f° r shaft and q w j a = 9 yca for casing 85 

3.11 Bearing Rf wba ,F xga ( 03 o) and mount Rf wva ,F xsa ( t0 o) transmissibility spectra 

for example case II as predicted by our formulation and White's model [20]. 

(a) Bearing force transmissibility with f wba = F xba . (b) Mount force 
transmissibility with = F xva . (c) Force transmissibility with f wba = F zba 
for bearing and f wva = F wa for mount, (d) Moment transmissibility with 
fwba = M y ba for bearing and f wva = Myv a for mount 86 


xi 


3.12 Effect of the mean radial bearing force F xbm on the bearing transmissibility 
spectra Rf wba ,F xsa (®o) ^ or example case n. (a) Force transmissibility with 

f wba = F xba . (b) Force transmissibility with f wba = F zba . (c) Moment 
transmissibility with f wba = M yba 87 

3.13 Example case m. (a) Schematic of the experimental setup consisting of an 
overhung shaft, 2 ball bearings and a rectangular plate [17]. (b) Static model 
of shaft and bearings used for computing [K] bm . (c) Finite element model 
of the experimental system shown in (a). The generalized stiffness matrix 
elements for the bearing model are shown by dashed lines, connecting 1 node 


on the shaft to 4 nodes on the plate 89 

3.14 Driving point accelerance spectra yielded by the proposed 

model, simple model and experiment for example case HI 93 

3.15 Cross point accelerance spectra p (a> 0 ) yielded by the proposed 


model and experiment for example case III. In this case, the simple model 
predicts zero response u zca at the plate, (a) Point Cl (see Figure 3.13c). 

(b) Point C2 (see Figure 3.13c) 94 

4. 1 Schematic of a geared rotor system with flexible shafts and rigid spur gear 
and pinion, and supported by 4 rolling element bearings on a flexibly 
mounted casing. The shaft at the motor end is driven at mean speed Q. ^ 

which in turns drives the load at mean speed £2 zL - This physical system is 

symmetric about a plane intersecting both the driving and driven shafts. 

The y axis is assumed parallel . to the gear mesh line of action in the pressure 
angle direction at mesh point. This schematic is used for example case II 
and III 98 

4.2 Schematic of the single stage rotor system with flexible shaft rotating at 
mean speed Q z and supported by 2 identical rolling element bearings on a 



flexibly mounted rigid casing. This physical system is symmetric about a 
plane intersecting the shaft axial axis and parallel to the x-z plane. This 
schematic is used for example case I 100 

4.3 Lumped parameter model of the flexible shafts with rigid rotors. Each shaft 
is supported by 2 identical bearings on both ends, (a) Internal rotating system 
of a typical geared rotor system of Figure 1. (b) Internal rotating system of a 


typical single-stage rotor system of Figure 2. (c) Discretization of a shaft 
segment and the degrees of freedom associated with each lumped mass ... 108 

4.4 Natural frequencies 004 and CO 5 versus the number of shaft segments, n s 

for example case I 120 

4.5 First 6 mode shapes of example case I: single-stage rotor system including 

rigid casing and flexible mounts 121 


4.6 Effect of casing mass n^. normalized with respect to shaft and rotor masses, 
and mount stiffness k vx normalized with respect to k bxx on system natural 


frequencies coj for example case I 123 

4.7 Effect of bearing preload F bzm on (£>j for example case I 125 


4.8 Bearing transmissibility for example case I. (a) Force transmissibility 

Rp xba> F xsa (^z)- O 5 ) Moment transmissibility RM yba ,F xsa (^z) which is not 


predicted by simple model 127 

4.9 Mount transmissibility for example case I. (a) Force transmissibility 

Rp p (fl z ). (b) Moment transmissibility R M p (£2 Z ) 128 

xva* xsa yva> xsa 

4.10 Effect of mount stiffness ky* on the force transmissibilities for example case I. 

(a) Bearing, (b) Mount 129 


xin 


4.11 Effect of casing mass on the force transmissibilities for example case I. 

(a) Bearing, (b) Mount :... 130 

4.12 Effect of casing mass m c (normalized with respect to shaft and gear masses) 

and mount stiffness k vx (normalized with respect to k bxx ) on the system 
natural frequencies <Oj for example case II 139 

4.13 Comparison of transmissibility and normalized dynamic transmission error 
spectra for example case n. (a) Bearing force transmissibility Rp^F^ 0 ^)- 

(b) Bearing moment transmissibility RM xba ,F h ( ffl h) which is not predicted by 
simple model, (c) Mount transmissibility RM zba ,F h ( 0) h)- (d) Normalized 

dynamic transmission error p(coj,) 141 

4.14 Mount transmissibility studies for example case n. Effect of (a) casing 
mass m c , (b) mount stiffness and (c) mean bearing displacement 5^ 

0n R M zba ,F h (®h) 142 

4.15 Effect of mean bearing displacement 5^ on (a) bearing force transmissibility 
Rpy ba ,F h (®h)> (b) bearing moment transmissibility RM xba ,F h ( <0 h)» and 

(c) normalized dynamic transmission error pfo^) for example case II 143 

4.16 Example case EH: (a) Schematic of the NASA gear test facility, (b) NASA 
spur gear pair system supported on rigidly mounted flexible casing [54] . . . 145 

4.17 Measurement locations on the casing plate for example case HI 148 

4.18 Comparison of casing flexural vibrations for example case IH with 
$zm = 0.08 mm. (a) Side plate supporting the bearings and facing the 
dynamometer (location 18a). (b) Top plate cover (location 18b). 

(c) Side plate with no bearings (location 18c) 149 

4.19 Sensitivity of predicted casing acceleration to mean bearing displacement for 
example case III. (a) Side plate supporting the bearings and facing the 


xiv 


dynamometer (location 18a). (b) Top plate cover (location 18b). 

(c) Side plate with no bearings (location 18c) 152 

5.1 Gearbox with rigid and flexible mounts 157 

5 . 2 NASA gear test facility [55] 159 

5.3 Mode shapes of the rigidly mounted, stiffened NASA gearbox as 
predicted by FEM (EMA results are given in parenthesis). 

Bold sign implies larger amplitude 161 

5 .4 Number of modes within each one-third octave bands for the NASA gearbox. 

(a) Predicted and measured number of modes for rigidly mounted casing. 

(b) Number of modes for rigid and flexible mounting conditions 162 

5 . 5 Rigid body modes of the flexibly mounted gearbox as predicted by FEM ..164 

5.6 Schematic of example case I: a plate-cantilevered beam, (a) Experimental finite 

system, (b) Theoretical infinite system 168 

5.7 Comparison between Lyon and Eichler's approximation [56,58] and our 
proposed formulation for example case I. (a) Comparison of IZ C I/IZ C +Z S I. 

(b) Comparison of predicted Tj sc with experimental results given by 

Lyon and Eichler [56,58] 170 

5.8 Boundary conditions for example case II: semi-infinite shaft-bearing 

system 171 

5.9 Coupling loss factor T| sc used in example case II for various bearing stiffness 

coefficients listed in Table 5.2 174 

5.10 Example case HI: (a) schematic of the circular shaft-bearing-plate-mount 

system and (b) A SEA model of the shaft-bearing-plate-mount system 176 

5.11 Boundary conditions for example case III: finite shaft-bearing system. 

(a) Moment applied at the bearing end. (b) Force applied at the free end.... 179 


xv 



5.12 Comparison between theory and experiment [17] for example case HI 

2 

with a finite shaft. Here plate mobility level in dB is 101ogio<V c > re 
<V^>=1.0(m/s) 2 for F ysa =1.0N 183 

5.13 Comparison between theory and experiment [17] for example case HI. 

Here theory considers a semi-infinite shaft 184 

5.14 Predicted coupling loss factors T| sc for a semi-infinite and a finite shaft 

in example case HI 186 

5.15 Effect of bearing preload on the normalized casing plate response 

(example case DI) 187 

5.16 Example case IV: (a) a generic geared rotor system with casing and mounts 

and (b) boundary conditions for a finite shaft-bearing system 189 

5.17 Comparison between theory and experiment for example case IV 
(Y=10, Q=2, cr c =c C 2 , §zm = 0 04mm). (a) Sound pressure level Lp. 

(b) Spatially averaged casing acceleration 195 

5.18 Effect of Y^n s /r)c on predicted Lp for example case IV 

(H c =0.02, Q=2, cr c =a cl , 8 zm =0.04mm) 197 

5.19 Sound pressure level prediction with two O c formulations for example case IV : 
a cl = 1.0 (ideal radiator) and a c 2 for simply-supported plate 

(y=10, Q=2, t| c =0.02, 8 zm =0.04rnrn) 198 

5 .20 Effect of bearing axial preload 8^ on Lp for example case IV (rj c =0.02, 

Q=2, 7=10) 199 


XVI 


LIST OF SYMBOLS 


A 0 unloaded distance between the inner and outer raceway groove curvature 
centers 

Aj loaded distance between the inner and outer raceway groove curvature centers 
ofj-thball 

A(to) accelerance transfer function 

<A 2 > mean square acceleration spectra 

Ac casing plate surface area 

Ag shaft cross sectional area 

a 0 ,aj locations of the outer and inner raceway groove curvature centers 

[C] system damping matrix 

[C]b bearing damping matrix 

c wave speed 

DOF degree of freedom 

dg gear diameter 

dp pinion diameter 

djjj inner raceway diameter of bearing 


xvu 


bearing pitch diameter 


^bm 

djj 0 outer raceway diameter of bearing 

d s shaft diameter 

El flexural rigidity of the shaft 

Ec Rayleigh dissipation function 

Ej kinetic energy 

E U potential energy 

E s ,Ec total vibratory energy level in the internal (s) subsystem and external (c) 
subsystem (Chapter V) 

e(t) static transmission error 

F w generalized force in Lagrangian formulation 

Fjb a (t) alternating bearing force in the j = x, y, z or r direction 

Fjb m mean bearing force in the j = x, y or z direction 

Fjsa(t) applied alternating force on the shaft, j = x, y, z or r 

Fjs a p complex Fourier coefficient of Fjsa^)» P - 1* 2, 3, .... 

Fj sm applied mean force on the shaft, j = x, y or z 

{ f(t) } a generalized alternating applied load vector 

{ f } ap complex Fourier coefficient of { f(t) } a , p = 1, 2, 3, .... 

{ f } bm mean bearing load vector 


xvui 


{ f(t) } ca alternating casing load vector 


{ f(t) } s total shaft load vector 
{ f(t) } sa alternating shaft load vector 
{ f } sm mean shaft load vector 
G ' bearing outer ring geometrical center 

H k nonlinear functions, k = 1,2,3,...,V, given by equations (2.12a,b) and (2.13) 
h geometrical thickness 

IjG mass moment of inertia of the gear about j = x, y, z axis 

IzL mass moment of inertia of the load about z axis 

IzM mass moment of inertia of the motor about z axis 

Ijp mass moment of inertia of the pinion about j = x, y, z axis 

Ij c mass moment of inertia of the casing about j = x, y, z axis 

Ij s mass moment of inertia of the shaft about j = x, y, z axis 

{ I } Lumped inertia row vector 

K n rolling element load-deflection stiffness constant 

[K] system stiffness matrix 

[K] v mount stiffness matrix 

[KJbm proposed bearing stiffness matrix of dimension 6 


xix 


[K]bms a matrix of dimension 5 as a subset of [K] bm with last row and column 
excluded 

[K] b gear mesh coupling stiffness matrix 

[K] s shaft stiffness matrix 

[K] e shaft segment stiffiiess matrix 

s 

kjM effective torsional stiffiiess of the driving shaft coupling 
k^L effective torsional stiffiiess of the driven shaft coupling 
k biw bearing stiffiiess coefficient, i, w = x, y, z, 0 X , 0 y , 0 Z 

k biw estimated bearing stiffness coefficient, i, w = x, y, z, 0 X , 0 y , 0 Z 

k s ,k c shaft (s) or plate (c) wavenumber (Chapter V) 

kj, gear mesh stiffness 

k v j mount stiffiiess coefficient, j = x, y, z, 0 X , 0 y , 0 Z 

L a frequency bandwidth-averaged mean square acceleration in dB re <A 2 > re pl.0g 2 

Ly spatially and frequency bandwidth averaged mean square mobility level (dB re 
<V 2 > re f=1.0m 2 /N 2 s 2 ) 

Lg shaft segment length 

L s shaft length 

Mj ba alternating bearing moment about j = x or y direction 
Mjbm mean bearing moment about j = x or y direction 


xx 


Mjsa alternating shaft bending moment about j = x or y direction 

[M] system mass matrix 

[M] s shaft mass matrix 

[M] c casing mass matrix 

m^ rotor mass 


rr^ casing mass 

mg shaft mass 

rrigej shaft lumped mass 

( m } lumped mass row vector 

N g number of gear teeth 

n rolling element load-deflection exponent 

casing plate modal density (see Chapter V) 
n s number of shaft segments (Chapter IV) or shaft modal density (Chapter V) 
p(t) normalized dynamic transmission error 

Qj resultant normal load on the j-th rolling element 
Q acoustic source directivity (Chapter V) 

q w generalized displacement vector in Lagrangian formulation 
{ q(t) } a generalized alternating displacement vector 
{q)ap complex Fourier coefficient of {q(t)) a , p = 1, 2, 3, ... 


xxi 



{ q } b total bearing displacement vector 

{ q } b a alternating bearing displacement vector 

{qlbm mean bearing displacement vector 
{ q(t) } ca alternating casing displacement vector 
{q(t)} sa alternating shaft displacement vector 
{q} sm mean shaft displacement vector 

R(to) force or moment transmissibility transfer function or room constant (see 
Chapter V) 

Re real part of a complex number 
— > 

R j j-th bearing coordinate vector, j = 1, n s +l, n s +2, 2n s +2 
tl bearing radial clearance 

rj pitch radius (roller) or radii of inner raceway groove curvature centers (ball) 

S room surface area (Chapter V) 

S c total surface area of the casing plate in example case m (Chapter HI) 

[T]b bearing field transfer matrix 

Tj sa (t) applied alternating torque on the shaft, j = x, y or z 
Tjsm applied mean torque on the shaft, j = x, y or z 
t time 

Uj ca (t) alternating casing translational displacement in the j = x,y or z direction 


xxu 



uj sa (t) alternating shaft translational displacement in the j = x,y or z direction 
u jsm mean shaft translational displacement in the j = x, y or z direction 
{ u(t) } ca casing alternating displacement vector 
{ u (t) } sa shaft alternating displacement vector 
V(oo) mobility transfer function 

<V 2 > spatially and frequency bandwidth averaged mean square mobility 
v jsa’ v jca shaft (s) or plate(c) alternating velocity in the j=x,y,z direction 
X solution to the problem H k = 0, k = 1, 2, 3, V 

Z total number of rolling element 

Z s driving point shaft impedance (Chapter V) 

Z c driving point plate impedance (Chapter V) 
z number of loaded rolling element 

z D characteristic impedance of the surrounding medium (Chapter V) 
a c unloaded bearing contact angle 

otj loaded j-th rolling element contact angle 

a average absorption coefficient (Chapter V) 

|3j a alternating bearing angular displacement about the j = x or y direction 
Pjm mean bearing rotational displacement about the j = x or y direction 
AFrb m incremental mean radial bearing force 


XXUI 


ASru, incremental mean radial bearing displacement 

A© frequency bandwidth 

Sgj resultant elastic deformation of the j-th ball element 

5r j resultant elastic deformation of the j-th roller element 

5X error vector 

Sjm mean bearing translational displacement in the j = x, y or z direction 

5j resultant elastic deformation of the j-th rolling element 

^ja alternating bearing translational displacement in the j = x, y or z direction 

(5) r j effective j-th rolling element displacement in the radial direction 

(5) z j effective j-th rolling element displacement in the axial direction 

d>j stability functions given by equation (3.6) and Table 3.1 

{ (j) ) j mode shape vector Corresponding to the j-th natural frequency 

r| s ,q c dissipation loss factor for the internal (s) or external (c) subsystems 

T| sc coupling loss factor 

K radius of gyration 

{ fi } p complex Fourier coefficient of the principal coordinate vector, p = 1 , 2, 3, ... 
n s external power input to the internal (s) subsystem 

Fisc net power transfer from the internal (s) subsystem to external (c) subsystems 
p material density 
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0j ca (t) alternating casing angular displacement about the j = x, y or z coordinate 

0j sa (t) alternating shaft angular displacement about the j = x, y or z coordinate 

0jsm mean shaft angular displacement about the i-axis; j = x, y, z 

1 6(0 ) Ra rotor alternating angular displacement vector 
{ 0(t) } ca casing alternating angular displacement vector 
{ 0(0) sa shaft alternating angular displacement vector 
a Rayleigh damping matrix proportionality constant 
<J C casing plate radiation efficiency (Chapter V) 

£2 z mean rotational speed of the shaft 

£2 zm mean rotational speed of the driving shaft 

£2 zl mean rotational speed of the driven shaft 

to bandwidth center frequency (Chapter V) 

co h gear mesh frequency 

coj undamped natural frequency, j = 1, 2, 3, ... 

coj d damped natural frequency, j = 1, 2, 3, ... 

to 0 fundamental frequency 

(Op excitation frequency, p = 1, 2 , 3, ... 

\|/ / bearing load angle 

\|/j angular distance of j-th rolling element from the x-axis 
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modal damping ratio 
[ ] transpose of a matrix or vector 
I I magnitude or determinant 
( ) first time derivative 

( ) second time derivative 

( ) complex conjugate 

( A ) estimation based on simple models 
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ABSTRACT 


A new mathematical model is proposed to examine the vibration transmission 
through rolling element bearings in geared rotor systems. Current bearing models, 
based on either ideal boundary conditions for the shaft or purely translational stiffness 
element description, cannot explain how the vibratory motion may be transmitted from 
the rotating shaft to the casing. For example, a vibration model based upon the simple 
bearing formulations can only predict purely in-plane type motion on the flexible casing 
plate given only bending motion on the shaft. However, experimental results have 
shown that the casing plate motion is primarily flexural. This study clarifies this issue 
qualitatively and quantitatively by developing a comprehensive bearing stiffness matrix 
of dimension 6 model for the precision rolling element bearings from basic principles. 
The proposed bearing stiffness matrix is verified partially using available analytical and 
experimental data, and is completely characterized. 

This study extends the proposed bearing formulation to analyze the overall geared 
rotor system dynamics including casing and mounts. The bearing stiffness matrix is 
included in discrete system models using lumped parameter and/or dynamic finite 
element techniques. Eigensolution and forced harmonic response due to rotating mass 
unbalance or kinematic transmission error excitation for the following examples are 
computed: I. single-stage rotor system with flexible shaft supported by two bearings on 
rigid casing and flexible mounts, II. spur gear pair system with motor and load inertias 
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attached to two flexible shafts supported by four bearings on flexibly mounted rigid 
casing, and III. case II with flexible casing and rigid mounts. In several of these 
examples, analytical predictions match with measured data which validate the proposed 
formulation. Numerical predictions show that the proposed theory is capable of 
predicting bearing moment transmissibility in addition to force transmissibility. 

A statistical energy analysis model combined with the bearing stiffness matrix is 
developed to predict the high frequency asymptotic dynamic behavior. The mean-square 
spatially averaged vibroacoustic responses are calculated for several example cases. A 
physical gearbox is also analyzed to demonstrate the salient features of the proposed 
technique. Good agreement is found between theory and experiment. 
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CHAPTER I 


INTRODUCTION 


1.1 INTRODUCTION 

Noise and vibration generated by the rotating mechanical equipment including 
geared drives have always been a problem in the implementation of new technology in 
automobiles, rotorcrafts and industrial machines [1-6]. Recently, the need for reliable 
vibration/noise prediction methods have been found to be crucial as faster and lighter 
machines are being designed [6-8]. In most of these rotating systems, structure-borne 
noise paths through bearings, which support the rotating shafts on flexible or rigid 
casings, are dominant [5,6,9]. Hence, in order to obtain reliable mathematical 
prediction of the overall dynamic system, a complete understanding of the vibration 
transmission mechanism through bearings, and the role of bearings as a dynamic 
coupler between the shaft and casing, is critical. 

Current bearing models, based on ideal boundary condition or purely translational 
stiffness element description, cannot explain how the vibratory motion may be 
transmitted from the rotating shaft to the flexible casing and other connecting structures 
in rotating mechanical equipment [10-15]. These simple models are only adequate for 
the free and forced vibration analyses of the rotor dynamic system enclosed in a rigid 
casing. For example, a vibrational model of a rotating system based upon the existing 
bearing models can only predict purely in-plane type motion on the flexible casing plate, 
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given only the bending motion on the shaft. However, experimental results have shown 
that the casing plate motion is primarily flexural or out-of-plane type [9,16,17]. This 
paradox is essentially due to an incomplete understanding of the bearing as vibratory 
motion transmitter in rotating mechanical equipment 

The main focus of this dissertation is to clarify this issue quantitatively and 
qualitatively by developing a new mathematical model for the precision rolling element 
bearings, and extends the proposed bearing formulation to examine vibration 
transmissibility in rotating mechanical equipment through several example cases of 
bearing systems and geared drives. The superiority of the proposed model compared to 
simple models is also demonstrated in these example cases. A typical shaft-bearing- 
casing-mount system is shown in Figure 1.1. The rigid or flexible shaft may be 
subjected to forces and/or torques and supported by a bearing on a flexibly or rigidly 
mounted casing. Here, the vibration transmission is from the shaft to the casing and 
mount through the bearing system. Figure 1.2 shows a typical rolling element bearing 
subjected to forces and moments due to the rolling element deformation. The bearing is 
free io rotate about the axis perpendicular to the bearing plane, and hence does not 
transmit any dynamic moment about this axis. However, dynamic moments about the 
other two orthogonal axis exist which have not been considered in simple bearing 
models. Finally, a generic geared rotor system consists of a motor, spur gear pair, 
flexible shafts, load, flexible couplings, rolling element bearings, flexible casing and 
mounts as shown in Figure 1.3a is also considered. The system is excited by the gear 
kinematic transmission error at gear mesh frequency. A single-stage rotor system with 
rotating mass unbalance excitation as shown in Figure 1.3b is treated as a special case of 
Figure 1.3a. Further description of each system will be presented in later chapters. 
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Figure 1.1 Schematic of a typical shaft-bearing-plate-mount system 
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Figure 1.3 Generic rotating mechanical systems, (a) A geared rotor system, (b) A 
single-stage rotor system 
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1.2 LITERATURE REVIEW 

Simple bearing models either assume ideal boundary conditions for the shaft or 
time-invariant translational springs in the axial and radial directions [10-15]. The ideal 
boundary conditions for the shaft are typically simply-supported for short bearings, 
clamped for long bearings or free (for torsional motion only) [10-12]. Formulas for the 
stiffness coefficients derived from the force-displacement relation commonly used by 
bearing designers [18,19] are given by Harris [18], Gargiulo [14] and White [20]. In 
1982 Rajab [21,22] realized the limitations of the simple models and philosophically 
proposed two additional stiffness coefficients which couple the radial and rotational 
bearing degrees of freedom, given radial and moment about the axis transverse to the 
radial line of action. However, in Chapter II, this study has shown that his analyses of 
the ball and roller bearings were erroneous which resulted in an incorrect bearing 
stiffness matrix of dimension 2. In 1988, Young [23] extended Rajab's [21] analyses to 
include mean axial force while retaining other features of Rajab's model. This resulted 
in bearing stiffness matrix of dimension 3 which this study has again found to be 
incorrect 

Simple bearing models are widely used in vibration models of rotor dynamic 
systems which typically exclude casing and mount dynamics, to calculate critical speeds, 
dynamic stability, and responses due to shaft excitations such as mass unbalance and 
gear transmission error [10-15]. In most cases, the vibration transmission through 
bearings is not the primary issue, and moreover such vibration models are clearly 
inadequate in- this regard. A satisfactory study of the vibration transmission through 
bearings in systems similar to Figure 1.1 has yet to be reported even though several 
researchers have considered it with limited success [20-22,24]. 
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Publications on the modal analyses of geared drives and single or multiple stage 
rotor systems indicate that very little has been reported directly on the force 
transmissibility through bearings, and the dynamic effects of bearing, casing and 
mounts on the internal rotating system [10-13,15,25-29]. A comprehensive review of 
the available vibration studies of casing and mounts has been given by Lim and Singh in 
1989 [6]. In these studies, the dynamic interaction between the casing-mount system 
and the internal rotating system is often not incorporated, and in few instances when this 
interaction is modeled, only purely radial and/or axial force on the bearings are included. 
Such models are still inadequate in explaining how the vibration is transmitted from the 
shaft to the casing. A more comprehensive review of the relevant studies will be given 
in each chapter. 

1.3 SCOPE AND OBJECTIVES 

A new mathematical model for the precision rolling element bearing in Figure 1.2 
is developed and incorporated in linear time-invariant discrete and broad band vibration 
models of Figures 1.1 and 1.3. This study proposes a comprehensive bearing stiffness 
matrix of dimension 6 which explains the vibratory motion transmission through the 
bearings and allows for the study of overall geared rotor system dynamics. The lumped 
parameter and dynamic finite element techniques are used to develop the discrete 
vibration models while statistical energy analysis method is used for the broad band 
vibration models. Experimental validation is also included in each chapter, where the 
driving point and cross point accelerance and mobility levels predicted by theory are 
compared to experiments. 
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The specific objectives of this dissertation are grouped as follows: bearing 
stiffness formulation, bearing system studies, geared rotor system studies, and statistical 
energy analysis. Each chapter is self sufficient since it is written in a journal paper style. 
Accordingly, a detailed problem statement is also included in each chapter. 

a. Bearing stiffness formulation : Since simple bearing models used in rotor dynamic 
analyses are inadequate in explaining the role of bearing as a vibratory motion 
transmitter, this study resolves this issue by proposing and developing a new rolling 
element bearing stiffness matrix which is suitable for the analysis of the vibration 
transmission through either ball or roller bearing. A numerical scheme is also developed 
to compute the stiffness coefficients and indicate the existence of solutions to the 
nonlinear algebraic bearing equations describing the bearing load-displacement 
relationships. The proposed bearing model is partially verified by comparing with 
published analytical and experimental results. In addition, the character of the bearing 
stiffness matrix and its sensitivity to various bearing parameters will be discussed. 
(Chapter II) 

b. Bearing system studies : The specific objectives of this chapter are to incorporate the 
proposed bearing matrix developed in Chapter II in linear discrete vibration models of 
the bearing system as shown in Figure 1.1 using lumped parameter and dynamic finite 
element methods to compute the eigensolution and forced harmonic response, and to 
evaluate the dynamic stability. The vibration transmission through bearing is also 
predicted for several example cases considered previously [14,20,24] and an 
experimental setup [17]. The advantages of the proposed formulation compared to the 
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simple models is demonstrated by comparing their predicted transfer functions. The 
theory is also validated by comparing analytical predictions with experimental data on a 
shaft-bearing-plate system. (Chapter HI) 

c. Geared rotor system studies : Current geared rotor system vibration models exclude 
the effects of casing and mounts, and do not address the overall system behavior. The 
objectives of this chapter are to incorporate the proposed bearing matrix in the discrete 
vibration model of the generic geared rotor system of Figure 1.3 and conduct overall 
system studies by calculating eigensolutions and forced harmonic responses with 
emphasis on the prediction of vibration transmission through rolling element bearings. 
The effects of casing and mount dynamics on the internal rotating system is also 
evaluated. Example cases which include a single-stage rotor system with flexible shafts 
supported by two identical rolling element bearings on rigid casing and flexible 
mounts.and a spur gear pair with motor and load inertias attached to two flexible shafts 
supported by four rolling element bearings on rigid or flexible casing and compliant or 
massive mounts will be studied analytically and/or experimentally. Also, the advantages 
of the proposed formulation as compared to simple models of geared drives will be 
demonstrated. (Chapter IV) 

d. Statistical energy analysis : At very high frequencies, the narrow band approach using 
the lumped parameter or dynamic finite element model may not be adequate due to the 
high structural modal density. To overcome this problem, statistical energy analysis 
method is used to predict the vibratory energy transmission in and noise radiation from a 
geared rotor system as illustrated in Figure 1.4a. The proposed bearing matrix is again 
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Figure 1.4 Vibration energy transmission in a geared rotor system, (a) Structure-borne 
noise paths, (b) Statistical energy analysis model. 
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incorporated in the vibratory energy model. In this method, only the mean-square 
spatially averaged response over one third octave frequency bandwidths are predicted 
which involves solution to a set of algebraic equations obtained through the vibratory 
energy balance of each subsystem shown in Figure 1.4b. Several example cases 
including a plate-cantilevered rectangular beam, circular shaft-bearing-plate system and a 
geared rotor system are chosen to demonstrate the salient features of this technique. 
(Chapter V) 



CHAPTER 2 


BEARING STIFFNESS FORMULATION 
2.1 INTRODUCTION 

Current rotor dynamic models describe precision rolling element bearings either as 
ideal boundary conditions for the shafts [10-12], or as purely translational stiffness 
elements [13-15]. Such simple bearing models may be adequate for the free and forced 
vibration analyses of the rotor dynamic system enclosed in a rigid casing. But these 
mathematical models cannot explain how the vibratory motion may be transmitted from 
the rotating shaft to the flexible or rigid casing and other connecting structures. For 
example, a vibration model of a system similar to Figure 2.1, based upon the existing 
bearing models, can only predict purely in-plane type motion on the flexible casing plate 
given only the bending motion on the shaft. However, experimental results have shown 
that the casing plate motion is primarily flexural or out-of-plane type [9,16,17]. This 
paradox is essentially due to an incomplete understanding of the bearing as vibratory 
motion transmitter in rotating mechanical equipment including geared drives where 
structure-borne noise paths through bearings are often dominant. 

This chapter clarifies this issue qualitatively and quantitatively by developing a 
new mathematical model for precision rolling element bearings. A schematic of a 
generic system with a flexible shaft rotating at speed Cl z and subjected to mean load 
vector { f } sm = ( ^ wsm ,T wsm } , w — x, y, z, flexible casing and mount is shown in Figure 
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x 



Figure 2.1 Schematic representation of the vibration transmission problem. Here the 
flexible shaft is subjected to mean forces F wsm and torques T wsm where w 
— x, y or z, is the direction and subscript m and s implies mean and shaft 

respectively. Also, 0 is the angular displacement and u is the translational 
displacement. 
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2.1; the shaft is supported on one of the following bearings: deep groove ball bearing, 
angular contact ball bearing, thrust ball bearing, straight roller bearing or taper roller 
bearing. A new bearing stiffness matrix [K] bm will be proposed which is expected to 
demonstrate a coupling between the shaft bending motion and the flexural motion of the 
casing plate. It will be shown that the translational bearing stiffness coefficients 
currently used in rotor dynamic models are a small subset of the proposed [K] bm . 
Several example cases are employed to validate our theory. Our bearing model can be 
easily incorporated in analytical or numerical models typically used for the dynamic 
analyses - this will be the basis of Chapters ID and IV of this dissertation. 

2.2 LITERATURE REVIEW 

The ideal boundary conditions for the shaft have typically been assumed to be 
simply-supported for short bearings, clamped for long bearings or free (in the torsional 
mode only) [10-12]. In other cases, researchers describe the bearing as time-invariant 
translational springs with stiffness coefficients k brT and/or k bzz in the radial and axial 
directions, respectively [13-15,20]. Formulas for such nonlinear stiffness coefficients 
are given by Harris [18] and Gargiulo [14]; these are derived from the radial or axial 
mean force-displacement equation commonly used by the precision rolling element 
bearing designers [18,19]. Their derivations neglect the effects of radial clearance and 
mean bearing force vector {f} bm on the load distribution and hence are applicable only 
for constant load angle v/i of 180 degrees. White refined these formulations by using a 
finite difference approximation for the computation of stiffness coefficients for radial 
ball and roller bearings, and by including the effects of radial clearance and force on the 
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load angle Y/ [20]. Even with these refinements the mathematical model is still 
inadequate and incapable of predicting the vibration transmission across bearings. 

In 1982 Rajab [21], realized the limitations of the current simple theory and 
philosophically proposed two additional stiffness terms k br Q and k^Q which couple the 
relative radial and rotational bearing displacements between the inner and outer rings, 
given the mean radial load and moment about the axis transverse to the radial line of 
action. However our investigation, as demonstrated later in Sections 2.4 and 2.5, has 
shown that his analyses of the ball and taper roller bearings were erroneous which 
resulted in an incorrect 2x2 bearing stiffness matrix. In 1988 Young [23] extended 
Rajab's [21] analyses to include the mean axial force F zbm , and then used a discrete 
summation over all of the loaded rolling elements to obtain bearing forces and moment 
instead of the integral form while still retaining other features of Rajab's model. This 
resulted in a 3x3 bearing stiffness matrix which our investigation has again found to be 
incorrect. Some of the salient features of Rajab's [21] and Young's [23] models are 
summarized in Reference [22]. 

Experimental determination of the bearing stiffness coefficients has been strictly 
limited to the translational coefficients k brr and k bzz . A method for the measurement of 
in situ bearing stiffness under oscillating loading conditions has been given by Walford 
and Stone [30]. Recently, Kraus et al. [24] designed an in situ measurement test stand 
to determine the translational bearing stiffness from measured vibration spectra, in 
conjunction with the single degree of freedom system theory. They determined the 
effect of preload, bearing release and rotational speed £2 z on k bnr and k bzz . Their results 
show that k brr and k bzz are essentially linear and the effect of £2 Z is negligible when a 
high preload is applied on the bearing. 
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2.3 ASSUMPTIONS AND OBJECTIVES 

Due to the following key differences, a separate formulation of [K]b m for both ball 
and roller type rolling element bearings is required: (i) ball bearings have elliptical 
contacts and roller types have rectangular contacts between the inner race, rolling 
elements and outer race when loaded, and (ii) the loaded contact angles ocj of the ball 
types may change but (Xj in the roller type remains relatively constant [31]. Each bearing 

is characterized by its kinematic and design parameters such as unloaded contact angle 
a Q , radial clearance r L , effective stiffness coefficient K n for inner ring-single rolling 

element-outer ring contacts, angular misalignment, preloads, radius of inner raceway 
groove curvature center for ball type and bearing pitch radius for roller type [18,19,31], 
It is expected that [K] bm is given in terms of these parameters. 

The mean bearing displacements [q] bm as shown in Figure 2.2 are given by the 
relative rigid body motions between the inner and outer rings. The total bearing 
displacement vector is given as {q]b={q}bm + { £ l( t )}ba where {q(t)} ba is the fluctuation 
about the mean point {q]b m during the steady state rotation. Accordingly one must 
consider time varying bearing stiffness coefficients. However in our analysis, such time 
varying bearing stiffness coefficients are neglected by assuming very small vibratory 
motions i.e. fq} ba « {q)bm> and high bearing preloads. Consequently, only the mean 
bearing loads and displacements are included in the derivation of [K]b m . The basic 
load-deflection relation for each elastic rolling element is defined by the Hertzian contact 
stress theory [18,19,32], and the load experienced by each rolling element is described 
by its relative location in the bearing raceway. Further it is assumed that the angular 
position of each rolling element relative to one another is always maintained due to the 
rigid cages and pin retainers. Secondary effects such as centrifugal forces and 
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Figure 2.2 Rolling element bearing kinematics and coordinate system. Here the 
following nomenclature is used: d^ is the outer raceway diameter, d bm is 
the bearing pitch diameter, d b j is the inner raceway diameter, \j/ is the 
angular position of rolling element, 8 wm is the mean translational 
displacement, p pm is the mean angular displacement, F wbm is the mean 
bearing force, and M pbm is the mean bearing moment where w = x, y, z, 
and p=x, y, are the directions. 
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gyroscopic moments on the bearing are ignored as these effects are evident only at 
extremely high rotational speeds. Tribological issues [32,33] are beyond the scope of 
this study and hence our analysis assumes bearings to be unlubricated. 

The specific objectives of this chapter are to: (i) propose and develop a new rolling 
element bearing stiffness matrix [K] bm which is suitable for the analysis of the vibration 
transmission through either ball or roller bearing, (ii) develop a numerical scheme to 
compute [K] bm and discuss the existence of solutions to the nonlinear algebraic bearing 
equations describing load-displacement relationships, (iii) verify our proposed model by 
comparing its predictions with published analytical and experimental results [14,20,24] 
for the translational stiffness coefficients k bxx , k byy and k bzz , (iv) relate [K] bm to 
various kinematic and design parameters, and perform parametric studies to investigate 
the effect of unloaded contact angle a Q and preloads, and (v) characterize the nature of 
[K] bm and recommend its usage. Finally it should be noted that dimensionless 
parameters will not be used here as the metric units are invariably employed to specify 
bearings [32]. 

2.4 BEARING LOAD-DISPLACEMENT RELATIONS 

In this section, the relationships between the bearing forces (F xbm , F ybm , F zbm ) 
and moments {M xbm , Myb m } transmitted through the rolling element bearing, and the 
bearing displacements {q)bm as given in Figure 2.2 will be derived for both ball and 
roller bearings. The mean applied loads [f} sm at the shaft as given in Figure 2.1 and 
bearing preloads generate the mean bearing displacements [q} bm and loads [f} bm . 
These displacements {q)bm are use£ l to derive the resultant elastic deformation 5(\}/j) of 
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the j-th rolling element located at angle \j/j from the x-axis. From the ball bearing 
kinematics shown in Figure 2.3, S B (\j/j) is 


A( V )-A 0 . 8 Bj >0 
8 B ( Vj )= |o , 8 Bj S0 

A<Vj) - V< 5 '^ j - < S * 

(S*) zj = A o sin «o+( 5 ) zj ; (S*) rj = A 0 cos a 0 + (8) rj 


(2.1a) 


(2.1b) 


(2.1c) 


where A 0 and A are the unloaded and loaded relative distances between the inner and 
outer a^, raceway groove curvature centers. Similarly for the roller bearing kinematics 
shown in Figure 2.4 for (Xj=a 0 , 8 R (\|/j) is 



(8) . cos a. + (8) . sin a. 
n j zj j 

0 


5 Rj>° 

S Ri S0 


( 2 . 2 ) 


Note that in equations (2.1) and (2.2) 8 B j S 0 or 6 r j S 0 implies that the j-th rolling 
element is stress free. In both equations (2.1) and (2.2), the effective j-th rolling 
element displacements in the axial (8) Z j and radial (8) r j directions are given in Figure 2.5 
in terms of the bearing displacements {q)bm* 

(8) . = 8 + r . { P sin (\|/ .) - (3 cos(y .) ) (2.3a) 

zj zm j ir xm VT y 'ym VY j' ' v ' 

(8) rj = S xm cos Vj + 8 ym sin y . - r L 


(2.3b) 
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Figure 2.3 Elastic deformation of rolling element for non-constant contact angle a 
given by the change in the distance between the inner aj and outer a 0 
raceway groove radius curvature centers due to the mean bearing loads or 
displacements. 
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Figure 2.4 Elastic deformation of rolling element for constant contact angle ccj = a D 
given by the change in the relative position of the inner and outer raceways 
due to the mean bearing loads or displacements. 
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where rj is the radial distance of the inner raceway groove curvature center for the ball 
type or is the pitch bearing radius for roller type. Equations (2.1)-(2.3) in conjunction 
with the Hertzian contact stress principle [18,19,32] stated as follows yield the load- 
deflection relationships for a single rolling element . 

Qj = K n Sj (2.4) 


where Qj is the resultant normal load on the rolling element, and K n is the effective 
stiffness constant for the inner race-rolling element-outer race contacts and it is a 
function of the bearing geometry and material properties [18,19,31]. Note that the 
exponent n is equal to 3/2 for ball type with elliptical contacts and 10/9 for roller type 
with rectangular contacts. Previously, we have mentioned that the loaded contact angle 
ocj for the roller bearing remains unchanged from the unloaded position a G , but on the 
other hand (Xj may alter in the ball bearing case. The sign convention is such that otj is 
positive when measured from the bearing x-y plane towards the axial z-axis as shown in 
Figures 2.3 and 2.4, and negative otherwise. For the ball bearing of Figure 2.3, the 
loaded contact angle ctj is 


tan(<Xj) 


A 0 sin a Q + (8) zj . 
A 0 cos a Q + (8) rj 


(2.5) 


where (8) Z j and (8) r j are given by equations (2.3a) and (2.3b). It is appropriate here to 

note that Rajab [21,22] and Young [22,23] in their derivation of the bearing stiffness 
model used an expression similar to equation (2.2) but with 5 xm =py m =8 zm =r L =0 in 
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Rajab' analysis and 8 xm =P ym =r L =0 in Young's analysis for both ball and roller 
bearings. Since always (Xj is given by equation (2.5) irrespective of the formulation and 
since equation (2.2) is valid only if (Xj=a 0 , their ball and roller bearings analyses are 

obviously incorrect. Expressions similar to equation (2.1)-(2.5) with minor differences 
have also been used by Eschmann et al. [31], Jones [34] and Davis [35], but their 
intentions were to calculate static bearing forces rather than to derive the bearing 
stiffness models for vibration transmission analysis. 


2.5 DEVELOPMENT OF BEARING STIFFNESS MATRIX [K] bm 

Our proposed bearing stiffness matrix [K] bm is a global representation of the 
bearing kinematic and elastic characteristics as it combines the effects of z number of 
loaded rolling element stiffnesses in parallel given by 8j > 0. First, we need to relate the 
resultant bearing mean load vector {f} bm to * e bearing displacement vector {q} bm . 
This can be achieved through vectorial sums Qj (5 wm ,P pm ; w = x, y, z and p = x, y) in 
equation (2.4) for all of the loaded rolling elements which lead to the following bearing 
moments [M w b m ] and forces [F w b m ] as follows 


f 
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Replacing Qj and <Xj in equation (2.6) in terms of {S wm ,P pm } yields the following 
explicit relationships between {fJbm and (q}bm f° r ball bearings 


M 
. M 
M 


xbra 

ybm 


zbm 


] * ■[ _y[ A ° sin a° + ( 5> Zj ] 2 + [ A 0 cos a 0 + (8) ^l 2 - A 0 } 
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(2.7a) 


f F xbm 
{ F ybm 
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[A 0 cos a 0 + (5) r .] cos \jr ' 
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[A q sin a 0 + (8) ] 


(2.7b) 


and similarly for roller bearings 


[ M xbm] 

z 

sin v ) 

j M ybm 

- = K n Sin a 0 X r j { (8) ^ cos oc 0 + (8) zj sina 0 } n - 

1 

- cosy . j 

Kbm. 

j 

0 J 


f F xbm 

z 

' cosa 0 cos\jr^ ' 

j F ybm 

’ = K „ X {( 5 ) rj cos a o + < 5 ) zj sin a o)" • 

cos a sin v • • 

O Y J 

l F zbm . 

j 

sin a_ 

L J 


(2.8b) 


26 


where (5) rj - and (5) ZJ - are functions of {8 wm ,P pm } as defined by equation (2.3). 
Approximate integral forms of equations (2.7a, b) and (2.8a,b) are often used instead of 
the summation forms to eliminate explicit dependence on Y), especially in the case of 


only one or two degrees of freedom bearings [18,32]. For instance Rajab [21,22] chose 
the integral form representation but made a mathematical error in constructing the 
integrand which led to further errors in his analysis. 

Now we define a symmetric bearing stiffness matrix [K]b m of dimension 6 from 


equations (2.7a,b) and (2.8a,b) and by assuming that {q)ba « (qlbi 


m 



w, l = x, y, z 


(2.9) 


Here each stiffness coefficient must be evaluated at the mean point (q}bm- Explicit 
expressions for the ball bearing stiffness are as follows; note that [K]b m is symmetric 
i.e. kbi w =kjj W j. 


k bxx K n X- 
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(Aj-A 0 ) cos 2 V j j 


n A. (8*) 2 . 

- A J ~ r rj + a 2 -(S*) 2 

Aj-A 0 J 'tj 

T 3 1 

j 


(2.10a) 
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A 3 
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- 1 
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(2.10c) 


1 

— (2.10d) 


(2.10e) 


(2.1 Of) 


(2.10g) 


(2.1 Oh) 
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k b0 0 K n X" 

y y j 


r j (A j “ A Q ) n cos 2 Vj j 


” A i (S >'i A 2 ,*.,2 

VV +A j- (8 


(2.10o) 


k bi0 = k b0.o =0 ; i = x >y» z 


(2.1 Op) 




where (5 ) z j, (8 ) r j and Aj are defined by equation (2.1). And the roller bearing 
stiffness coefficients k^i w = kb w j are given explicitly as 


k, = nK n cos 2 
bxx n 


k bxy 2 


.n - 1 


J 

z 


Rj 


k b„=| K n sin 2a 0 Z,5" R . ' cos y . 


n - 1 


k byy = nK n COs2<l oS S Rj 

j 


COS 2 

(2.11a) 

1 . - 
sin 2y. 

(2.11b) 

l 

COS \)/j 

(2.11c) 
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sin 2\i/. 
Rj Y ] 

(2.1 Id) 

<> ^ 1 0 
.o n . cos 
j Rj y j 

(2. lie) 
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k byz = 2 K n sin 2a ot «Rj~ * sin y 
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(2.11g) 
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k by0 = 2 K » sin 2 “o? r j 5 Rj ' si " 2 'l' j 


k. = - -7- K n sin 2a n V r. 8" 1 sin 2 w. 
by 0 4 n °“ J Rj Y J 


k b Z z =nK n sin 2 a oZ 5 R i 1 
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k b 0 e =nK n sin2a oS r j 2 8R j 1 cos 2 v j 
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k b,9 = k b0.e =° : i=x * z 

2 | Z 


(2. 1 lh) 

(2. lli) 

(2.11j) 

(2.11k) 

(2.11/) 

(2.11m) 

(2.1 In) 

(2.1 lo) 
(2. lip) 


where 8 R j is defined in equation ( 2 . 2 ). It should be noted that all stiffness terms 
associated with the torsional degree of freedom p^ are zero due to the fact that an ideal 
bearing allows free rotation about the z-direction. Also, the translational stiffness 
coefficients k bii , i=x,y,z for 8 ym =8 zm =p xm =p ym =0, 5 xm = 5 zm =P xm =p ym =0 or 
Sxm = 5 ym = Pxm = Pym = 0 are equivalent to the bearing stiffness coefficients commonly 
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used by investigators [14,15,20]. The nature of these and other features of [K] bm will 
be discussed later in Section 2.9. 

2.6 NUMERICAL ESTIMATION OF [K] bln 

The coefficients k b j w can be computed by one of the following two methods: I. 
directly compute k b j w given mean bearing displacement vector {q)bm employing 
equations (2.10a-p) and (2.11a-p), or II. numerically solve the nonlinear algebraic 
equations described by equations (2.7a, b) and (2.8a, b) to obtain {q}bm from {f} bm , 
and then evaluate k b ; w per method I. Note that {f} bm may be functions of the mean 
shaft loads, bearing preloads, and shaft and casing compliances depending on the 
configuration and flexibility of the rotating mechanical system. If the bearing system is 
statically determinate, then {f Jbm ma y be computed explicitly in terms of {f} sm and 
preloads using the force and moment equilibrium equations. Conversely for an 
indeterminate system, appropriate field equations for the shaft and casing plate are 
needed in addition to the equilibrium equations to obtain {f Jbm which must also include 
shaft and casing compliances. Calculations of [f} bm and {q} bm in this case are 
simultaneous, which may be extensive especially when the system is very flexible, and 
may even require discretization using finite element or lumped mass technique. 
However, in many real machines the in-plane stiffness of the casing plate which 
supports most of the mean bearing load is much higher than the bending stiffness of the 
shaft. Hence the casing in-plane stiffness term may be neglected without contributing 
any large error to [f} bm [18,31]. And only the Euler’s beam equation for a statically 
indeterminate shaft is used along with the nonlinear bearing load-displacement equations 
(2.7a, b) and (2.8a, b). 
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Method I is computationally direct and needs no discussion. But method II deals 
with as many as 10 N nonlinear algebraic equations for N bearings if the casing 
flexibility is neglected. One must choose an appropriate numerical method as the 
nonlinear algebraic equations must be solved iteratively [36,37]. In addition, the 
available numerical methods need a prior knowledge of the approximate location of the 
solution vector being sought and hence one must be careful in interpreting the numerical 
results. In this study, we adopted the Newton-Raphson method for its good 
convergence characteristic [36,37]. To implement this method, equation (2.6) for each 
bearing is rearranged as 



(2.12a) 


(2.12b) 


where H^,H2,...,H5 are functions defined for computational reasons. For an 
indeterminate system, there are additional functions H6,H7,...,Hv from the field 
equations. Using Taylor’s series, any function in equations (2.12a,b) can be 
expanded about the solution vector X = [q]i, m for a statically determinate system andX 
= [ [q] bm , {f } . ] for a statically indeterminate system as follows by neglecting 


second and higher order terms. 
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V 9H. 

H k (X+8X) - H k (X) + 5X. *» k=l,2,3,...V (2.13) 

j J J 

The solution for the incremental vector 8X can be obtained by setting H k (X+8X) = 0 
per equations (2.12) and (2.13) which yields a set of linear algebraic equations. This 
vector 8X is added to the previously computed vector X given by H k (X) = 0 for the 
next iteration until the convergence criterion, say that 5X is within a specified tolerance, 
is satisfied. Our proposed numerical scheme can be summarized as follows: (i) guess 
bearing displacement vector {q)bm and/or load vector {f) bm , (ii) compute 8X and check 
against a specified tolerance, (iii) add SX to the previous solution vector X and repeat 
steps (i) and (ii) until the convergence criterion is satisfied. We have found that a few 
initial guess trials are required in most cases to obtain reasonable results. 

2.7 VALIDATION OF PROPOSED MODEL 

In order to validate our theory we compare the translational stiffness coefficients of 
the proposed bearing matrix [K] bm with published analytical and experimental results 
[14,20,24]. First we apply our theory to predict the nonlinear axial k bzz = k bzz (5 zm ) 
and radial k brr = k brr (5 rm ) stiffnesses as shown in Figure 2.6. Our predictions are 
found to be within 2 % of Gargiulo's [14] formulas which are commonly used for both 
ball and roller bearings. 

For the second example case, we consider the ball bearings used by Kraus et 
al.[24] for an in-situ determination of the bearing stiffness. Using their bearing design 
parameters, we compute radial stiffness coefficient k brr as a function of the axial preload 
F zbm . Excellent comparison between theory and experiment is seen in Figure 2.7. 



0.025 

Radial Deflection (mm) 


i 


Figure 2.6 Comparison between the proposed theory and Gargiulo's formulas [14] for 
axial and radial stiffness coefficients of ball and roller bearings. 
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Axial Preload (kg) 


Figure 2.7 Comparison between the proposed theory and the experimental results of 
Kraus et al.[24] for as a function of the mean axial preload. 
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Finally, we compare our results for the nonlinear radial stiffness k brr with those 
reported earlier by White [20] for both ball and roller bearings. We note discrepancies 
in Figure 2.8 between our theory and White's results. In order to explain these we now 
define k brr using the finite difference approximation which was also used by White: 
k brr = AF brm / ASjm ~ F brm / (5 rm -r L ). Now a good match is evident in Figure 2.8 

A 

between our k brr values and the data given by White. However, the correct formulation 
is obviously given by our proposed theory which is based on the analytical partial 
derivatives k brr = dF brm / 35m, as the displacement 8,™ may be large. 

2.8 PARAMETRIC STUDIES 

The proposed matrix [K] bm includes a coupling between the casing flexural 

motion and shaft bending motion which is reflected by some of the dominant off- 
diagonal, k bx0y , k by0x , k bz0x and k bz0y , and rotational diagonal, k b0x0x and k b0y0y , 

stiffness coefficients; these are labeled as 'coupling coefficients’ for discussion 
purposes. Such stiffness coefficients are investigated further by varying preloading 
conditions and unloaded contact angle a ^ for both ball (set A) and roller (set B) bearings 
whose design data are listed in Table 2.1. 

The coupling coefficients given a constant mean radial displacement 8 rm (radial 
preload), as shown in Figures 2.9 and 2.10 for both ball and roller bearings 
respectively, are found to increase as a Q increases and reach a maximum when a Q is 
near 90°. On the other hand, the radial translational stiffness coefficients in the x and y 
directions are found to decrease as <Xq increases. These observations imply that for deep 
groove ball type or straight roller type bearing (a 0 « 0°) the radial stiffness coefficients 
k brT are dominant, but for angular contact ball type or taper roller type bearing (Oq > 0°) 
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Table 2. 1 Design parameters for typical ball and roller bearings used for parametric 
• . studies 


Parameters 

Set A (ball type) 

Set B (roller type) 

Load-deflection exponent n 

3/2 

10/9 

Load-deflection constant K n (N/m n ) 

8.5 E9 

3.0 E8 

Number of rolling element Z 

12 

14 

Radial clearance r L (mm) 

0.00005 

0.00175 

Pitch radius ^(mm) 

19.65 

21.25 

Ao (mm) t 

0.05 

- 


t Unloaded distance between inner and outer raceway groove curvature centers (see 
Figure 2.3) 

ft Equivalent to rj for roller bearings and for ball bearings given in equation 

(2.3) 





Bearing Stiffness (GN/rad) Bearing Stiffness (GNm/rad) Bearing Forces (MN/mm) 
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Unloaded Contact Angle (deg) 


Figure 2.10 Dominant stiffness coefficients of roller bearing set B for 0° < a 0 < 90° 
and given a constant mean radial bearing displacement S xm = 0.025 mm. 
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the coupling terms are more significant. Note that in Figure 2.10, all the stiffness 
coefficients are zero at a 0 = 90° for the roller type. This is due to the fact that in the 
thrust roller bearing, radial flanges are included to resist the roller motion in this 
. direction which is not modeled here, and hence these stiffness coefficients must vanish. 
In addition, thrust roller bearings are designed to carry axial loads [18,31]. On the other 
hand, ball bearings have finite stiffness coefficients at Oq = 90° due to the curvature of 
the raceway which provide some resistance to the radial preloads. In general, the trends 
in both ball and roller bearing stiffness properties are similar when each is subjected to 
mean radial displacement or preload. 

In the case when the bearings are subjected to mean axial displacement (axial 
preload), as shown in Figure 2.11 for the ball type and Figure 2.12 for the roller type, 
the number of nonzero stiffness coefficients are less than those seen for the radial 
preload only. Again, it is observed that both ball and roller bearings display similar 
trends. Over mid to high oc 0 values, the coupling coefficients are found to be 
significant. The translational stiffness coefficients are relatively constant except for the 
axial stiffness which increases as Oq increases. This is expected due to the inclination of 
the rolling element line of contact from the x-y plane which increases elastic support in 
the z-direction. At a Q = 0°, all the stiffness coefficients for roller bearings are zero as 
there is no constraint in the axial direction. In real bearings such a constraint is provided 
by the axial flanges [18,31], however this bearing is not designed to carry any axial 
preload. 

Results for the misalignment in ball and roller bearings simulated by specifying a 
mean bearing angular displacement P ym are shown in Figures 2.13 and 2.14 
respectively. The dominant stiffness coefficients are the same as those seen for the 


Bearing Stiffness (MN/rad) Bearing Stiffness (MNm/rad) Bearing Stiffness (kN/mm) 
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Figure 2.1 1 Dominant stiffness coefficients of ball bearing set A for 0° < a 0 < 90° and 
given a constant mean axial bearing displacement 5 zm = 0.025 mm. 




Bearing Stiffness (GN/rad) Bearing Stiffness (GNm/rad) Bearing Stiffness (MN/mm) 
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Figure 2.12 Dominant stiffness coefficients of roller bearing set B for 0° < a 0 < 90° 
and given a constant mean axial bearing displacement 8 zm = 0.025 mm. 





ing stiffness (MN/rad) Bearing Stiffness (MNm/rad) Bearing Stiffness (MN/mm) 
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Unloaded Contact Angle (deg) 


Figure 2. 1 3 Dominant stiffness coefficients of ball bearing set A for 0° < a 0 < 90° and 
given a constant misalignment p ym = 0.015 rad. 



Bearing Stiffness (GN/rad) Bearing Stiffness (GNm/rad) Bearing Stiffness (MN/mm) 
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radial preload case. For ball bearings, most of the stiffness coefficients remain constant 
for 0° < a c < 90°. On the other hand, the stiffness coefficients for roller bearing have 
trends similar to those found for the radial preload cases. 

From the detailed parametric studies, it is concluded that the nature of [K] bm is 
dictated by the bearing type, a 0 and preloads. Also, the coupling coefficients are not 
negligible in most cases as assumed previously by many investigators. 

2.9 CONCLUSIONS 

Results of Section 2.8, which show similar trends for some of the cases, imply 
that there may be a systematic approach to characterize the proposed bearing stiffness 
matrix [K] bm . From the kinematic and geometrical considerations, it is always possible 
to impose any bearing displacement vector {q)bm which denotes relative rigid body 
motions between the inner and outer rings as long as the rolling element is still within 
the elastic deformation regime. On the other hand, an arbitrary application of {f} bm may 
not produce a singular displacement response from the bearing due to its kinematic and 
geometrical constraints. Hence, we compute [K] bm and {f} bm by systematically 
varying {q} bm . The results of all possible forms of [K] bm are listed in Table 2.2 and 
2.3 for ball and roller bearings respectively. Also included here are the current bearing 
models which are based on the translational spring descriptions; these models do not 
show any coupling. Note that the exact values of the stiffness coefficients are not given 
as these depend on specific parameters; therefore only the dominant k b jj terms are listed 
for all possible bearing load configurations along with the corresponding {q} bm and Oq. 
Also, note that not all combinations of the bearing loads are possible which complicates 
bearing stiffness calculations further, especially for the numerical method H Tables 2.2 
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Table 2.2 Comparison between the proposed and current ball bearing stiffness 
coefficients.(p = x, y; i = x, y but i* p) 


Mean 

Mean bearing displacement 

Dominant stiffness coefficients 

bearing^ 

loads 






a o ®0° 

0°«x o <90° 

a o «90° 

current^ 

proposed^ 

Ppm 

5pm 

- 

- 

k PP 

^xx> k yy’^zz>^0x0x’ k 0y0y ,k z9i 

Pzm 

5zm 

5zm 

— 

^pp»^zz 

^xx>^yy>^zz> k 0x0x ,k 0y0y’ 

^X0y>^y0x 

Pzm 

- 

- 

5zm 

kzz 

^xx’^yy^zz’^OxOx^OyOy 

Mpm 

Ppm 

- 

- 

- 

^xx»^yy»^zz>^0x0x» k 0y0y» k iz 

Pzm'Mpm 

— 

- 

5zm>Ppm 

kzz 

^xx»^yy»^zz»^0x0x> k 0y0y’ k z0p 

Pxm>Pym 

5xm>5ym 

— 

— 

k pp 

^xx»^yy»^zz»^0x0x ,k 0y0y’ 

^xy’^0xOy’^z0x ,k z0y 

F M 
1 pm» lvl pm 

5pm»Ppm 

— 

— 

k pp 

^xx«^yy’^zz’^0x0x’ k 0y0y’ 

^x0x’^y0y’^pz’^z0i 

^xm>^ym 

5pm»Ppm 

— 

— 

— 

^xx’^yy’^zz’^OxOx’^QyOy’ 

^xy»^xz»^yz,^0x0y 

Ppm>Pzm> 

^im 

”pra>®zm' 

Pirn 

5pm>5zm» 

Pirn 

5pm>5zm» 

Pim 

1c Ic 
^pp^zz 

^XX>kyy,k zz ,k0 x 9 x > k 0y0y> 

^X0y’ky0 X ’kp Z> k z 0i 

Pzm>^xm» 

Mym 

5zm>Pxm> 

Pym 

5zm*Pxm» 

Pym 

5zm>Pxm> 

Pym 

^PP’^ZZ 

kxX’kyy»kzz>k0x0 x , k 0y9y, 

kxy>k0 x 0y.k z 0 x ,k z 0 y 


combinations of {q} m 

^PP’^ZZ 

all non- zero except 0 Z terms 


t Ideal boundary condition models used to describe the bearing are not tabulated. 

tt Here the subscript b which implies bearing has been omitted for brevity. 

ttt All terms associated with 9 Z are zero because of the free rotation about the z axis. 
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Table 2.3 Comparison between the proposed and current roller bearing stiffness 
coefficients.(p = x, y; i = x, y but i * p) 


Mean 

bearing^ 

loads 

Mean bearing displacement 

Dominant stiffness coefficients 

a „~0° 

0°<a o <90° 

a o =90° 

current^ 

proposed^ 

Fpm 

^pm 

- 

- 

k PP 

^bpp 

F Z m 

— 

^zm 

— 

kppkzz 

^xx’^yy’^zz’^GxBx’^QyQy’ 

^X0y>^y0x 

Fzm 

- 

— 

^zm 

^zz 

^zz>^0x0x’ k 0y0y 

F M 

r zm’ m pm 

- 

- 

^zm>Ppm 

^zz 

^zz^OxOx^QyOy’^Op 

Fxm^Fym 

^xm»^ym 

- 

- 

k pp 

|r |r V 

K xx» K yy> A xy 

Fpm’F zm 

^im 


Spm>^zm 

Pim 

— 

k PP’ k ZZ 

^xx>^yy’^zz’^0x0x’^6yOy’ 

^X0y’^-y9x’ k pz’^z0i 

F M 
1 zm» lvl xm> 

Mym 

— 

“ 

^zm’Pxim 

Pym 


^zz’^OxOx^OyOy’^xOy’ 

^z9x ,k z0y 

{f)m 

combinations of 


k pp» k zz 

all non-zero except 0 Z terms 


t Ideal boundary condition models used to describe the bearing are not tabulated. 

tt Here the subscript b which implies bearing has been omitted for brevity. 

ttt All terms associated with 0 Z are zero because of the free rotation about the z axis. 
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and 2.3 should provide some insight to the solution of the nonlinear algebraic bearing 
load-deflection equations which requires a priori knowledge of the type of solution 
being sought as outlined earlier. In most practical problems, mean bearing loads are 
typically known. This knowledge can be combined with Table 2.2 or 2.3 to formulate 
the nonlinear load-deflection equations in the simplest form by deleting all of the zero 
displacement terms. 

Tables 2.2 and 2.3 show that the coupling coefficients kt, x0y , k by0x , k bz0x , k bz0y , 
k b0x0x and k b0y0y are found to be dominant in most of the ball bearing cases, and only 

in some of the roller bearing cases. This is essentially due to the curvature of the 
raceway in ball bearing which invariably causes the rolling element to orient itself such 
that 0° < otj < 90° which generates ball loads in the z direction as well. However, in the 
roller bearing case where <Xj = a 0 , the same phenomenon does not occur when a 0 = 0° 
or 90°, and the coupling coefficients are seen only when * 0° or 90°. In fact for the 

0° and 90° unloaded contact angle cases, the stiffness coefficients associated with x and 
y directions and those associated with the z, 0 X , and 0 y directions do not exist 
simultaneously; the former is dominant when a 0 = 0 and the latter prevails when a 0 = 
90°. Another case of interest here is the case when bearing loads are complex as given 
by the last row in Tables 2.2 and 2.3 where all of the bearing stiffness coefficients 
unrelated to the rotational degree of freedom 0 Z exist. Solution to these cases may 
require a large number of iterations. 

In summary, we have developed a comprehensive bearing stiffness matrix from 
the basic principles which includes all possible rigid body degrees of freedom of a 
bearing system. This matrix has been validated partially using several analytical and 
experimental examples. Further validation of [K] bm is not possible as coupling 
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coefficients are never measured [24,30]. Nonetheless, our theory is general in nature 
and is applicable to even those configurations which may be different from the generic 
case shown in Figure 2.1. Further research is required to incorporate tribological issues 
[32,33] in this formulation. However the proposed stiffness matrix in its present form, 
unlike the current models, is clearly capable of explaining the nature of vibration 
transmission through bearings - this is the subject of Chapters III and IV of this 
dissertation which will also include further comparisons between theory and experiment. 


CHAPTER III 


BEARING SYSTEM STUDIES 
3.1 INTRODUCTION 

Current bearing models [10-15] can not explain how the vibratory motion may be 
transmitted from the rotating shaft to the casing and other connecting structures in 
rotating mechanical equipment. For instance, experimental results [9,16,17] have 
shown that casing plate motion for a system similar to Figure 3.1 is primarily flexural or 
out-of-plane type given only the bending motion on the shaft. Using existing vibration 
models, only in-plane type motions on the casing plate are obtained. Such limitations 
associated with current bearing models have been discussed thoroughly in Chapter II of 
this dissertation. Also in Chapter II, a new mathematical model for the precision rolling 
element bearings has been developed in order to clarify this issue qualitatively and 
quantitatively. 

This study extends the proposed bearing formulation and demonstrates its 
superiority over the existing models in vibration transmission analyses. A schematic of 
a generic system with a flexible shaft rotating at constant speed II z , flexible casing and 

mount is shown in Figure 3.1. The shaft is supported by a rolling element bearing 
which is modeled by a stiffness matrix [K]jj m of dimension 6 as proposed in Chapter H 
The excitations at the rotating shaft are given in terms of an alternating load vector 
[f(t))sa = {FjsaW.TjsaW} 7 = { f (t) } s — [f) sm ; j=x,y,z, where Fj^t) and Tj sa (t) are the 
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x 



Figure 3.1 Schematic representation of the vibration transmission problem. Here the 
flexible shaft is subjected to alternating forces Fj sa (t) and torques TjsaW 
where j = x, y or z, is the direction and subscript a implies alternating. 
Also, 0 is the angular displacement and u is the translational displacement 
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alternating force and torque respectively, {f(t)} s is the total load vector of dimension 6, 
{f)sm represents the mean load vector, and superscript T implies the transpose. In the 
vibration analysis, {f} sm and bearing preloads are not included as they do not appear in 
the governing equations of the linear vibration model but are used for computing [K] bm . 
The effect of bearing coupling coefficients, which are off-diagonal and rotational 
diagonal terms of [K] bm as described in Chapter II, on the eigensolution, forced 
vibration, and vibration transmission through bearings is evaluated. Our theory will be 
illustrated and validated through 3 physical system example cases; experimental 
verification is also included. 

3.2 LITERATURE REVIEW 

The existing bearing models which assume either ideal boundary conditions [10- 
12] for the shaft or translational stiffness elements [13-15] have already been discussed 
in Chapter II. Various formulas for estimating translational stiffness coefficients 
commonly used by researchers have been compared with our proposed [K] bm 
formulation. These simple bearing models are widely used in vibration models of the 
rotor dynamic systems, which typically exclude casing and mount dynamics, to calculate 
critical speeds, responses due to shaft excitations such as mass unbalance and gear 
transmission error, and dynamic stability [10-15]. In most of these cases, the vibration 
transmission through bearings is never or not the primary issue, and moreover such 
vibration models are clearly inadequate in this regard. A satisfactory study of the 
vibration transmission through bearings in rotating system similar to Figure 3.1 is yet to 
be reported even though several researchers have considered it with limited success [20- 
22,24]. In 1979 White [20] evaluated the rolling element bearing vibration transfer 
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characteristics using a two degrees of freedom (DOF) vibration model of the system 
shown in Figure 3.1. His formulation is based on only the radial bearing stiffness 
coefficient . He concluded that an increase in preload increases and system 
natural frequencies. He also found that the effect of bearing nonlinearity is negligible at 
higher preloads. In 1987 Kraus et al. [24] proposed a single degree of freedom model 
for a similar physical system (with a very compliant mount) to estimate k^ from 
measured vibration transmission spectra. In both of these studies, the coupling 
coefficients of [K]b m are not included. 

In 1982 Rajab [21] philosophically proposed a bearing stiffness matrix which 
consists of kbn-, kbr© and kbee coefficients. Some of the key features of his model are 
also summarized in Reference [22]. This model is in fact a subset of our [K]b m as 
shown in Chapter II of this dissertation. Several errors in his analytical development 
have already been pointed out. In addition, he incorporated his bearing model in a 
system study using a commercial structural synthesis program [38]. However, based 
on our study we have inferred that he incorrectly synthesized the system model given the 
plate experimental modal data, shaft finite element model and analytical bearing model. 
Moreover, an error in this case was committed when he converted kb r e and kbee 
coefficients to "effective stiffness coefficients" which he claimed to couple the shaft 
bending motion to the plate out-of-plane motion. This method excludes the bearing 
rotational degree of freedom which obviously is not correct 

3.3 ASSUMPTIONS AND OBJECTIVES 

Linear discrete vibration models of the generic system shown in Figure 3.1 are 
used to incorporate [K]b m and to characterize the vibration transmission through rolling 
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element bearings. The stiffness coefficients of [K] bm are evaluated using the analytical 
expressions presented in Chapter II of this dissertation. Effect of the gyroscopic 
moment on the shaft dynamics is not included. Since the bearing system is statically 
indeterminate, the direct stiffness formulation technique is used to obtain the system 
governing equations as opposed to the flexibility formulation. The governing equations 
for the system vibration model can be given in the matrix form as 

[M]{q(t)} a +[C]{q(t)} a + [K]{q(t)} a = {f(t)} a (3.1) 

where [M], [C] and [K] are the system mass, damping and stiffness matrices 
respectively, and {q(t)} a and {f(t)} a are defined as the generalized alternating 
displacement and applied load vectors respectively. Due to the linearity of the vibrating 
system, mean shaft loads (f)bm an< * preloads do not directly affect the dynamic 
response of the rotating system and hence are excluded from equation (3.1). However, 
{ f)bm and bearing preloads are assumed to be constant to ensure a time-invariant [K] bm 
matrix which depends only on these mean loads or on the mean deflection operating 
points. Accordingly, only the alternating shaft loads {f(t)}sa in Figure 3.1 which 
represent typical machine excitation due to the kinematic errors, mass unbalances and 
torque fluctuations are included in the forced vibration problem. The energy dissipation 
associated with the rolling element bearings is assumed to be an energy equivalent 
viscous damping matrix [C] b = a [K] bm where a is the Rayleigh damping matrix 
proportionality constant. Dynamic instabilities due to the oil whirl phenomenon and 
asymmetry of rotating elements [1 1,12] are clearly beyond the scope of this study and 
hence are not considered here. 
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The specific objectives of this chapter are to: (i) incorporate the proposed bearing 
matrix [K] bm , developed in Chapter II of this dissertation, in the linear discrete vibration 
model of the rotating mechanical equipment as described by equation (3.1) using both 
the lumped parameter and dynamic finite element methods, (ii) evaluate the dynamic 
stability of the proposed bearing system model using the Liapunov's second method, 
(iii) calculate eigensolution and forced harmonic responses, and predict vibration 
transmission through rolling element bearings for three example cases, (iv) demonstrate 
the advantages of our formulation over the existing models by Kraus et al. [24] and 
White [20], and (v) validate the proposed theory by comparing analytical prediction with 
experimental data on an analogous system. 

3.4 SYSTEM GOVERNING EQUATIONS 

3.4.1 Method A: Lumped Parameter Model 

The proposed bearing matrix [K] bm can be easily implemented in equation (3.1). 

Note that, the coupling coefficients of [K] bm provide the capability to predict casing 

rigid body angular 0j ca (t), j=x,y,z, and translational uj ca (t) motions given only the 

unidirectional transverse shaft forces. Hence we can couple the shaft motions to the 

motions of a casing of a system similar to Figure 3.1 but with rigid shaft and rigid 

casing using a lumped parameter model. The bearing preloads can now be included in 

the mean shaft load vector [f} bm by a direct vector addition as the rigid shaft can be 

assumed to be a single lumped mass for this purpose. An alternating displacement 

vector { q(t) } a = {{q(t) } J a , (q(t)} ^ ) T , is defined where [q(t)} sa = {uj sa (t),0 jsa (t)} T 

T 

and { q(t) } ca = { Uj ca (t),0j ca (t) } , j=x,y,z, are the shaft and casing alternating 
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displacement vectors respectively. The governing equations of motion for this generic 
vibration model with DOF = 12 are given by equation (3.1) with 


[M] = 


[M] 

[ 0 ] 


[ 0 ] 

[M] 


c J 


[C] = O [K] 


[K] = 

bm 


k. k 
bxx 


[K] = _ 

[K] 

bm 

[K3 k 

bm 

-[K] 

[ K ] + 

bm 


{f(D> 

={ {f(,)} 4 

a l {0} J 

bxy 

k bxz 

k bx0 x 

k. 0 01 

bx0 y 

byy 

k 

byz 

k by0x 

k , * 0 

by0 y 


k 

bzz 

k bz0 x 

o 

>* 

1 



k b0 x 9x 

k b0 x 0 y ° 


bm 


symmetric 


k b0 y 0 y 


0 

0 


(3.2a, b) 


(3.2c, d) 


(3.2e) 


where the stiffness matrices [K] bm and [K] v pertain to the bearing and mount 
respectively, and the matrices [M] s and [M] c are diagonal shaft and casing mass matrices 
respectively; each matrix of dimension 6. Specific examples of this method along with 
the eigensolution and forced response studies will be presented in Sections 3.7 and 3.8. 


3.4.2 Method B: Dynamic Finite Element Formulation 

Consider the dynamic finite element method of incorporating [K] bm in equation 
(3.1) especially when shaft and casing plate are elastically deformable over the 
frequency range of interest. This method is different from the lumped parameter 
formulation of Section 3.4.1 which assumes non-compliant shaft and casing. For 
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example, if the flexible casing plate is considered to be very large compared to the 
bearing dimensions, then the bearing nodal point on the shaft can be coupled to only one 
bearing nodal point on the plate as shown in Figure 3.2a. Accordingly, the present form 
of [K] bm is implemented in the finite element model as a generalized stiffness matrix like 
the lumped parameter model. On the other hand, when the flexible casing plate 
dimensions are finite and of the order of bearing dimensions, then several bearing nodal 
points are considered as shown in Figure 3.2b. The discretization philosophy here 
assumes that a relative displacement vector, given by the difference between the 
averaged displacement vector of bearing nodal points on the plate and the displacement 
vector of a bearing nodal point on the shaft, is equivalent to the actual rigid body bearing 
motion. Accordingly, we divide the bearing stiffness coefficients equally among all the 
generalized stiffness elements connecting the bearing nodal points on the plate to a single 
bearing nodal point on the shaft. In the limit, where all the bearing nodal points on the 
plate are collapsed to a single nodal point, [K] bm is recovered as in the first method. 
Our finite element formulation uses conventional structural elements typically available 
in commercial software programs [39] - this will be illustrated in Section 3.9. Other 
features of this method are similar to those discussed earlier in Section 3.4.1. 

3.4.3 Other Methods 

Alternate methods of incorporating [K] bm in equation (3.1) such as finite 
difference which is similar to method B, flexibility, component mode synthesis and 
transfer matrix formulations are also possible. In the flexibility formulation, the bearing 
flexibility matrix can be obtained by inverting a subset of the bearing stiffness matrix 
[K] bms which excludes zeroes conresponding to 0 Z angular direction from [K] bm . In the 
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Flexible shaft 


(a) Single Bearing Nodal Point on the Plate 



Figure 3.2 Discretization method for implementation of [K] bm in finite dement model 
of a system similar to Figure 3.1. Here {q(t)} sa and {q(t)} ca are the 
alternating shaft and plate displacement vectors at bearing location 
respectively. Subscript a implies alternating component 
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transfer matrix method, the field matrix [T] b for a bearing can be easily related to 
[Klbms* 


[{q<0} 


= [ T] , 


f{q(0} 


ca 


[T]. = 


[I] 

[ 0 ] 


[K] 

[I] 


- 1 
' bms 


(3.3a,b) 


where { (q(t) (f(t) } T and { (q(t) }^, {f(t) } T are now the state vectors 

at bearing locations on the shaft and casing plate respectively. Equation (3.3) can now 
be integrated with transfer matrices of the shaft and plate which are well documented in 


References [40,41]. Direct application of these alternate methods [40-43] are beyond 
the scope of this paper and are left for further research. 


3.5 BEARING SYSTEM STABILITY 

The stability of the proposed linear, non-gyroscopic model of a bearing system 
similar to Figure 3.1, which is governed by equation (3.1) with [f(t)] a = [0], can be 
determined using several techniques such as Liapunov’s stability method, Routh- 
Hurwitz criteria, or from the direct evaluation of system eigenvalues. Here the 
Liapunov’s second method is used for its simplicity when applied to such a vibration 
model [44,45]. If the system matrices [M], [C] and [K] of equation (3.1) are always 
symmetric and positive definite, then the system is asymptotically stable per Liapunov. 
The first condition is directly satisfied since [M], [C] and [K] are symmetric. Further, 
since [M] is diagonal and consists of only positive entries, it is clearly positive definite. 
For [K], the positive definite test can be performed by evaluating its principal minor 
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determinants which is demonstrated here for the generic lumped parameter model. 
Consider the decomposition of [K] given by equation (3.2) into a product of 3 matrices 



where the square submatrix [0] of the appropriate dimension consists of only zero 
entries, [I] is an identity matrix of the same dimension, and [K] bm and [K] v have been 
defined in Section 3.4.1. The determinant of [K] is the product of the determinants of 
the three matrices on the right hand side of equation (3.4) 

|K| = |[K] V ||(K]J (3.5) 

If [ K] > 0, w j = 1,2,. ...P and P = 1,2 12, [K] is positive definite. Other 

wj 

principal minor determinant with P<12 can be obtained by excluding the stiffness 

coefficients which are not entries in the principal submatrix of equation (3.5). Since 
| [ K v ] | = k vx *k v y«k vz *k v e x *k v 0 y * k V 0 z > 0, it implies that equation (3.5) is positive 

only if [ K ] is positive. We may recall that [K] bm has zero entries in the last row 
bin 

and in the last column corresponding to torsional 0 Z angular direction which forces the 
bearing system to be semidefinite. Now define a new matrix [K] bms of dimension 5 as 
a subset of [K] bm with these zeroes excluded. If [K) bms has positive principal minor 
determinants, then this system is dynamically stable because it will consist of stable 
oscillations superimposed on the mean shaft rotational motion Q z * Q z (t). Further, it 
follows that [C] bms which is proportional to [K] bms is also positive definite if [K] bms is 
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positive definite. The resulting equation obtained from the expansion of the determinant 
of [K] bms in terms of its entries k bw j, wj=l,2,3,4,5, is given as follows in terms of the 
stability functions <l>j, j = 1,2, 3, 4. 

[ K] =2 0 4 (55,34,45,35) {k b24 0 4 (ll,23,12,13)+k bl4 0 4 (22,13, 12,23)} + 

bms 

2 0 4 (44, 35,45,34) {k b25 (D 4 (ll,23,12,13)+k bl5 0 4 (22,13,12,23)}+ 

2 0 4 (24,35,25,34) {k bl4 0 4 (25,13,23,35)+k bl 5 0 4 (12,34,24,13)} - 
k b i5 k b22 k b 3 4 0 3 (14, 35, 15,34) + k b24 k bl5 k b i 4 0 3 (23, 35,33, 25) + 

2 k bl4 0 4 (25,45,24,55) 0 4 (13, 23,12,33) + { k bl4 k b25 } 2 k b33 + 
Oi(l,2) <D 2 (45,34,35, 44, 55) +02(1,3) 0 2 (45,24,25, 44,55) + 
k b i 5 Oi(2,3) 0 3 (14, 45, 15,44) +k bll { 0 4 (24, 35, 25,34) } 2 + 


{k b i 4 } 2 { k b23 0 3 (23, 55,35, 25) - k b22 02(3,5) } + 

Oi(4,5) { k b ii Oi(2,3) + O z (23, 13,12, 33,22) } + 

2 k bl5 0 4 (24, 45,25, 44) 0 4 (23, 13,12, 33) + 

2 k b i 5 0 4 (14, 25, 15,24) 0 4 (23,34,24,33) (3.6a) 

Oi (wi,w n ) = k bwiWl k bw ^ w ^ — { k bWlWn } (3.6b) 

0 2 (wi,wn,wiii,wiv,wv) = 2k bwi k bwn k bw ^— k bw ^ { k bw ^ } ~k bw ^, { k bw ^ } (3.6c) 

0 3 (Wj.wn.wm.wjv ) = 2k bWl k bwn - k bwm k bWlv (3.6d) 

^4 (wi,w n ,w in ,w IV ) = k bwi k bwn - k bWnI k bW][V (3.6e) 


where wj, j = I,II,III,IV,V, are the dummy variables and each wj may represent either a 
single number or a set of two numbers in equation (3.6a). The principal minor 
determinants of [K] bms can also be derived from equation (3.6) by excluding the 
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appropriate stiffness coefficients which are not entries of the particular principal 
submatrix. Hence, the bearing system is stable if each principal minor determinants 
derived from equation (3.6) is positive. The stability of the proposed bearing model 
given in Chapter II of this dissertation can now be verified using these conditions. 
Inequalities associated with the stability criteria for these models are summarized in 
Tables 3.1 and 3.2 for ball and roller bearings respectively. These inequalities arise due 
to the fact that we are yet to impose any restrictions on these stiffness coefficients. In 
Chapter II of this dissertation, these stiffness coefficients are given as functions of 
bearing kinematic and design parameters, and hence any coefficient can not assume an 
arbitrary value as it is related uniquely to other coefficients through these parameters. 
Extensive numerical studies performed over a wide range of these parameters have 
indicated that the bearing models proposed in Chapter II are indeed stable provided the 
preloads are sufficiently large to avoid the clearance non-linearity. Figures 3.3-3.5 
illustrate examples of these bearing system stability studies for precision rolling element 
bearings whose design data are given in Table 3.3. In all of these figures, the stability 
functions Oj given in Tables 3.1 and 3.2 are found to lie within the stable region. 

3.6 SYSTEM RESPONSE 

The eigensolution of the linear, non-gyroscopic undamped system, formulated by 
setting { f (t) } a = {0} and [C] = [0] in equation (3.1) given by method A or B, yields real 
valued natural frequencies C0j, j=l,2,3,..., and the modal matrix [U] = {<t>} 2 > 

.... ...] for the stable system. Since the system is proportionally damped, the 

modal damping ratio is t,. = o co? /2 and the damped natural frequency is 



64 


Table 3. 1 Bearing system stability criteria for the proposed ball bearing model 
(j =xor l,y or2; p = xorl,yor2 but p^j) 


Proposed bearing model^ 
see Table 2.2 

Stability criteria^ 

kxX’kyy>k zz ,k0 x Q x ,kQyQy,k z 0p 

<J>!(3,p+3) > 0 

kxx>kyy,k zz ,k0 x 0 x ,k0y0y,k x 0y,ky0 x 

Oj(l,5) > 0 ; Oj(2,4) > 0 

kxx»kyy.k zz ,k0 x 0 x ,kQy0y 

always stable 

kxx»kyy>k zz .k0 x 0 x ,k0y0y,kp z 

<l>i(p,3) > 0 

kxx>kyy>k zz ,k0 x 0 x ,k0y0y,k X y, 

k0x0y*k Z 9x’kz0y 

0 1 (1,2)>0;0 1 (3,4)>0; 
k 33 <I> 1 (4,5)+<D 2 (4 5,3 4,3 5,4 4,5 5) > 0 

kxx»kyy,k zz ,k0 x 0 x ,k0y0 y ,k x 0 x , 

ky0y,kj z ,k z 0p 

ki ^ 22 ^ 1 ( 3 , 4)+kj3 2 kppk44+k 1 4 2 k 2 2k33 > 0 ; 

<J>!(j,3)>0; 

{k 33 0 1 (l,4)<I> 1 (2,5)0 5 (p p+3,j+3 j+3,p p,4 4,5 5)} 
{kj 3 k 3 p+3) 2 • {O s ( j j+3,p p,l 1,2 2,j+3 j+3)} > 0 

kxx»kyy>k zz ,k0 x 0 x ,k0y0y,k X y, 
kxz’ky Z k0 x 0y 

Oj(4,5)>0; 0^1,2) >0; 
k n 0 j(1,2)+0 2 (2 3,1 2,1 3,2 2,3 3) > 0 

kxx’kyy>k zz ,k0 x 0 x ,k0y0y,k x 0y, 

ky0>;>kjz,kz0p 

«> 1 (p,j+3)>0;^ 1 (j,3)>0; 
kjjOi(3,p+3)+<I > 2(3 p+3j 3j p+3,3 3,p+3 p+3) > 0 

all non- zero except 0 Z terms 

\[K] |>0 

1 bmsl 


t Here the subscript b which implies bearing has been omitted for brevity, 
tf Stability functions Oj and <J> 2 are defined by equations (3.6b) and (3.6c) and is 
given by <D 5 ( v/ h 'w n ,vr m ,v/ l v,vrv ) = {k bWl ) 2 k b w n - k bWm k 5wjv k bWy 
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Radial Deflection (mm) 


S3(ball) 
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Radial Deflection (mm) 
(b) 


Figure 3.3 Plot of stability criteria functions for ball and roller bearings subjected to 
mean bearing radial deflection 5 xm . (a) SI = Oj(2,4) and S2 = 0^1, 3). 
(b) S3 =k 511 O^) + <D 2 (35, 13,15,33,44). 
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Axial Deflection (mm) 


Figure 3.4 Plot of stability criteria functions for ball and roller bearings subjected to 
mean bearing axial deflection 8 zm . Here, SI = <I> 1 (2,4) and S4 = OjCl.S). 




Stability Function (MN /m) Stability Function (MN 2 ) 



Angular Deflection (rad) 

(b) 

Figure 3.5 Plot of stability criteria functions for ball and roller bearings subjected to 
mean angular misalignment p ym . (a) SI = 0^2,4) and S2 = 0^1,3). (b) 
S3=k bll <D 1 (3,5)+d> 2 (35,13,15,33,44) 


cT - £ 
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Table 3.3 Design parameters for typical ball and roller bearings used for system studies 


Parameters 

Ball type 

Roller type 

Load-deflection exponent n 

3/2 

10/9 

Load-deflection constant K n (N/m n ) 

1.0 E9 

1.0 E8 

Number of rolling element Z 

12 

14 

Radial clearance r L (mm) 

0.00005 

0.00175 

Pitch diameter (mm) 

40.05 

38.00 

Ao (mm) t 

0.05 

- 

Unloaded contact angle Oq 

40° 

15° 


t Unloaded distance between inner and outer raceway groove curvature centers 
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given by (0 ^ = (0 ^1 - Z,. . Free vibration response due to the initial conditions is 

not considered as only the steady-state particular solution corresponding to sinusoidal or 

periodic load vector {f(t)} a is of primary interest Define excitation by the Fourier series 

expansion as {f(t)} a = X {f) a D el(0pt where 0) D = pco 0 > co 0 is the fundamental 

P 

frequency, and {f} ap is the complex Fourier coefficient load vector. The steady-state 
particular solution {q(t)} a is given by the normal mode expansion technique [40,42,43] 
as 


{q(t) } a = [U] 


iw P t 


’ ^jP 


{<M T {f} 

lY 1 j 1 1 ap 


to?- coJ+ i2C .ci).co n 

J p J p 


; p=l,2... (3.7a, b) 


An alternate approach would be to assume the harmonic solution for the alternating 

displacement as { q(t) } a = X {q} a n el0>pt - Substituting this and {f(t)} a definition into 

P 

equation (3.1), we get 

Adjf- 0)2[M] + X[K\] 

{q,a -> = |[-^[M1^[K)]| (f ’^ ; " =1+iO “0 <3 ' 8a ' b) 

where the operator Adj refers to the adjoint of the dynamic stiffness matrix. Features of 
this method are summarized in References [40,42,43]. Since the vibration transmission 
across the bearing is the primary issue, we now define sinusoidal load transmissibility 
R(pco) terms between two arbitrary locations I and II as 
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R f f (® P ) = 
*wU’ jiU 


f wIa(°V 


f jna^ M P^ 


w,j = 1,2 6 


(3.9) 


where f wIa and fj IIa are components of the dynamic load vectors at two arbitrary 
locations I and II respectively. The accelerance A(a) p ) and mobility V(co p ) transfer 
functions with motion at location I due to an alternating force or torque fju a applied at 
location II on the shaft are 


A 

q wU 



) = 





f jfla 


V(co n ) = 
v P ICO, 


A(co p ) 


w,j = 1,2,. ..,6 


(3.10a) 

(3.10b) 


where q wIa is a component of the acceleration vector at location I. Other frequency 
response functions can also be defined in a similar manner [40,42,43]. 


3.7 EXAMPLE CASE I: RIGID SHAFT AND PLATE SYSTEM 
3.7.1 Vibration Models 

Consider the mechanical system shown in Figure 3.6a which is assumed to be 
freely suspended or softly mounted such that [K] v = [0]. A ball bearings (see Table 
3.3) with constant axial preload is supporting a short rigid shaft subjected to a mean 
torque T zsm * T zsm (t) and a sinusoidal radial force F rsa (t) = F rsal e ia) ° l applied very 
close to the bearing. A lumped parameter model with DOF = 12 is proposed in Figure 
3.6b. Conversely, the same system has also been analyzed by Kraus et al. [24] using a 
simple vibration model with DOF = 1 as shown in Figure 3.6c with only kb rr 


Rigid Casing 


(a) 



Ball Bearing 


Rigid Shaft 


Free boundary condition 


(b) 


Rigid Shaft { 
Ball Bearing { 

Rigid Casing { 


Rigid Shaft | 
(c) Ball Bearing j 

Rigid Casing { 



DOF = 6 


DOF = 6 


Figure 3.6 Example case I: freely suspended rigid shaft, ball bearing and rigid plate 
system subjected to alternating radial force F rsa (t) applied at the shaft, (a) 
Physical system, (b) Proposed multi-degree of freedom vibration model 
with DOF = 12. (c) Simple model by Kraus et al. [24] with DOF = 1. 
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coefficient The bearing stiffness matrix for an axially preloaded ball bearing in Figure 
3.6b has non-negligible stiffness coefficients k bxx , k byy , k bzz , k b e x0x , k b0y0y , k bx0y 
and k by0x which are functions of the mean axial preload as given in Chapter II of this 

dissertation. The system matrices of equation (3.2) can be modified for this case by 

suppressing other bearing stiffness coefficients and [K] v . It can be easily observed 

from equations (3.1) and (3.2) that 5 sets of uncoupled differential equations exist The 

simplest three sets are homogeneous and pertain to the rigid body torsional motions 

0 ^( 1 ) of the shaft and casing, and axial vibration u za (t) of the shaft-casing system which 

are of no interest here. The remaining two sets are almost identical and associated with 

either {u xa (t),0 ya (t)} T or {u ya (t),0 xa (t)} T degrees of freedom for rigid shaft and casing. 

If the coordinate system is chosen such that F^t) line of action coincides with the x- 

axis, then the steady-state solution to the set of differential equations in terms of 
T 

(u ya (t),0 xa (t)} is trivial. Hence, the problem reduces to a semi-definite vibration 
system with DOF = 4. Accordingly, rewrite [M] s , [M] c and [K] bm in equation (3.2) in 
terms of the displacement vector {q(t)} a = { u xsa(*)>® ysa(*)> u xcaW»®ycaft)}^ as 





f m_ 0 I 

[M] = 
s 

[o\] 

; [M] = 
c 

[o I C J ; [W bm = 


k k 

bxx bx0 y 
k bx8 y k b0 y 6 y 


(3.11) 


First two eigenvalues corresponding to the rigid body motions in x and 0 y directions are 
zero. The dimension of equation (3.11) is further reduced to DOF = 2 by defining 
relative motions 8 xa (t) = u xsa (t) - u xca (t) and (3 ya (t) =0 ysa (t) - 0 yca (t) which turn out to 
be the bearing rigid body motions. 
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7 

3 

CO 

0 

f* 

xa 

L 0 

P , 

L yaj 


7m ~ m s + m c 



(3.12a) 


(3.12b) 


It may be noted that purely translational 5 xa (t) model by Kraus et al. [24] shown in 
Figure 3.6c constitutes a subset of equation (3.12) with k|j X Q y = 0. Eigensolution of 
equation (3.12) with [C]=[0] yields the following natural frequencies C0j and modes 


B 1 ±B 2 

C ° 1 . 2 J2y m ™ s yi l s ’ Bl Ymmsk b0 y e y + Y i I s k bxx 


(3.13a,b) 


B 2 J m m s k b 0 y 0 y Y I I s k bxx) +4 7 m m s 7 I I s k 


2 

bx0. 


T 

{<D> = 


1,2 yj m s + Yj’s 8 3 


{ 1 • B 3 } 


(3.13c) 

(3.14a) 


B 3 = 


7 m m,,k LA . - y I k ±B„ 

' m s b0 y 0 y ' I s bxx 2 


2y r k 


I bx0. 


(3.14b) 


On the other hand, the eigensolution of the single degree of freedom system is given by 

r . T 


A 

CO = 
2 


bxx 
Ym m s 


and {0} = 


l { VYm m s } 


where subscript 2 and superscript A are 


chosen to indicate that this solution essentially estimates C0 2 aud W 2 * n equations 
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(3.13) and (3.14). Since a >2 does not include k bx Q y and k b Q y 0 y , ©2 < w 2 as evident 
from Table 3.4 for 3 different axial preloads. This natural mode is dominated by MO 
as indicated by equation (3.14). And, the first mode {<(>} 1 which is predominantly 
P ya (t) is also affected by the axial preload. 


Table 3.4 Bearing stiffness coefficients and undamped natural frequencies (Hz) of 
example case 1 ^ 


Axial 

preload 

Fzbm(N) 

Bearing stiffness coefficients 

Proposed model 
(DOF = 2) 

Simple 

model 

(DOF=l) 

A 

C0 2 

k b xx (N/m) 

k b x0 y (N) 

k b e y e y (Nm) 

tOj 

to 2 

115 

1.84 E7 

-3.05 E 6 

1.36 E4 

156 

372 

341 


2.13 E7 

-3.12 E 6 

1.70 E5 

191 

395 

367 

285 

2.43 E7 

-3.09 E 6 

2.02 E5 

221 

416 

392 


t Other system parameters are: m s = 10.0kg, Is=0.025kgm 2 , Y m =0.4,y s =0.3,a=lE-6s. 


3.7.2 Bearing Transmissibility 

The forced harmonic response of equation (3.12) can be obtained using the 
dynamic stiffness approach given by equation (3.8) 
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5 (t) = 

xa v ' 


( k b0 y 0 y ^ “oY, Is ) Ym F xsal e 


i© 0 t 


( k b9 y e y ^ 0) o'Y I I s ^( k bxx^ ® oYm m s) ( k bx0 y k) 


(3.15a) 


P (t) = 

r ya v ’ 


- k X 7 F e it0ot 
bx0 y *m r xsai e 


( k b0 y 0 y ^ “oVs^bxx^ ®o Ymin s) ( k bx0 y ^) 


(3.15b) 


It can be seen from equation (3.15) that F xsa (t) not only excites 8^(0 but P ya (t) as well 
which is obviously not predicted by Kraus et al. [24]. The steady-state solution for this 
simple model is given by the following; compare it with equation (3.15a). 



(3.16) 


Both models are used to determine the load transmissibility magnitude terms R(co 0 ) 
which are computed using (S^O^py^t)} 7 , [K] and [C]. Dynamic bearing force F xba (t) 
and moment M yba (t) magnitudes excited by the shaft force F xsa (t) are given by force 
transmissibility ^F xba ,F xsa (° ) o) an d moment transmissibility ^M yba ,F xsa ( <B o) 

respectively. 


R 


F F o) 

xba’ xsa 


.J* 


Y mA /( 1 + a2t0 o) 2B 4 + ( 1+<y2c0 o) B 5 


V (1+a2(0 o 


) 2 B 4 + (l + a 2 coJ)(B 5 + B 6 ) + B 7 


(3.17a) 
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R M ,F 

yba xsa 


YmVs^o J 1 


(l + oV 0 )k 2 bxe ^ 


^/(1 + c 2 co 2 ) 2 B 4 +(1 + ct 2 coJ)(B 5 +B 6 ) + B 7 


B = [ k 


k - k 
b0 y 9 y bxx bx0 y 


/• 

)■ 


- 2k k 
, y bxx be y e y 


) 


(3.17b) 

(3.17c) 

(3.17d) 


B 5 = ® oYi I s k bxx( k bxx Y I I s a) o + 2k bx 0 
B 6 = CO o'Ym m s k be y e y ( k bG y 0 y Yn > m s co o + 2k b x0 y " 2k bxx k b0 y 0 y j ( 3 - 17e ) 
B 7 = 2co oYm m { 2k bxx k b0 0 ' C °o(Yi I s k bxx +Y m m ^bGyBy)} + 


03 o Ym m ^ s Ym m s 7 j ^ s® o 2k bx0 y ^ 


(3.17f) 


Only the force transmissibility, as given below, is predicted by the simple model [24]; 
compare it with equation (3.17a). 


R f ,f (®o) 

xba xsa 


k b„^ V( i + o2<a o) 


V( k bxx - ®0Tn,™ S ) +( k bxX C “) 


(3.18) 


Figure 3.7 compares equations (3.17) and (3.18). Our model predicts higher 
Rp^b a> F xsa (®°) and o >2 than the simple model due to the additional constraints imposed 

by [K] bm . Also, it is clear that the simple model can not predict dynamic moment 
transfer through the bearing. The bearing transmissibility functions R(co 0 ) predicted by 
our model for 3 different axial preloads are shown in Figure 3.8. Note that the resonant 
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Figure 3.7 Bearing transmissibility spectra Rf wba F xsa (0) o ) for example case I. Here, 
Rl: force transmissibility with f w ^ = F x ^ and R2: moment transmissibility 
with f w b a = Myb a , as predicted by our model with DOF = 2 and the simple 
model by Kraus et al. [24] with DOF = 1. 



Moment Transmissibility (m) Force Transmissibility 
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(a) 



Figure 3.8 Effect of mean axial preload on the bearing transmissibility spectra 
Rf . p (©o) defined by equation (3.17) for example case I. (a) Force 

wD«i* xsa 

transmissibility with f w b a = F x j, a . (b) Moment transmissibility with f w b a = 

Myb a . 
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amplitudes vary depending on the amount of axial preload and the resonant frequencies 
increase with increasing preloads as expected. 

3.8 EXAMPLE CASE II: RIGID SHAFT AND PLATE SUPPORTED ON 
FLEXIBLE MOUNTS 

3.8.1 Vibration Models 

The physical system of example case I is modified to include flexible mounts, [K] v 
* [0], and mean radial shaft force F rsm * F rsm(0 and F rsa (t) = F^e^ot as shown in 
Figure 3.9a. The ball bearing (see Table 3.3) is also preloaded in the axial direction. 
This is modeled using lumped parameter theory with finite mount stiffness coefficients 
k v j, j = x, y, z, 0 X , 0 y or 0 Z as illustrated in Figure 3.9b. Since the rigid shaft 
assumption still holds, F^ = F rbm is applied directly on the bearing in a manner similar 
to the axial preload. The bearing matrix [K] bm in Figure 3.9b has non-negligible 
stiffness coefficients k bxx , k b yy, k bzz , k b 0 x e x » k b 0 y 0y, k bx Q y , k b yQ x k bxz and k bz Q x 

which are functions of the mean bearing load vector {f} bm = {f} sm or the mean bearing 
displacement vector {q} bm = {q) sm - {q} C m as given by Chapter II of this dissertation. 
Conversely, White [20] has investigated this problem using a simple model with DOF = 
2 as shown in Figure 3.9c. It may be noted that his bearing system model did not 
include the effect of axial preload. 

Like example case I, the governing equations (3.1) and (3.2) can be modified and 
reduced to 4 uncoupled sets of differential equations. The first two sets associated with 
0 ysa (t) and 0 yca (t) are homogeneous. The third set is similar to equation (3. 1 1) but with 
x and y subscripts interchanged, and two mount stiffness coefficients k vy and k v g, 

included. However, it is still homogeneous and therefore has only trivial steady-state 
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Figure 3.9 Example case II: rigid shaft, ball bearing and rigid plate system supported 
by flexible mounts and subjected to alternating radial force F xsa (t) applied at 
the shaft, (a) Physical system, (b) Proposed multi-degree of freedom 
vibration model with DOF = 12. (c) White's vibration model [20] with 
DOF = 2. 
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solution. The final set corresponds to {q] a (t) = (u xsa (t), u zsa (t), 0 ysa (t), u xca (t), u zca (t), 
e yca (t)} T with [M] s , [M] c , [K]b m and [K] v of equation (3.2) reduced to 




f m s 

0 0 1 




'm 

c 0 

0 ‘ 




0 m s 0 


; [Ml, 

* 

0 

m c 

0 




Lo 

0 I.J 




0 

0 

h. 



' k bxx 

k bxz 

^bX0y 


'k 

vx 

0 

0 

-] 

[K] = 

k bxx 

k 

bzz 

k bz0 v 


; 1KI V = 


0 

k vz 

0 


bm 



y 




0 

0 

k.. Q 


k 

k 

1c 






VO 

y J 


bx9 y 

bz0 

y J 








(3.19a,b) 


(3.19c, d) 


The vibration model by White [20] may be formulated by retaining only two equations 
corresponding to u xsa (t) and u xca (t) and excluding all bearing and mount stiffness 
coefficients except for k^ xx and k vx . 

Analytical eigensolution of the undamped system is not possible since it requires 

solving for the zeroes of a 6-th order polynomial in co 2 . Therefore, this problem is 

solved numerically using an eigenvalue routine [42]. Using the same system parameters 

as in example case I with 3 different F xsm , natural frequencies and modes are found as 

given in Tables 3.5b and 3.5c for both our and White's models. Corresponding bearing 

T 

mean loads {F x bm’Fzbm>Mybm} and relevant bearing stiffness coefficients computed 
using the method derived in Chapter II are listed in Table 3.5a. Tables 3.5b and 3.5c 
indicate that only the first and fifth modes of our model are predicted by White's model; 
here superscript A is again used to denote estimation based on the simple model. The 
first natural frequency C0j predictions by both models are very similar. But 0)5 
prediction, whose mode is similar to the second mode of example case I, indicates a few 
discrepancies. White's model also underestimates this natural frequency due to the 
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Table 3.5 Results of example case II t 


(a) Computed bearing mean loads and stiffness coefficients 


Shaft 
mean load 

F 

r xsm 

(N) 

Bearing loads ^ 


Bearing stiffness coefficients 


Fzbm 

(N) 

^4ybm 

(Nm) 

k bxx 

(N/m) 

k bxz k bx8 y 

(N/m) (N) 

kbzz k bz9y 

(N/m) (N) 

kbeyOy 

(Nm) 

44 

198 

0.62 

2.26E7 

5.84E7 -2.90E5 

8.50E7 -3.10E3 

1.71E4 

94 

216 

1.16 

2.64E7 

1.14E7 -2.33E5 

8.55E7 -1.36E4 

1.72E4 

122 

229 

1.38 

2.93E7 

1.40E7 -1.98E5 

8.60E7 -2.19E4 

1.74E4 


(b) Undamped natural frequencies (Hz) 


Mean 
shaft load 

Fxbm (N) 


Proposed model (DOF=6) 


Simple model (DOF=2) 

COi 

co 2 

co 3 

cd 4 

© 5 . 

©6 

A 

©1 

A 

“5 

44 

93 

100 

111 

287 

350 

607 

97 

321 

94 

96 

100 

115 

300 

351 

614 

97 

345 

122 

97 

100 

117 

308 

364 

620 

98 

362 


(b) Modes of vibration 


Proposed model (DOF=6) 


Simple model (DOF=2) 


(<f>h 

{<t>h 


{*}4 

Ws 

(0} 6 

($h 

{$>5 

0.217 

0.031 

-0.030 

-0.130 

0.180 

0.044 

0.218 

0.229 

-0.035 

0.203 

0.003 

0.022 

-0.036 

0.236 



1.648 

0.084 

5.966 

-0.354 

-1.245 

-0.068 



0.171 

0.084 

-0.071 

0.107 

-0.139 

-0.031 

0.187 

-0.178 

-0.027 

0.194 

0.008 

-0.018 

0.028 

-0.165 



0.167 

0.009 

0.917 

4.648 

3.295 

0.073 




t Other system parameters are: m s =10.0 kg, I s =0.025 kgm 2 , m c =15.0 kg, 
I c =0.03 kgm 2 , a=lE-6 s, k vx =lE7 N/m, k vz =lE7 N/m, k x e y =l E5 Nm. 
^xbm = Pxsm 

ttt These are for mean shaft load Fxsm ~ Fxbm - 94 N. 
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incomplete bearing stiffness model employed. Here again our predicted modes include 
5 xa , 8 za and P ya displacements of shaft and casing components which are not 
considered by White's model. 

3.8.2 Frequency Response 

The forced harmonic response solution is also determined numerically using the 

dynamic stiffness approach outlined in Section 3.6. Driving point (with j = s) and cross 
point (with j = c) accelerance spectra A q ^ a F xsa (co 0 ) are given in Figure 3.10 with q WJ - a = 

u xja> u zja or Qyja f° r both models. Here, we observe that White's [20] model 
overestimates the magnitudes of the accelerance and can not predict, unlike our model, 
axial u za (t) and 0 ya (t) angular motions on the shaft and casing. Figure 3.11 shows the 
bearing transmissibility spectra Rf wba ,F xsa (®o) for f w b a = F x j> a , F zxa or M y b a which 

indicate that transmissibilities corresponding to ^zxa and M yba are not predicted using 
White's model. Such loads also serve as mechanisms for vibration transmission 
through the bearing to the casing, in addition to F x t, a . The mount transmissibility 
spectra Rf wva ,F xsa ( C0 o) as shown in Figure 3.11 indicate that F zva and M yva are also 

transmitted to the mounts in addition to F xva due to casing motions in x, z and 0 y 
directions. The effect of mean radial bearing force F x b m on the load transmissibilities 
R(co 0 ) through the bearing is shown in Figure 3.12 for 3 different mean loads. We 
again observe that the resonant amplitudes and frequencies are mean load dependent 
through [K]b m . Similar trends show that an increase in mean load raises the resonant 
frequencies although the effects are not as pronounced as those found in example case I. 


Magnitude (1/kg) Magnitude (1/kg) 
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Figure 3.10 Accelerance spectra A q wja ,F xsa ( w o) for example case II as predicted by our 

formulation and White's model [20]. (a) Driving point accelerance with 
q w ja = u xsa- (*>) Cross point accelerance with q w j a = u xca . (c) Accelerance 
with q wja = u zsa for shaft and q^ a = u zca for casing, (d) Accelerance with 
Qwja = Qysa f° r shaft and q w j a = 6 yca for casing. 






Figure 3.11 Bearing Rf F (co 0 ) and mount Rf p (co 0 ) transmissibility spectra 

wca* xsa wva* xsa 

for example case II as predicted by our formulation and White's model 
[20]. (a) Bearing force transmissibility with f wba = F xba . (b) Mount force 
transmissibility with f wva = F xva . (c) Force transmissibility with f wba = 
F zba f° r bearing and f wva = F zva for mount, (d) Moment transmissibility 
with f wba = M yba for bearing and f wva = M yva for mount 






Moment Transmissibility (m) Force Transmissibility Force Transmissibility 
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Figure 3.12 Effect of the mean radial bearing force F xbm on the bearing transmissibility 
spectra F f wba ,F xsa ( w o) f° r example case n. (a) Force transmissibility with 

fwba = F xba- ( b ) Force transmissibility with f wba = F zba . (c) Moment 
transmissibility with f wba = M yba . 
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3.9 EXAMPLE CASE III: EXPERIMENTAL STUDY 

3.9.1 Physical setup 

The final example case examines the experimental setup of Lin [17] as shown in 
Figure 3.13a. This system is similar to Figure 3.1 and consists of a 159mm long x 
25mm diameter non-rotating shaft supported by two rolling element bearings of 25mm 
bore x 51mm outer diameter. One is supported on a rectangular plate of approximate 
dimensions 762mm x 457mm x 9mm and the second is rigidly connected to the base. 
The plate is also bolted to a massive base structure. Excitation force F ys which consists 
of a mean F ysm = 445N via a preloaded spring and an alternating F ysa (t) component 
applied transversely at the free end of the shaft using a vibrating shaker. Driving and 
cross point accelerance spectra are measured at the shaft and on the plate respectively. 
Further details of this experiment are summarized in Reference [17]. 

3.9.2 Bearing Analysis 

Initially, only the static analysis is performed to obtain [K] bm for this experimental 

system using the method proposed in Chapter II of this dissertation. The static analysis 

neglects plate flexibility; this assumption is valid since the bearing mean loads are 

sufficiently low and do not deflect the plate. The shaft-bearing system is statically 

indeterminate as shown in Figure 3.13b. The mean force F ysm on the shaft produces 

T 

mean bearing load vector [f] bm = {0,F y bm>Fzbm> M xbm>0} wk i ck depends on the 

T 

mean bearing displacement vector {q)b m = {0,8 ym ,8 zm ,P xm ,0} . The proposed 
bearing matrix [K] bm includes stiffness coefficients k bxx , k byy , k bzz , k b0x0x , k b0y0y , 

k bx0y , k by0 x > k byz and k bz0 x which are direct functions of {q} bm . In contrast, the 





^ybjm fyb 2 m 

Applied Force 



Measured Response 


(C) 


Figure 3.13 Example case HI. (a) Schematic of the experimental setup consisting of an 
overhung shaft, 2 ball bearings and a rectangular plate [17]. (b) Static 
model of shaft and bearings used for computing [K] bm . (c) Finite element 
model of the experimental system shown in (a). The generalized stiffness 
matrix elements for the bearing model are shown by dashed lines, 
connecting 1 node on the shaft to 4 nodes on the plate. 
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conventional bearing models include only k bxx , k byy and k bzz coefficients. From Figure 
3. 1 3b, the force and moment equilibrium equations for this system are 

Fybim + Fyb 2 m — ^ysm = 0 (3.20a) 

^ybim + ^yb2m — ^yb2m^l + ^ysm(^l + ^ 2 ) “ 0 . (3.20b) 

^zbim + E zb2 m = ® (3.20c) 

Since the shaft-bearing system is statically indeterminate, bending theory for the shaft 
and rigid body motion constraint in the z direction are used to estimate stiffness 
coefficients 


El Sy 2 m El 8yj m + Fy b j m (/^/6) + My bim (/^/2) - 0 

(3.21a) 

El Px 2 m — El Pxjtn + Ey bim (^lV^ + Mybim^l = 0 

(3.21b) 

^zjin ~ 8 Z2 m = 0 

(3.21c) 


Additionally, 6 nonlinear algebraic equations defined by the mean bearing load- 
displacement relations as given in Chapter II are required. These nonlinear algebraic 
equations are solved using Newton-Raphson method [36,37]. Since mean loads on 
each bearing are sufficiently large, the bearing stiffness coefficients for both bearings are 
almost identical, as listed in Table 3.6 along with other system parameters. 
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Table 3.6 Design and estimated parameters for two identical rolling element bearings 
used in example case HI 


Load-deflection exponent n=3/2 

Ao^=0.05 mm 

k bzz=1.72E8 N/m 

Load-deflection constant K n =6.92E9 N/m n 

/l=41 mm 

k bx0y =-2.56E5 N 

Number of rolling element Z=10 

72=84 mm 

k by e x =3.52E5 N 

Radial clearance r L =5.0E-5 mm 

kbxx =1.44E8 N/m 

k bze x =4.02E5 N 

Pitch diameter=38.1 mm 

kjjyy =3.69E8 N/m 

k b0x6x =4.19E4 Nm 

Unloaded contact angle ao=0° 

k byz =2.04E8 N/m 

k b0y0y =l.O2E4 Nm 


t Unloaded distance between inner and outer raceway groove curvature centers 


3.9.3 System Study 

We incorporate the proposed rolling element bearing stiffness matrix [K] bm in a 
finite element model which includes shaft and plate dynamics, using the formulation 
given in Section 3.4.2. The finite element model shown in Figure 3.13c is implemented 
with a commercial software [39]. The shaft component is modeled using 2 noded 
Timoshenko beam elements with axial degrees of freedom in addition to the bending 
motion. The plate model is constructed using 4 noded quadrilateral plate elements with 
shear deformation and rotary inertia effects. Each node has 3 translational and 3 
rotational degrees of freedom. Four generalized stiffness matrices corresponding to the 
first bearing, each matrix being equivalent to 1/4 [K] bm , are used to couple the single 
shaft node to 4 plate nodes. The second bearing connects one end of the shaft to a 
grounded node. The boundary conditions for the plate along the perimeter are chosen to 
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be a combination of ideal clamps, u xa(t) — u ya(t) — u zaW — ®xa(0 — ®ya(0 - ®zaW ~ 
and simple- supports u za (t) = 0 as shown in Figure 3.13c in order to represent the 
physical model as much as possible. Here, the energy dissipation is assumed to be 
given by the modal damping ratio £ = 0.03. A sinusoidal force F ysa(0 = F ysal el(Dot is 
applied at one end of the shaft to simulate the experiment 

Over the frequency range of 400Hz to 2000Hz, Figure 3.14 compares the driving 
point accelerance spectra A Uysa> F ysa (G) 0 ) = ii ysa /F ysa . The simple theory shown here 

represents the conventional way of modeling bearings while other the features are 
exactly the same as in proposed model. Our predictions match measured spectra very 
well. Conversely, the simple model predicts slightly higher accelerance amplitude and 

lower resonant frequency in the vicinity of 800Hz due to the incomplete bearing model 
used. Cross point accelerance spectra A U2cai F yga (C0 o ) = Uzca/F ysa are shown in Figures 

3. 15a and 3.15b where u is measured for 2 different locations on the plate as shown 
in Figure 3.13c, and excitation F ysa (t) is once again applied transversely at the shaft. 
Here, each predicted accelerance spectrum has been averaged over 4 points in the 
immediate vicinity of the measured location. Reasonable comparisons between the 
proposed model and experiment are seen. Here, the discrepancies are primarily due to 
physical setup complexities and the limitations associated with the finite element model 
in describing some of these. In Figure 3.15, the simple model is not included because it 
predicts exactly zero out-of-plane or flexural motion of the plate. Next, the cross point 
mobility level Ly is defined by averaging mean square mobility spatially over the entire 
plate and over a frequency bandwidth Aco. This level is directly related to the structure- 
borne noise or vibratory energy transmitted through the bearing. 
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Frequency (kHz) 


Figure 3.14 Driving point accelerance spectra A u (co Q ) yielded by the proposed 
model, simple model and experiment lay Lm [17] for example case HI. 




Magnitude (1/kg) Magnitude (1/kg) 




./ \! > 


v\ k 


Proposed Model 
Experiment [Lin] 


Frequency (kHz) 



Proposed Model 
Experiment [Lin] 


Frequency (kHz) 


Figure 3.15 Cross point accelerance spectra A Uzca p sa (co 0 ) yielded by the proposed 
model and experiment by Lin [17] for example case HI. In this case, the 

simple model predicts zero response u zca at the plate, (a) Point Cl (see 
Figure 3. 13c). (b) Point C2 (see Figure 3.13c). 
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(3.22) 


where ( ) implies the complex conjugate, S c is the plate surface area and Re{ } implies 
the real part of the complex number. Table 3.7 compares Ly predictions by proposed 
and simple models with experimental data. It can be seen from this table that the 
proposed model predicts the experiment quite well and the simple model fails to predict 
any plate vibration. We can therefore conclude that our model is indeed valid for 
vibration transmission analyses. 


Table 3.7 Predicted and measured cross point mobility level Ly as defined by equation 
(3.22) 


1/3 Octave band 
center frequency (Hz) 

Experiment (dB) 
[Lin] 

Proposed 
model (dB) 

Simple 
model (dB) 

400 

-102 

-105 

-oo 

500 

-92 

-96 

-oo 

630 

-95 

-94 

-oo 

800 

-88 

-97 

-oo 

1000 

-87 

-95 

-oo 

1250 

-97 

-108 

-oo 

1600 

-108 

-115 

-oo 

2000 

-106 

-107 

-oo 
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3.10 CONCLUDING REMARKS 

A new mathematical model for precision rolling element bearing has been 
developed and incorporated in linear system dynamic models using lumped parameter 
and finite element modeling techniques for the vibration transmission studies of a 
generic single shaft-bearing-plate system. Stability studies indicate that the bearing 
system is dynamically stable for most of the practical designs. Through 3 example cases 
including one experimental study, we have shown that our proposed vibration model is 
clearly superior to the models currently available in the literature. The current models 
tend to underestimate the resonant frequencies and force/moment transmissibilities, and 
overestimate the accelerance amplitudes as compared to our models. The proposed 
model also predicts how the vibratory bending motion on the shaft is transmitted to the 
casing, illustrated through coupling coefficients of the proposed bearing stiffness matrix 
[K] bm . Finally, the forced response trends indicated that increase in the mean bearing 
loads increases system resonant frequencies. We are extending this model to predict 
vibration transmission in rotating equipment with multiple shafts, bearings and gears. 
Other applications are evident as our theory is general in nature. However, it is 
restricted to linear systems. Bearing non-linearities are being examined in a parallel 
study [46]. 


CHAPTER IV 


GEARED ROTOR SYSTEM STUDIES 
4.1 INTRODUCTION 

The focus of this chapter is on the dynamic analysis of geared rotor system shown 
in Figure 4. 1 which includes a spur gear pair, shafts, rolling element bearings, motor, 
load, casing and flexible or rigid mounts. For this purpose, discrete vibration models 
are developed and used to predict vibration transmission through bearings. Also, the 
effects of bearing, casing and mount dynamics on the internal rotating system dynamic 
characteristics are investigated. Traditional analysis approaches [15,25-29] in the gear 
dynamics area concentrate on the internal rotating system and exclude dynamic effects of 
casing and flexible mounts. Moreover, simple bearing models are typically used which 
assume either ideal boundary conditions on the shaft or translational spring elements. 
Some of the limitations associated with current bearing models have been discussed 
thoroughly in Chapters II and III of this dissertation. For instance, simple bearing 
formulations can not explain how the vibratory motion may be transmitted from the 
rotating shafts to casing and support structures, and moreover can not predict the effects 
of bearing, casing and mount dynamics on the internal rotating system adequately. 
Chapter II of this dissertation also presents a new mathematical model for the precision 
rolling element bearings to clarify this issue qualitatively and quantitatively, and Chapter 
III utilizes the proposed bearing stiffness matrix [K]t, m formulation to analyze the 
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Figure 4. 1 Schematic of a geared rotor system with flexible shafts and rigid spur gear 
and pinion, and supported by 4 rolling element bearings on a flexibly 
mounted casing. The shaft at the motor end is driven at mean speed Cl 
which in turns drives the load at mean speed This physical system is 

symmetric about a plane intersecting both the driving and driven shafts. 
The y axis is assumed parallel to the gear mesh line of action in the pressure 
angle direction at mesh point. This schematic is used for example case II 
and III. 
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vibration transmission problem in a generic shaft-bearing-casing-mount system. Our 
bearing formulation has already been validated through several example cases. 

This study extends the proposed bearing formulation of Chapters II and III to 
examine vibration transmissibility in rotating mechanical equipment through two generic 
systems which are a spur gear pair drive and a single-stage rotor system with rolling 
element bearings as shown in Figures 4. 1 and 4.2 respectively. It should be noted that 

Figure 4.2 can be treated as a special case of Figure 4.1- this will be illustrated later. 
The physical systems are assumed to be driven at mean speed O 2 = in Figure 

4. 1 and O z in Figure 4.2; subscript M refers to the motor end. Each shaft is supported 
by two identical rolling element bearings which are modeled as stiffness matrices [K] bm 
of dimension 6 as proposed in Chapter II of this dissertation. These generic system may 
be excited by motor and/or load torque fluctuations, rotating mass unbalances or gear 
kinematic transmission error e(t) [29,47]; here e(t) is defined as the deviation of the 
relative gear-pinion angular position from its relative ideal location. Typical excitation at 
the shaft is defined by an alternating load vector { f(t ) } 53 = l Fj^ft) ,Tj t) } T = {f(t)} s - 
{f} S m» j=x,y,z, where Fj^t) and Tj^t) are the alternating force and torque respectively, 
{f(t)} s is the total load vector of dimension 6 , {f} sm represents the mean load vector, 
and superscript T implies the transpose. In the case of the geared drive, {fft)}^ consists 
of e(t) at the mesh point. Additionally, in the case of the single-stage rotor system, 
transverse forces F m (t) and F ysa (t) due to mass unbalances are of interest. Note that in 
the dynamic analysis {f} sm and bearing preloads are not included as they do not appear 
in the governing equations of the linear vibration model but are used for computing 
[K] bm . Other effects such as bearing coupling coefficients, motor and load inertia, and 
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x 



Figure 4.2 Schematic of the single stage rotor system with flexible shaft rotating at 
mean speed Q z and supported by 2 identical rolling element bearings on a 
flexibly mounted rigid casing. This physical system is symmetric about a 
plane intersecting the shaft axial axis and parallel to the x-z plane. This 
schematic is used for example case I. 
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casing and mount dynamics are considered as the emphasis is on the overall system 
behavior. 

4.2 LITERATURE REVIEW 

The available literature on bearing models commonly used for internal rotor and 
gear dynamic analyses has already been discussed in Chapters II and III of this 
dissertation. Although there are numerous publications [10-13,15,25-29] on the modal 
analyses of geared drives and single or multiple-stage rotor systems, very little has been 
reported directly on the force transmissibility through bearings, and the dynamic effects 
of bearing, casing and mounts on the internal rotating system. 

4.2.1 Casing and Mount Dynamics 

A comprehensive review of the available vibration studies of casing and mounts, 
without the internal rotating system, has been given by Lim and Singh in 1987 [6]. 
Other related publication not included in Reference [6] are by Gaul and Mahrenholtz [48] 
in 1984, and Smith [49] in 1988 who have developed lumped parameter models of 
flexible machine foundations, excluding the internal rotating system dynamics. Their 
studies report rigid body translational and rotational vibration responses of the casing 
when excited by forces and/or moments. However, these studies are not specifically 
directed towards geared rotor systems, and dynamic interactions between internal 
rotating system and casing-mount system are not incorporated. In addition, Lim et al. 
[50] in 1989 have also performed a dynamic finite element analysis of a helicopter 
gearbox excluding the internal rotating system. Results again show that casing rigid 
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body vibration modes are dominant over the lower frequencies for a flexibly mounted 
casing, but numerous casing plate elastic modes are observed at higher frequencies. 

4.2.2 Gear Dynamics 

Current gear dynamic models include only the internal rotating system and simple 
bearing models, and typically exclude casing and mount flexibilities [15,25-29]. In few 
instances casing and mounts have also been included , but restricted to cases where only 
purely radial and/or axial force on the bearing, and unidirectional vertical/horizontal rigid 
body motion on the casing are modeled, as summarized in Reference [6]. Such models 
are still inadequate in explaining how the vibration is transmitted from the gear mesh to 
the casing and into the machine foundation as witnessed in previous experiments 
[9,16,17]. 

% 

4.2.3 Rotor Dynamics 

Similar to the geared rotor models, most of the existing rotor dynamic models 
concentrate on the internal rotating system and address issues related to stability, critical 
speeds and rotating mass unbalance response [10-13]. But a few investigators have 
included support flexibility models [51-53]. For instance Lund and Wang [51] in 1986 
proposed an impedance matching approach to reduce the large degrees of freedom 
(DOF) required in such models. Using one example case, they reported that support 
foundation has little influence on the internal rotor resonances, but then cautioned that 
this might not be true for other systems. Vance et al. [52] in 1987 incorporated 
measured support foundation parameters in a transfer matrix model of the rotor dynamic 
system, and concluded that the omission of support flexibility may miss some of the 
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rotor critical speeds of interest. Earles et al. [53] in 1988 reported, based on a finite 
element model of the support foundation, that the rotor critical speeds shift by including 
support flexibility. In all of the above mentioned studies, the overall dynamic behavior 
of the rotor system is only investigated partially. Moreover, only a subset of journal 
bearing [K] bm coefficients is included which can not predict completely the vibration 
transmission through bearings. 

4.3 ASSUMPTIONS AND OBJECTIVES 

Linear discrete vibration models of the geared and single-stage rotor systems 
shown in Figures 4.1 and 4.2 are used to incorporate [K] bm similar to Chapter III. The 
stiffness coefficients of [K] bm are evaluated using the analytical expressions presented 
in Chapter n. Each rotating shaft is modeled as an Euler beam in the lumped parameter 
model and as a Timoshenko beam in the dynamic finite element model but the effect of 
the gyroscopic moment on the shaft dynamics is not included. The casing is assumed to 
be rigid in the lumped parameter model and flexible in the dynamic finite element model. 
The rigid casing assumption is valid for many practical designs as it may be massive and 
the rotating speeds may be sufficiently low to avoid significant elastic deformation. The 
governing equations for both discrete vibration models can be given in the general 
matrix form as 

[M] (q(t) } a + [C] {q(t) } a + [K] (q(t) } a = (f(t) } a (4.1) 

% 

where [M], [C] and [K] are the system mass, damping and stiffness matrices 
respectively, and {q(t)} a and {f(t)} a are defined as the generalized alternating 
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displacement and applied load vectors respectively. Like Chapter HI of this dissertation, 
{fJbm and bearing preloads are also assumed constant to ensure a time-invariant [K]b m 
matrix. Energy equivalent viscous damping matrix [C]*, = a [K]b m is assumed for the 
energy dissipation mechanism in the bearing where a is the Rayleigh damping 
proportionality constant. Other features of the proposed theory given in Chapter HI of 
this dissertation are retained. 

The specific objectives of this study are to: (i) incorporate the proposed bearing 
matrix [K]j, m , developed in Chapter II, in discrete vibration models of generic geared 
system of Figure 4.1 as described by equation (4.1) using lumped parameter and/or 
dynamic finite element method, (ii) conduct overall system studies by calculating 
eigensolutions and forced harmonic responses, (iii) evaluate the effects of casing and 
mount dynamics on the internal rotating system, and predict vibration transmission 
through bearings, (iv) examine the following 3 example cases: I. single-stage rotor 
system with flexible shaft supported by two identical rolling element bearings on rigid 
casing and flexible mounts as shown in Figure 4.2, II. spur gear pair system with motor 
and load inertias attached to two flexible shafts supported by four rolling element 
bearings on flexibly and rigidly mounted rigid casing as shown in Figure 4.1, and III. 
case II with flexible casing and rigid mounts, and (v) demonstrate the advantages of our 
formulation over the existing vibration models. 
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4.4 THEORY 

4.4.1 Method A: Lumped Parameter Model 

4.4.1. 1 Equations of Morion 

Both physical systems shown in Figures 4.1 and 4.2 are discretized using lumped 
parameter technique to yield equation (4.1) through the Lagrange's equations of motion 
[45]. 


_d 

dt 


( 3E, 


9E 


u 


3E 


9E 


3E 


V 3 <iw 3 < i w 


Z- J -+-T- iL +^ Le -=F w ; w = 1,2,3,... (4.2) 

3q w Bq w 3q w w 


where E T and Ey is the kinetic and potential energies respectively, 
E c = y{q} [C]{q} is the Rayleigh’s dissipation function and F w is the 
generalized force. The total system potential energy Ey and kinetic energy E T are 
obtained by adding the energy of each system component which are derived in 
subsequent sections. Equation (4.2) when applied to Figure 4.1 will lead to equation 
(4.1). It is then rearranged in terms of the following {q} to obtain the partitioned 
system matrices 


{ q (0} a = "^{ u (0}^{ u (0} 
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[M] 

= DIAG ■ 

[[“'I 



[Kill 

[Kl,2 

[K] 1 
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IDr f»c 




[C] = o [K] 


(4.3a) 

(4.3b) 

(4.3c) 
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where {0(t)} Ra (subscript R indicates rotor) consists of the alternating angular 

displacement of all rotors including gears, motor and load; (m) and {1} are the 

lumped mass and inertia row vectors respectively; and the operator DIAG { ) transform 

the row vector into a diagonal matrix with components of the vector corresponding to 

T 

the diagonal elements. The stiffness sub-matrices [K] w ^ = [K]j w , w,j = 1, 2, 3, 

consist of the appropriate terms corresponding to the partitioned {q(t)} a . For example, 

T T T 

[K] n couples { (u(t) ) (u(t) } } degrees of freedom to itself, and [K] 12 provides 

^ * sa ca 

a coupling between { (u(t) (u(t) } T and { (QCt)}^ (0(t) } T degrees of 

freedom. 

Now assume that the rotary inertia of each shaft lumped mass is negligible. The 
system matrices given in equation (4.3) are therefore simplified by eliminating {0(0}^ 
degrees of freedom from {q(t)} a in equation (4.3a). The resulting mass and stiffness 
matrices are 


[M] = DIAG -T (m)^ (in) 

« <»>!} 

(4.4a) 

r (k] 11 -[k) 13 [kj-; ( k] 31 

[ rK] 21 -[K) 23 [K]-*[K] 31 

[ K ) 22 -[ K ] 23 [ K ] 3 1 3 [K] 32 _ 

(4.4b) 


This analytical lumped parameter formulation will be used for example cases I and II. It 
may be noted that if the effect of the rotary inertias is required, equation (4.3) can be 
used instead. 



107 


4. 4.1. 2 System Kinetic Energ y 

A flexible shaft of length L s with rigid rotor of mass m R is subdivided into n s 

number of segments of equal length L e = Ls/n s with lumped masses at both ends of each 

segment as shown in Figures 4.3a and 4.3b for a single-spur gear mesh geared drive 

and single-stage rotor system respectively: Each lumped mass has 3 translational and 3 

rotational degrees of freedom as illustrated in Figure 4.3c. The total system kinetic 

T 

energy Ej is given by {q} [M] {q}. The mass matrix [M] is diagonal and consists 
of the lumped masses and inertias. 

4.4. 1.3 Shaft Stiffness Matrix 

e 

Using the direct stiffness approach, the stiffness matrix [K] of dimension 12 

s 

c 

corresponding to the alternating displacement vector (q(t) } ^ = (u x j(t), u x j + j(t), u y j(t), 

u yj+l(0» tt z j(t), ti z j+i (t), 0 yj (t) , ® y j+i(0> QxjW, ®xj+l(0. 0zj(t)» ©zj+lW ^ > of 4 generic 
shaft segment is given by 

m*. 


t K „/ s tK„ e i' 
I K e„l' [Keel'd 


(4.5) 


6 6 6 "A Q 

and the non-zero elements of [K uu ] s , [ K u Q ] s = [ K Q u ] s and [K QQ ] s of dimension 
6 are given by the following equations (4.6), (4.7) and (4.8) respectively. 


3 

k ll = k 22 = k 33 = k 44 = -k 12 = -k 21 = ~ k 34 = -k 43 = 12 EI / L e .*» 
k 55 = ^66 = — k 56 = _k 65 = ^ E / E e 


(4.6) 
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(a) 

Rotor 



Figure 4.3 Lumped parameter model of the flexible shafts with rigid rotors. Each shaft 
is supported by 2 identical bearings on both ends, (a) Internal rotating 
system of a typical geared rotor system of Figure 4.1. (b) Internal rotating 
system of a typical single-stage rotor system of Figure 4.2. (c) 

Discretization of a shaft segment and the degrees of freedom associated with 
each lumped mass. 
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2 

k l 1 = -k 22 = k 33 = -k 44 = k 12 = _k 21 = k 34 = -k 43 = 6. El/ L e (4.7) 

k ll =k 22 =sk 33 = k 44 = 4El/L e I 

k 12 = k 21 = k 34 = k43 = 2 EI Ae i 

k 55 ~ ^ = “ k 56 = _k 65 = GJ/L e (4.8) 


where E is the modulus of elasticity, G is the shear modulus of elasticity, I is the 
moment of inertia, and J is the polar moment of inertia of the shaft. The lumped 
stiffness matrix [K] s corresponding to {q(t)} sa = {u xj (t), u y j(t), u z j(t), 0 yj (t), 0 x j(t), 

m g 

0 ZJ (t)} , j = 1,2,3,... for the shaft is constructed by the superposition of all [K] 
matrices and merging terms associated with each degree of freedom. 



[t K uJs 

L [K eJs IK ee ] s J 


(4.9) 


The non-zero elements of [K uu ] s , [K uQ ] s = [K 0u ]^ and [K 00 ] s of dimension 3(n s +l) 
are similarly given by the following equations (4.10), (4.11) and (4.12) respectively. 


k ll = k 22 = k 3n s +l,3n s +l = k 3n s +2,3n s +2 = ^ 2 E V E e * 

k 33 = k 3n s +l,3n s +l = AE/L e (4.10a) 

k wj+3,wj = k wj,wj+3 = ~ l 2 E V E e » W I = 1 ,2,4,5, 7,8,.. .,3n s —l (4.10b) 

k wjj,wjj = 24 El/L e ; Wjj = 4,5,7,8,10,1 l,...,3n s — 1 (4.10c) 

k w m +3,w m — k wm,wnr+3 — — ^ E / E e * ^III — 3,6,9,...,3ng (4. lOd) 

^wjy.wjy = 2 ^ E / E e » W IV = 6,8,12,...,3n s (4. lOe) 
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2 

kn = k 2 2 = - k 3n s +l,3n s +l = “ k 3n s +2,3n s +2 = 6 EI / L e (4.11a) 

~ k wj+3,wi = k wj,wpt-3 = ^ El/L e ; wj = l,2,4,5,7,8,...,3n s — 1 (4.11b) 

k l 1 = k 22 = k 3n s +l,3n s +l = k 3ns+23n s +2 = 4 EI A^e * 

k 33 = k 3n s +l,3n s +l = GJ/L e (4.12a) 

k wj+3,wj ~ k wj,wj+3 = ^ El/L e ; Wj = l,2,4,5,7,8,...,3n s — 1 (4.12b) 

k Wn .w n = 8 El/L e ; w„ = 4,5, 7, 8, 10, 11 3n s -l (4.12c) 

kw m+3,wm = k w ra ,w m +3 = ~ GJ/L e ; wjjj = 3,6,9,...,3n s (4.12d) 

k wjv, w iv = ^ » w iv = 6,8,1 2,. ..,3n s (4.12e) 


Note that here [K] s is only for a shaft, of dimension 6(n s +l). The same formulation 
should be applicable to a geared system with multiple shafts where each shaft stiffness 
matrix is still given by [K] s . Accordingly, the shaft potential energy Eu s is given as 

E ur2t‘>iL [K U‘i> Sil < 4 - 13 > 


Formulations for the motor and load rotary inertias, and flexible coupling torsional 
stiffness will be given in the later sections on example cases. 

4.4. 1.4 Gear Mesh Stiffness Matrix 

For the generic geared rotor system with a spur gear pair, the driving and driven 
shafts are coupled via a linear, time-invariant spur gear mesh stiffness k h * kj,(t). Now 
we define the gear mesh coupling stiffness matrix [K] h as follows where coupled 


Ill 


torsional and translational motions of the gear and pinion are defined by {uyQ a (t), 

UyPa(t), 9zGa(0» ®zPaW ) 


k. 


- k. 


k h d G 


k h d P 


- k, 


[K]„ = 


k h d G 


k h d G 


k h d G 


k h d G 


k h d P 


k h d P d G 


(4.14) 


k h d P 


k h d P 


k h d P d G 


Vp 


where (Iq and dp are the gear and pinion diameters respectively; subscript G and P refer 
to gear and pinion respectively. Here, the y axis is taken to be parallel to the gear mesh 
force which is along the pressure angle direction at the mesh point. The potential energy 
E Uh in this case is 


u 


Uh 2{ U yGa’ u yPa’ G zGa ,0 zPa} [Kl h i 


u 


yGa 
yPa 
zGa 
zPa . 


(4.15) 


4.4. 1.5 Flexible Mount Stiffness Matrix 

The flexible mounts are represented by a diagonal stiffness matrix [K] v 
corresponding to {q(t)} ca which consist of effective stiffness coefficients k vw , w = x, y, 
z, 0 X , 0 y , 0 Z . This modeling procedure assumes that no coupling exists between the 
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casing rigid body degrees of freedom due to the flexibility of the mounts. Accordingly, 
one has 


E Uv=f{1>I tKl v{‘l> <a (4 - 16) 

4.4.1. 6 Bearing Stiffness Matrix 

The proposed bearing stiffness matrix [K] bm corresponding to the bearing degrees 
of freedom {S wa (t),f} wa (t)} T , w = x, y, z, has been formulated in Chapter II. The 
corresponding potential energy Ey^ is 

Ein.= T«T5 .,(3 .T [K] 

Ub 2 1 waj’ K wajj bm 


waj 


waj 


(4.17) 


Here, {5 wa j(t),|3 wa j(t)} T may be expressed in terms of {qC0)saj (q(O)ca through a 
coordinate transformation for the j-th bearing located at R j = {xj, yj, zj) from the 
casing center of mass. For the geared rotor system, four bearings are denoted by j = 1, 
n s +l, n s +2 and 2n s +2 as shown in Figure 4.3a, while for the single-stage rotor system 
two bearings are denoted by j = 1 and n s +l as shown in Figure 4.3b. Assuming the 
small angle approximation, {5 wa j(t),(3 wa j(t)} T for the j-th bearing is 
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(4.18) 
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4.4.2 Method B: Dynamic Finite Element Formulation 

The dynamic finite element technique of incorporating [K]t> m in equation (4.1) is 
employed especially when the casing plate is elastically deformable over the frequency 
range of interest, as discussed previously in Chapter III. In this method, [K] bm is 
implemented in the dynamic finite element model as a generalized stiffness matrix. This 
stiffness element provides the only coupling between the flexible shaft and casing at the 
appropriate nodal points on both components corresponding to the bearing locations. 
Additionally, the gear mesh stiffness kj, in Figure 4.1 behaves like a linear translational 
spring. Since the coupling between torsional and transverse shaft vibrations is due to k h 
only, the corresponding portion of the stiffness matrix is also given by equation (4.14). 
Our finite element formulation uses conventional structural elements typically available 
in commercial software programs [39] for the shaft, casing plate and mounts - this will 
be illustrated for example cases I, II and HI. Other features of this method are similar to 
those discussed earlier in Section 4.4.1. 

4.4.3 Other Methods 

Alternate techniques of formulating the geared rotor system dynamic problems 
such as finite difference (which is similar to method B), flexibility, component mode 
synthesis and transfer matrix formulations are also possible. It may be noted that 
Berman [54] has philosophically proposed application of the component mode synthesis 
to geared problems. Other researchers [40,41] have utilized the transfer matrix method 
in the rotor dynamic problems. However, such models need to be modified or extended 
in order to solve the overall dynamic problems discussed in this dissertation. Direct 
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application of these alternate methods [40-43,54] are beyond the scope of this 
dissertation and are left for further research. 

4.5 EXAMPLE CASE I: SINGLE-STAGE ROTOR SYSTEM WITH 
RIGID CASING AND FLEXIBLE MOUNTS 

4.5.1 Vibration Models 

Consider the single-stage rotor system shown in Figure 4.2 where two axially 
preloaded identical deep groove ball bearings support a flexible rotating shaft; see Table 

4. 1 for bearing parameters. The shaft is assumed to be subjected to a mean torque T^n, 
* T zsm (t) and harmonic excitation forces (F xsa (t) = F^^OJ e^* 1 , F ysa (t) = F ysal (Q z ) 
e i(£l z t-K/ 2 ))T j ue tQ rotat j n g rotor mass unbalance. The rotor is fixed to the center of 
the shaft which coincides with the rigid casing centroid. The shaft is assumed to be 
transversely decoupled from the motor and load due to the flexible torsional couplings. 

A lumped parameter model with n s +2 lumped masses is developed according to 

the theory given in Section 4.4.1. The bearing stiffness matrix for an axially preloaded 
ball bearing has significant stiffness coefficients k bxx , k byy , k bzz , k b0x0x , k b0y0y , k bx0y 

and k by0x which depend on the mean axial preload. Since there is no external Fza(t) 

force on the system and no coupling between the torsional and transverse motions of the 
shaft, the dynamics associated with { u zsa j(t) ,u zca (t),0 zRa (t) } T degrees of freedom (here 
subscript j is a dummy index to identify the shaft lumped masses) are decoupled from 
the others and have trivial steady-state particular solution. The system matrices of 
equation (4) can now be simplified by suppressing such degrees of freedom. Also, it 
can be easily observed from equations (1) and (4) that 2 sets of uncoupled differential 
equations exist. One set has {q(t)} a = {u xsa j(t),u xca (t),0 yRa (t),9 yca (t)} T which are 
excited by F xsa (t) and the other set consists of {q(t)} a = {u ysa j(t), u yca (t), e xRa (t), 


Table 4. 1 Design parameters for deep groove ball bearings used for example cases t 


Load-deflection exponent n 

3/2 

Load-deflection constant K n (N/m n ) 

7.5 E9 

Number of rolling element Z 

12 

Radial clearance r L (mm) 

0.005 - 0.02 

Pitch diameter (mm) 

38.5 

Aq (mm) tt 

0.0625 

Unloaded contact angle a,. 

23° -47° 


t These bearings are currently being used in a NASA gear test facility [55]. Also see 
example cases II and HI. 

tt Unloaded distance between inner and outer raceway groove curvature centers. 
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flxcaG) } T with Fysa(t). The steady-state solution to these two sets will always result in 
the synchronous whirling motion of the internal rotating shaft However, since the two 
sets of equations are similar and independent only one set is considered in the following 
analysis. Hence, the problem reduces to a 2 dimensional vibration system with DOF = 
n s +4 which consists of a flexible shaft with rigid rotor vibrating transversely and 
coupled to a rigid casing and a flexible mounts through two bearings. Accordingly, [M] 
in equation (4.4) is rewritten in terms of the alternating displacement vector { q(t) } a = 
{ u xsaj(0> u xca(t)’®yRaW>®yca(t)) as 

[Ml = DIAG {m se .,m c ,I yR ,I yc } (4.19) 

where m and I y are the lumped mass and inertia about the y axis respectively, and the 

subscript e refers to the shaft segment. To define [K] in equation (4.4), the nonzero 
elements of [K] n , [K] 12 , [K] 2 2 , [K] 23 , [K] 13 and [K] 33 are given by equations (4.20), 

(4.21), (4.22), (4.23), (4.24) and (4.25) respectively. 


2 

k ll = k n s +l,n s +l = 12 EI/L e + kjj XX ; k ns+ 2^is+2 = k vx + 2 k bxx > 

k n s +2,l = k n s +2,n s +l = k l,n s +2 = k n s +l,n s +2 = ~ k bxx (4.20a) 

kw I+ l.wi = k w I ,w I+ l=-12El/L e ; wj = 1,2,3, ...,n s (4.20b) 

k wjj,wjj = 24 El/L e ; wjj = 2,3,. ..,n s (4.20c) 

2 

k ns/2,l = — k ns/2+2,l = ^ El/L e , k ns+22 = ( z l + z n s +l) k bxx 2 k bx0y » 

k 12 = -Z 1 k bxx — k bx9 y » k n s +l,2 = “ z n s +l k bxx — k bx0 y (4.21) 


117 

k 11 = 8 El/L e ; k 22 = k v0y +2 k b0y0y + ( z\ + z J s+1 )k bxx + (z^z,^) k bxQy (4.22) 

kl,ns/2 = ki, ns /2+l = 2 El/Lg ; 

k2l = ~ z i kbx0y ” k b0y0y ; k 2jI , s = - z„ s+ | k bx0y ~ khQ y0y (4.23) 

2 

kl 1 = ~kn s +l,n s = 6 El/L e + k bx0y ; k ns+2 j = k ns+2> „ s = -k bx0y ; 

2 

kns/2+1, ns/2+1 = = 6 El/L e (4.24a) 

2 

k i2 = -k ns +i,ns-l = -kns/2, ns/2-l =k ns /2+2jis/2+2 = 6 El/L e ; n s >2 (4.24b) 

2 

kwjjj-l.wjjj ~~k Wm -i iWm -2 =- k Wm+ns /2, Wm _ 2+ns /2 = k Wm+ns /2 Wni+n ^ = 6 El/L e ; 


n s >6, wni = 3,4,...,ns/2 (4.24c) 

ki 1 = k nstns = 4 El/L e + k^Qy (4.25a) 

kw IV ,w IV = 8 El/L e ; n s >2, wjy = 2,3,... ,n s — 1 (4.25b) 

kwy.wy— 1 — 2 El/L,, , k w ^_^ w ^ = 2 El/Lg ; 

n s >2, w v = 2,3,...,n s (w v * n s /2+i) (4.25c) 


Due to the symmetric nature of [K], [K] . = [K] . T . The generalized load vector 

wj jw 

{f(t)) a = (0 FxsaW’ 0.— } consists of only the x component of the rotor mass 

unbalance force. The other set of linear governing equation may be obtained by 
interchanging subscripts x and y in the above formulation. Our proposed vibration 
model can also be readily reduced to the conventional rotor dynamic models, excluding 
the gyroscopic moment, internal structural damping, dissimilar bearings and other 
secondary effects, by retaining only k bxx , k byy and k bzz in our bearing formulation. 
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4.5.2 Eigensolution 

The natural frequencies C 0 j and modes <j>j of Figure 4.2 with system parameters 
given in Table 4.2 are obtained numerically. Initially, we check for coj convergence rate 
as shown in Figure 4.4. through 034 and 035 by varying the number of shaft segments 
n s ; it may be noted that tDj, o>2 and CO3 behave like CO4, and cog trend is very similar to 
005 . Beyond n s = 4, there are no noticeable changes in the natural frequencies. These 
predictions are now compared with two simple lumped parameter vibration models and a 
dynamic finite element model (FEM) in Table 4.3. One of the simple models excludes 
the effect of casing and flexible mounts, and both employ a conventional bearing model 
with only k bxx incorporated as described earlier in Section 4.5.1. The dynamic finite 
element formulation models the flexible casing constructed with four noded quadrilateral 
plate elements - a description of the methodology will be presented later in example case 
III. Our proposed theory and FEM are in excellent agreement with each other as the 
deviations are within ±4% for first 6 natural frequencies. The corresponding mode 
shapes <)>j are shown in Figure 4.5. Note that for each mode in the x-z plane given here, 

there is a complementary mode in the y-z plane which will have the same natural 
frequency if k vy =k vx , k v Q x =k v Q y and I xc =I yc . Using the exactly same parameters, 

results obtained from two simple models are also compared in Table 4.3. The simple 
model with casing and mount dynamics included predicts lower natural frequencies as 
compared to our prediction except for co 5 ; here the symbol a implies estimation using a 
simple formulation. This model gives a rough approximation for most of the modes, as 
evident from Table 4.3. Further, the simple model without casing and mounts can not 
predict all the modes below 2 kHz. In addition, CO 3 shifts while other modes remain 
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Table 4.2 System parameters of example case I: single-stage rotor system shown in 
Figure 4.2 


Rotor mass m^ (kg) and inertia I y R (kgm 2 ) 

5.0, 1.24E-2 

Shaft mass (kg) and length L s (m) 

1.03, 4.2E-1 

Casing mass (kg) and inertia I yc (kgm 2 ) 

148, 4.27 

Shaft flexural rigidity El (Nm 2 ) 

1.626E3 

Bearing axial preload F z b m (N) 

238 

^bxx» ^bx0 y » ^b0 y 0 y ^ 

1.37E8, -1.96E6, 4.27E4 

k vx » k v Qy 

4.0E8, 1.8E7 


t Computed kj, w j as proposed in Chapter II. Only relevant bearing coefficients are 
listed. 


Table 4.3 Comparison of natural frequencies in Hz for example case I 



Proposed 

Theory 



Simple Models 


Mode 

FEM 

Aj% 

without casing 

with casing 

a 2 % 

1 

97.2 

97.2 

0.0 

69.5 

69.4 

28.6 

2 

262 

262 

0.0 

not predicted 

262 

0.0 

3 

309 

309 

0.0 

297 

294 

4.9 

4 

332 

334 

0.6 

. not predicted 

329 

1.5 

5 

1405 

1365 

-2.9 

1442 

1443 

-5.7 

6 

1755 

1690 

-3.8 

1452 

1453 

14.0 


A i% = 100 x (FEM - Proposed Theory) / FEM 
A 2 % = 100 x (FEM - Simple Model) / FEM 


Natural Frequency (kHz) 
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nearly the same. These observations are consistent with the findings reported earlier 
[51-53]. 

We examine the effect of casing mass (and inertia I yc « m^.) while retaining 
other parameters of Table 4.2. The fundamental frequency coj (first transverse mode of 
the shaft-rotor system) and 0)3 (second transverse shaft-rotor mode) are quite insensitive 
to me. On the other hand, 0)5 and (Og are proportional to (rric)' 1 ^ for a light-weight 
casing given by m c /(m s +m R ) < 1 , but are invariant for m c /(m s +m R )> 1 as shown in 
Figure 4.6a. The converse is seen for CO 2 and 0)4 which are nearly constant for 
m c /(m s +m R ) < 1, but are proportional to (m c )' 1/2 for m c /(m s +m R ) > 1. These results 
imply that the single degree of freedom approximation can be applied to 0)5 and co 6 in 

the region where m c /(m s +m R )< 1 , and similarly to CO 2 and CO 4 in the other region. 

2 

Next, we vary mount stiffness coefficients k vx (and = L s k vx ) to simulate 
the effect of mount flexibility on the system. Here 0)5 and cog are found to be unaffected 
by k vx due to the fact that the corresponding modes are predominantly shaft bending 
motion type as shown in Figure 4.5. Figure 4.6b indicates that 0)i and 0)3 are 
proportional to (k vx ) 1/2 for soft mounts (Kx/^-bxx « 1 ) but are insensitive to k vx for 
stiff mounts (k vx /kg xx » 1). Conversely, 0)2 and 0)4 are proportional to (k ^) 1 ^ 2 for 
stiff mounts but remain nearly invariant for soft mounts. Similarly, each mode in its 
appropriate region can be assumed to behave like a single degree of freedom system. In 
Figure 4.6b, unlike Figure 4.6a where only one transition point at m c =m s +m R exists, 2 
transition points are found. These are located below and above k vx /k bxx = 1, and 
increase with increasing shaft and bearing stiffnesses. 

Now, we examine the range of bearing preload F 2 5 m from 100 to 5000N on coj. It 
is seen that coj increases with increasing F zbm due to an increase in the magnitudes of the 


Natural Frequency (Hr) Natural Frequency (kHz) 
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(a) 



Figure 4.6 Hffect of casing mass m^ normalized with respect to shaft and rotor masses, 
and mount stiffness k vx normalized with respect to k bxx on system natural 
frequencies C0j for example case I. 
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bearing stiffness coefficients k bwj - as shown in Table 4.4 for the two extreme F zbm 
values; this observation is similar to those found earlier in the bearing system studies of 
Chapter III. Note that ©2 * s most insensitive to F zbm , and ©3 and ©4 show a 
moderate variation. But ©1 increases by about 45% due to larger bearing motions as 
compared to ©2, ©3 and ©4 modes. The remaining two modes of interest are also 
affected significantly as large motions across the bearings are again found in these 
modes. Unlike the trends associated with variations of m,, and k vx , here each ©j 
approaches an upper bound solution, as shown in Figure 4.7 for © lt if F zbm is 
sufficiendy high to stabilize magnitudes of k bw j terms. 


Table 4.4 Effect of bearing preload F zbtn on ©j (Hz) for example case I 


Preload F zbm (N) 

©l 

©2 

©3 

©4 

®5 

“6 

100 

89.4 

262 

305 

331 

1366 

1706 

5000 

130 

263 

321 

346 

1855 

1968 
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4.5.3 Transmissibility 

The forced harmonic response due to the rotating rotor mass unbalance force 
F xsa (t) = F xsa i(fi z ) e iflzt is obtained using the dynamic stiffness technique given earlier 
in Chapter III. Also, recall from Chapter III the definition of sinusoidal load 
transmissibility functions R(Q Z ) between two arbitrary locations I and II with load 
vectors f w i a (t) and fjn a (t) as R(Q Z ) = l f w.a< n 2 )|/l fjli a (& z )| where in this case 

fjila(t) = Fxsa(t). The bearing force RF xba ,F xsa (^z) and moment RM yba ,F xsa (O z ) 

transmissibility spectra are compared in Figure 4.8 for our formulation and simple 
model, given the same system of Table 4.2, both models differ only in [K] bm 

formulation. Note that the results presented are for only one bearing due to symmetry of 
the problem. We observe that the simple model predicts lower ^F xba ,F xsa (^z) ^an 

proposed model except in the vicinity of cop On the other hand, ^M yba ,F xsa (^z) is not 

predicted at all by the simple model. This component of the bearing transmissibility 
spectra is primarily due to coefficients k bx Q y and k b Q y Q y which are obviously not 

included in the simple model. Our predictions also show that higher modes over 250- 
350 Hz contribute significantly to the bearing transmissibility which are not seen in the 
force transmissibility spectra yielded by the simple model. 

Now we compare the mount load transmissibility in Figure 4.9. We observe that 
the trends of mount force transmissibility Rp p (£2 Z ) spectra predicted by both 

models are similar, although the simple model predicts slightly lower amplitudes above 
cop The simple model again can not predict the moment transmission RM yva ,F xsa (^z) 

across the mounts. Unlike the bearing transmissibility spectra, here the effects of $3 
and <J ) 4 on RM yva ,P xsa (^z) are more significant than <>4 and <j> 2 . This is mainly due to 

the large angular motions on the casing for <{>3 and $ 4 . 


Moment Transmissibility (m) Foice Transmissibility 
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Figure 4.8 Bearing transmissibility for example case I. (a) Force transmissibility 
RF xba> F xsa (Qz) - (b) Moment transmissibility R M ba ,F xsa (Q z ) which is not 
predicted by simple model. y a XSa 



Moment Transmissibility (m) Force Transmissibility 
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(a) 



(b) 


Figure 4.9 Mount transmissibility for example case I. (a) Force transmissibility 
R Fxva>Fxsa (Qz) - (b) Moment transmissibility R Myva ,F xsa (^ z ). 





Mount Transmissibility Bearing Transmissibility 



10 J 
10 * 

10 ” 

10 ' 

10 l 
10’ 

10 ^ 

10 5 

0 250 500 

Frequency (Hz) 

(b) 

Figure 4.10 Effect of mount stiffness k vx on the force transmissibilities for example 
case I. (a) Bearing, (b) Mount. 
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Next, we investigate the effects of rr^ and k vx on the bearing and mount 
transmissibilities. Figure 4.10a shows that essentially the same bearing transmissibilities 
are obtained for three k vx values. In the case of the mount transmissibility as shown in 
Figure 4. 10b, the amplitudes in general decrease with lower k^ except near resonances. 
The casing mass n^ also does not affect the bearing transmissibility significantly as 
shown in Figure 4.1 la. The mount transmissibility amplitudes increase with lower me 
as shown in Figure 4.11b except near resonances again. It may be noted that similar 
trends are observed in the bearing and mount moment transmissibility spectra. 

t 

4.6 EXAMPLE CASE II: GEARED ROTOR SYSTEM WITH RIGID 
CASING AND FLEXIBLE MOUNTS 

4.6.1 Bearing Analysis 

Now consider the geared rotor system in Figure 4.1 with flexibly mounted rigid 
casing whose bearing and other system parameters are given in Tables 4.1 and 4.5 
respectively. Four ball bearings which support the two shafts are subjected to identical 
mean axial displacement 5^. The spur gear pair drive is driven by a mean torque T^ 
* T zsm (t) which also generates mean radial bearing force F ybm j and moment M xbm j. 
The stiffness matrix [K] bm for each bearing under these loads has significant 
coefficients k bxx , k byy , k bzz , k b Q x Q x , k b g y 0 y , k bx Q y , k by 0 x , k byz and k bz Q x< A set of 

governing nonlinear algebraic equations, consisting of 3 bearing load-displacement 
relations for each bearing as given in Chapter II and from the shaft bending theory is 
given by 


Mxbmwj Tzsin^s/O? d(j) 2EI(2P xmWl +P xmWn +3(5y mWl 8 ymwjj)/Ls)/L s ; 
M xbmw n =- M xbm Wl ; for w II =n s+l if w P l, and w n =2n s +2 if wpn s +2. 


(4.26) 
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Table 4.5 


System parameters of example case II: geared rotor system shown in Figure 
4.1 + 


Gear and pinion masses mp = mg (kg) 

0.5 

Gear and pinion rotary inertias lyG^yP^zG^zP (kgm 2 ) 

1.5E-4, 3.0E-4 

Number of gear/pinion teeth 

28 

Gear/Pinion Pressure Angle (degrees) 

20° 

Shaft mass m s (kg) and length L s (m) 

2.8, 2.54E-1 

Casing mass m*. (kg) rotary inertias ly^I^ (kgm 2 ) 

77, 1.5, 1.9 

Shaft flexural rigidity El (Nm 2 ) 

1.25E4 

Effective Torsional Stiffness kj^k-pj^Nm/rad) 

6.05E3, 2.45E4 

Motor and Load inertia I^ 

1.00E-1, 3.35 

Gear mesh stiffness kjj (N/m) 

1.0E8 

Mean axial bearing displacement 5^ (N) 

6.00E-4 

Mean input torque T zbm (Nm) 

72 

k byy> k by0 x > k b0 x 0 x 

9.7E8, 6.0E5, 1.0E6 

k vy* k v0 x » k v0 z 

9.0E8, 1.4E7, 2.4E7 


t The internal rotating spur gear pair is currently being used in a NASA gear test 
facility [55]. 

tt Only relevant bearing stiffness coefficients are tabulated. 
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Due to the physical symmetry, we assume 8 ymWl = 8 ymWn and p xmWl = *Pxmw n which 
simphfy the nonlinear algebraic equation set for each bearing to 

M xbm = T zsm L s/ ( 8d G)“ 2EI Pxm/ L s 5 F ybm = T 7sml<k} ( 4 - 27 ) 

where M x b m and F bym are expressed in terms of a constant 8 zm , and variables 8 ym and 
Pxm- Solution to these 2 equations may then be used to compute kbjj(8 zm ,8 ym ,P xm ) 
directly. 

4.6.2 Vibration Models 

A lumped parameter model of Figure 4. 1 is developed in accordance with Section 
4.4.1. Each shaft is divided into n s number of segments as shown in Figure 4.3a. The 
system matrices of equation (4.4) are simplified by neglecting the longitudinal shaft 
motion and casing rigid body degrees of freedom u zca (t), u xca (t) and 0 yca (t). Both 
driving and driven shafts are assumed to be coupled to the motor and load respectively 
through flexible torsional couplings. Therefore, only the motor and load 1^ rotary 
inertias are considered, and the exterior portion of the shafts beyond the flexible 
couplings are modeled as purely torsional stiffness elements which are then combined 
with the flexible coupling stiffnesses. The system is excited by the static transmission 
error e(t) at the mesh point which generates gear and pinion force F h (t) = k h e(t), parallel 
to the line of action, and torque T^t) = l/2d(}k h e(t) or T^t) = l/2d P k h e(t), about the 
axial z-axis. The mass matrix [M] in equation (4.4a) in terms of (q(t)} = (u ysa j(t), 
u ycaW> ^xGaW’ ®zGaW’ ®xPa^’ ®zPa^* ®zMa^’ ®zLaW» ®xca^)’ ®zca^)) is 
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IM] — DIAG { m sej’ m xc>^xG>^zG»^xP>IzP^zM>^zL>^xc>^zc} (4.28) 

The matrices [K]jj, [K] 12 , [K] 13> [K] 22 , [K] 23 and [K] 33 define the symmetric [K] of 
dimension 2n s + 1 1 in equation (4.4b). Nonzero elements of [K]^ of dimension 2n s +3 
are 


^n s +2,n s +2 = k 2n s +2,2n s +2 = k ll = k n s +l,n s +l = 12 El/ L e + kj>yy 

3 

^wj+l.wj = ^wj.wj+ 1 = “12 El/L e ; 

3 

^ c w I +n s +2,wj+n s +l = ^wj+n s +l,wj+n s +2 = — ^ El/L e ; Wj = 1,2,3,.. .,2n s 
k w n , w n = 24 El / L e ; k Wn+ns+ 1 iWn+ n s +i = 24 El/L e ; 


wji = 2,3,. ,.,n s /2,n s /2+2,...,n s 

3 3 

k ns/2+Ms/2+l = 24El/L e +k h ; k 3ns/2+2i3ns/2+2 = 24 El/L c + k h ; 

k ns/2+l, 3 ns/2+2 — k 3ns/2+2,ns/2+l = —k h » k 2n s +3,2n s +3 = 4k byy + k vy » 
k w m ,2n s +3 k 2n s +3,wjjj k byy > — 1» 2n§+2 


(4.29a) 

(4.29b) 

(4.29c) 

(4.29d) 

(4.29e) 


The nonzero elements of [K] 12 of 2n s +3 rows and 8 columns are 

2 

k 3n s /2-l,l=- k 3n s /2+3,l = k n s /2,3 = - k n s /2+2 I 3 = 6EI / L e * k 2n s +3,7 = 4k by6 x i 
~ k ns/2+l,2 = k 3ns/2+2,2 = k h/< 2 ^g) i k ns /2+l,4 = - k 3ns/2+2,4 = k h/(2d P ) (4.30a) 
kw in> 7= “ k by e x +z w m k byy » k wjjj,8 =— x wqj k byy ’> wm=l,n s +l,n s +2,2n s +2 (4.30b) 

The nonzero elements of [K] 13 of 2n s +3 rows and 2n s columns are 
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2 2 
k l 1 = k n s +2,n s +l = 6 El/L e + kby6 x » k n s +l,n s = k 2n s +2,2n s = “ ^ El/L e + kby0 x 5 

2 

k 3ns/2+2,3ns/2+l = -k 3iis/2+l,3ns/2 = k ns/2+l,ns/2+l = “ k ns/2+Ms/2 = 6 EI / L e ; 

k 2n s +3,l = k 2n s +3,n s = k 2n s +3,n s +l = k 2n s +3,2n s = -k by9 x (4.31a) 

2 

k 12 = —k n s +l,n s -l = ^ 5 / 2 , 118 / 2-1 = k ns/2+2,ns/2+2 = k n s +2^is+2 = 6 El/L e ; 

2 

_k 2ns+2,2ns-l = - k 3ns/2+l,3ns/2-l = k 3ns/2+3,3ns/2+2 = 6 EI / L e ; n s >2 (4.31b) 

2 

_k W IV , w i v~ 1 =k w IV ,w IV + 1 =k Wiv+ n s/2+ 1 ,wjy 115 / 2 + 1 =-k wjy+n s +l ,Wjy- l+n s = 6 El/ L e ; 

2 

k w I y+ns+l,wiY+l+ns =-k wiv+ns/2+l,w I v+ns/2-l =k w I y+3n s /2+2,wiy+3n s /2+l := ^ E V E e * 

2 

■ k wjy+3ns/2+2,wiv+3ns/2-l — 6 E V E e * ^s— ^ * ^IV = 2,3,...,n s /2— 1 (4.31c) 

where nonzero elements given by equations (4.31b) and (4.31c) exist only if n s >2 and 
n s >6 respectively. The nonzero elements of symmetric [K ] 22 of dimension 8 are 

k n = k 33 = 8 El/L e ; k 22 = k TM + k h d G /4 ; k 25 = -k TM ; k 4 6 ■ -k^ ; 

k 44 = k TL + k h d p/4 ; k 2 4 = k42 = -k h d G d P /4 ; k 55 =k TM ; k^ = k^ ; 
k 77 = 4k b0 x 6 x + E (Zj) 2 k byy + k v0 x ; k 88 = k v0 z + £ (xj) 2 k byy (4.32) 

where k^L and ky^ are the effective torsional stiffnesses due to shaft and torsional 
coupling respectively at load and motor ends, and the summation X is over the 4 ball 
bearings denoted by j = 1 , n s +l, n s +2, 2n s +2. The nonzero elements of [K ] 2 3 of 8 
rows and 2 n s columns are 

k Ms/2 = k l, ns/2+1 = k 3,3ns/2 = k 3,3ns/2+l = 2 EI Ae 
k 7,w v = + z w m k by0 x “ k b0 x 0 x » k 8,w v = -x w ffl k by0 x I 


(4.33a) 


136 


w IH=l(wv=l),n s +l(w v =n s ),n s +2(wv=n s +l),2n s +2(wv=2n s ) (4.33b) 

Finally, the nonzero elements of symmetric [K ]33 of dimension 2n s are 

^11 = ^n s ,n s = ^n s +M s +l = ^ 2 n s , 2 n s = 4 ^A J e + ^b 0 y 0 y (4.34a) 

^wyj.wyj = ^n s +wyj,n s +wyj = ^ El/L e ; n s >2, wyj = 2,3,...,n s (4.34b) 

^■Wyn.Wyjj-l = ^Wyjj-l.Wyn = ^ns+Wyn^s+Wyn-1 = 2 El/Lg ', 

k ns+ w VII -Ms + w vn =2El/L e ; n s >2, wyjj = 2,3,...,n s (wyjj ^ n s /2+l) (4.34c) 


As discussed earlier, our proposed vibration model again differs from the conventional 
spur gear pair dynamic models essentially due to [K] bm formulation. The simple model 
with casing and mount dynamics can be obtained from our proposed model by retaining 
only k bxx , k b yy and k bzz in [K] bm . 


4.6.3 Eigensolution 

The natural frequencies toj and modes <j)j of Figure 4.1 with system parameters 
given in Tables 4.1 and 4.5 are computed using the proposed theory, a dynamic finite 
element model, and the simple theory with and without casing and mounts. The FEM 
model includes casing flexibility - this will be described later in example case II. Our 
proposed theory differs from FEM by less than +4% for the first six ©j and is within 
+ 10% for ©7 and ©g as shown in Table 4.6a. The corresponding natural modes are 
described in Table 4.6b. Simple models in general predict lower natural frequencies and 
deviate substantially except for ©j and © 3 . The simple model without any casing 
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Table 4.6 Eigensolution of example case II 


a. Comparison of natural frequencies coj in Hz 


Mode (j) 

Proposed 
j Theory 

FEM 

Aj% 

Simple Models 

. * 

without casing with casing 

a 2 % 

1 

0 

0 

0.0 

0 

0 

0.0 

2 

31.0 

31.0 

0.0 

29.8 

29.1 

6.1 

3 

480 

480 

0.0 

not predicted 

480 

0.0 

4 

526 

525 

-0.2 

not predicted 

510 

2.9 

5 

563 

564 

0.2 

not predicted 

560 

0.7 

6 

790 

760 

-3.9 

666 

650 

14.5 

7 

1147 

1093 

-4.9 

837 

809 

26.0 

8 

1208 

1100 

-9.8 

694 

681 

• 38.1 

A j % = 100 x (FEM 
A 2 % = 100 x (FEM 

- Proposed Theory) / FEM 

- Simple Model) / FEM 




b. Summary of mode shapes 

Mode <j>j 

Description 

1 

Motor-gear-pinion-load system rigid body torsional motion 

2 

Motor-gear-pinion-load system torsional motion 

3 

Casing rigid body rotational motion 0 xca 

4 

Casing rigid body translational motion u yca 

5 

Casing rigid body torsional motion 6 7ra 

6 

First shaft coupled transverse-torsional motion 

7 

Second shaft coupled transverse-torsional motion 

8 

Shaft transverse motion 
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dynamics as outlined here is in excellent agreement with the results reported by 
Kahraman et al. [ 15 ]. This provides further validation for our theory. 

The effect of casing mass (and inertias I xc , 1 ^ « n^) is shown in Figure 4 . 12 a 
for CO3, CO5 and 0)7. Here CO2 an< 3 (Og are unaffected by rr^. On the other hand, CO3 is 
always proportional to ( m c )*^ 2 except when nic^ms+mo+mp) « 1 . This 
proportionality feature, similar to the single degree of freedom theory, is also seen for 
CO7 and cog for a lightweight casing, and 0)4 and 0)5 for a massive casing. It may be 
noted that the trends of cog and CD4 are similar to CO7 and 0)5 respectively although they 
are not included in Figure 4 . 12 a. Unlike example case I where the transition point is at 
m c /(m s +ni(3+mp) = 1, two transition points are observed here due to the effects of the 
gear, pinion and casing rotary inertias. 

The effect of mount stiffness k vy (and k v0x , k v0z « k vy ) on coj is given in Figure 
4 . 12 b. We again find that 0)2 and co 6 are not affected by a variation in k vy . Also, two 
transition points are found which separate a region of constant coj from the region where 
C0j * (kyy) 1 / 2 for some of the natural frequencies. Here cog is similar to 007 which is 
proportional to (kyy) 1 ^ 2 for stiff mounts but is insensitive to k vy for softer mounts. The 
converse is seen for CO2 with a lower transition point. On the other hand, 035 « (kyy) 1 / 2 
in the region between the two transition points, and is constant elsewhere. Finally, ©3 
is seen to be always increasing with kyy. 

The effect of mean bearing displacement 8^ is summarized in Table 4.7 for two 
extreme 5 zm values. Only the natural frequencies associated with the shaft transverse 
and/or torsional motion are sensitive to 8 zm or preloads F z {, m (8 zm ). The rate of 
increase for coj is high for low 8 zm , but the rate decreases with a higher 8 zm as observed 
earlier in example case I due to the stabilization of kgjj. 


Natural Frequency (kHz) Natural Frequency (kHz) 
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(b) 

Figure 4.12 Effect of casing mass (normalized with respect to shaft and gear 
masses) and mount stiffness k vx (normalized with respect to k bxx ) on the 
system natural frequencies coj for example case II. 
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Table 4.7 Example case II: effect of mean bearing displacement 8 zm on coj (Hz) 


S z bm (mm) 

0)1 

0) 2 

©3 

©4 

®5 

©6 

©7 

© 8 



29.9 



562 

764 

852 

852 

1.00 




526 

563 

791 

1180 

1215 


4.6.4 Transmissibility Spectra 

The response due to the static transmission error excitation e(t) at the mesh point is 
computed using the same dynamic stiffness technique used for example case I. Here, 
only the fundamental harmonic of e(t) at mesh frequency © b is included. All 

transmissibility functions are normalized with respect to the magnitude of the gear mesh 
force Fj,(t) = kj, e(t). Figure 4.13a and 4.13b compare the bearing force RF yba5 F h ( w h) 

and moment RM xba ,F h ( co h) transmissibility spectra as predicted by our formulation and 

simple model. It may be noted that the transmissibility spectra for all four bearings are 

similar. Although simple model compares reasonably well with our proposed theory for 
R F y ba,Fh (co h). ft ft not capable of predicting RM xba ,F h (° ) h)* as also seen previously. In 

the case of the mount transmissibility shown in Figure 4.13c, only the net moment 
M zva (t) is transmitted. The resultant vertical force F yva (t) and moment M xva (t) are 
negligible due to the force and moment cancellations at the mount feet. The Fourier 
spectrum of the normalized dynamic transmission error p(t) = [y p (t)— 
y G (t)+(dp0p(t)/2)-(d G 9 G (t)/2)]/e(t) is shown in Figure 4.13d. Only two modes, <J> 2 
which is predominantly torsional vibration of the shafts and <t>g which is a coupled 




Mount Transmissibility (m) Force Transmissibility 
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Figure 4.13 Comparison of transmissibility and normalized dynamic transmission error 
spectra for example case II. (a) Bearing force transmissibility 
Rpyba.Fh^li)- Bearing moment transmissibility R M xba ,F h (®h) which 
is not predicted by simple model, (c) Mount transmissibility 
RM zba ,F h ( co h)- (d) Normalized dynamic transmission error p(coj,). 
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Figure 4.14 Mount transmissibility studies for example case II. Effect of (a) casing 
mass me, (b) mount stiffness k vx and (c) mean bearing displacement 5 zm 

on R M zba .Fh (c0 h ) - 








Moment Transmissibility (m) Force Transmissibility 


143 




Figure 4.15 Effect of mean bearing displacement 8 zm on (a) bearing force 
transmissibility R Fyba)Fh (oih)> (b) bearing moment transmissibility 

^M xba ,F h ( co h)> an ^ (c) normalized dynamic transmission error p(co h ) for 
example case II. 
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transverse-torsional mode, affect p(a>h) significandy. The sixth mode produces a higher 
p((%) amplitude at co^ than fa- Additionally, the simple theory underestimates p(cojj) as 
compared to the proposed model. 

The effects of m^., k vx and 8^ on the transmissibility spectra and p((%) are given 
in Figures 4.14 and 4.15. Mount transmissibility RM zvat F h ( (0 h) * s lower for a heavier 

casing, more compliant mount and higher 8 zm as shown in Figure 4.14. In addition, 
the effects of mg and k vx on RM zva ,F h ( co h) are more prominent than those of 8 zm . The 

bearing transmissibility spectra and p(C 0 h) are affected only by ^zm as shown in Figure 
4.15, and are virtually insensitive to mg and k vx for the parameters given in Tables 4.1 
and 4.5. Vibratory force and moment transmission through the bearing are slightly 
lower for a larger 8 zm due to the additional constraint provided by higher kyj values on 
the rotating shafts. Conversely, normalized dynamic transmission error p(cofj) increases 
if a larger bearing preload is specified. 

4.7 EXAMPLE CASE III: GEARED ROTOR SYSTEM WITH RIGIDLY 
MOUNTED FLEXIBLE CASING 

4.7.1 Physical Setup 

The final example case examines the NASA gear test facility as shown in Figure 
4.16 [55]. The system parameters are equivalent to those in Tables 4.1 and 4.5 except 
for the flexible steel casing of approximate dimensions 0.33 x 0.28 x 0.25 m and plate 
thickness of 0.006m. Rigid mount feet attach the four comers of the bottom casing plate 
to a massive foundation. High precision gear and pinion are used which are identical 
with 0.006m facewidth, 0.089m diameter and 1:1 ratio. Four axially preloaded high 
precision deep groove ball bearings are being used to support 0.03m diameter shafts of 
length 0.254m on the flexible casing. The input and output shafts are only coupled 
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Figure 4.16 Example case IE: (a) Schematic of the NASA gear test facility, (b) NASA 
spur gear pair system supported on rigidly mounted flexible casing [55]. 
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torsionally to the rest of the gear test facility through flexible couplings, and the geared 
system is driven by a 149kW DC motor. Vibrational level at various positions on the 
casing plate are measured using PCB 303 A and Endevco 2271 accelerometers over the 
operational speed range of 2250-5750 rpm which corresponds to the gear mesh 
frequency coj, range of 1050-2683 Hz. 

4.7.2 Vibration Models 

A dynamic finite element model of Figure 4.16b is developed using a commercial 
software [39] as outlined in Section 4.4.1. The magnitudes of k^y, which are included 
in the FEM model through a 6 dimensional generalized stiffness matrix element, are 
computed using the two nonlinear algebraic equations (4.27) which neglect the static 
elastic deformation of the casing plate. The shafts are modeled using 2 noded 
Timonshenko beam elements with axial degrees of freedom capability. Four noded 
quadrilateral plate elements with shear deformation and rotary inertia effects are used to 
construct the flexible casing. The shafts and torsional couplings exterior of the test 
gearbox are modeled as torsional stiffness elements only, and the gear, pinion, motor 
and load are incorporated in the FEM model using generalized mass and inertia 
elements. Gear mesh coupling between the gear and pinion is described by a generic 
stiffness matrix of dimension 6, similar to [K]j, m . Free rotational 0^(0 boundary 
condition is specified at the motor and load inertias, while ideally clamped u wa =0 wa =O, 
w=x,y,z, boundary points at the comer of the bottom plate are assumed to simulate the 
rigid mount feet. 

Natural frequencies and modes predicted by a similar FEM model have already 
been verified by comparison with experimental modal analysis in an earlier publication 
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by us [50]. Hence, eigensolutions will not be discussed here. But, it may be noted that 
the conventional geared rotor dynamic models with simple bearing models and without 
casing dynamics are not valid for high excitation frequencies (beyond the first casing 
plate elastic mode). For the forced response study, only the fundamental harmonic C 0 h 
of the static transmission error excitation e(t) = 3.5 sinCw^t) pm computed from the gear 
tooth profile and geometry [46]. About 110 dynamic degrees of freedom are selected to 
minimize computational effort while still maintaining sufficient accuracy. Modal 
damping ratio £ = 0.05 is assumed over the frequency range of interest 


4.7.3 Casing Response 

Several locations on 3 different casing plate surfaces as shown in Figure 4.17 have 
been chosen for the experimental validation of our theory. Figure 4.18 compares 
predicted and measured mean square acceleration spectra at ©h over 0 kHz to 4 kHz 
even though measurements have been conducted only between 1 kHz to 3 kHz. In 
general, good agreement is found between experiment and theory. Some discrepancy is 
observed above 2.4 kHz which is due to the limitation of the FEM model in the high 
modal density regime. It may be noted that the simple theory which utilizes only 
conventional bearing models cannot predict flexural vibration of the casing plate as 
shown in Figure 4.18. A broad band vibratory energy comparison is achieved by 
averaging the mean square value of the acceleration over one-third octave bands, in 
dB per unit frequency bandwidth Ag), which is defined as follows 


L ,dB = 10 log 
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(2Aco)S Re • 

A C0 
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; <A 2 > ref = 1.0 g 2 (4.35) 
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Figure 4.18 
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Comparison of casing flexural vibrations for example case III with 8 zm = 
0.08 mm. (a) Side plate supporting the bearings and facing the 
dynamometer (location 18a). (b) Top plate cover (location 18b). (c) Side 
plate with no bearings (location 18c). 
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where A* implies the complex conjugate of the acceleration. Our predictions are within 
+5.0 dB of the measured values as shown in Table 4.8. 

The sensitivity of results to [K] bm is evaluated next. Recall that [K] bm is 
computed based on the knowledge of the mean axial displacements 8 zm applied to the 
bearings. This parameter has been assumed to be a constant in our analysis although its 
exact value for the NASA experimental setup is not known. Using 2 reasonable values 
of 8 zm , we obtain a range of acceleration spectra in Figure 4.19. Almost all the 
experimental data are now within this range. This explains experimental scatter 
observed in Figure 4.18 as experimental 8 zm may vary slightly from one steady-state 
speed to another during testing. 

In example case II, we found that only a few coefficients in [K] bm are necessary 
depending on the excitation. Accordingly, we had used only k byy , k by Q x and k b Q x Q x . 
But in the present FEM model of case III, we must include the entire symmetric [K] bm 
matrix of dimension 6 with no simplification at all. Now if we omit the off-diagonal 
term k bx Q y which is not directly excited by e(t) in the y-direction of Figure 4.1 while still 

retaining other k b y coefficients, no significant changes in our predictions are seen. 
Similarly, if the off-diagonal terms k bzy and k bz Q x which are related to the shaft axial 

degrees of freedom are also neglected, predicted acceleration spectra vary slightly. 
However, if k by Q x and k b Q x Q x which constraint the bending motion of the shaft excited 

by e(t) are assigned zero values per simple theory, zero out-of-plane acceleration levels 
are observed. This confirms that the vibration transfer through the bearings, from the 
shaft bending motion to the casing flexural motion, is highly dependent on these two 


terms. 


Table 4.8 Comparison of (dB) for example case HI. 


Center 

Frequency (Hz) 


Location in Figure 4.17 

18(a) 

18(b) 

18(c) 

1250 

Theory 

-3.0 

-13 

-4.5 


Experiment 

-2.0 

-17 

-8.4 

1600 

Theory 

-5.5 

-7.0 

-11 


Experiment 

-10 

-8.1 

-14 

2000 

Theory 

-11 

-5.4 

2.1 


Experiment 

-10 

-7.4 

0.0 

2500 

Theory 

-3.0 

0.5 

5.9 


Experiment 

-7.3 

-1.0 

4.0 
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Figure 4. 19 Sensitivity of predicted casing acceleration to preload § zm for example case 
III. (a) Side plate supporting the bearings and facing the dynamometer 
(location 18a). (b) Top plate cover (location 18b). (c) Side plate with no 
bearings (location 18c). 
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4.8 CONCLUDING REMARKS 

A generic geared rotor system model has been developed using lumped parameter 
and dynamic finite element techniques which incorporates a new mathematical model for 
precision rolling element bearings proposed earlier in Chapter II. This system oriented 
model includes internal rotating system, rolling element bearings, flexible torsional 
couplings, motor and load, flexible or rigid casing, and compliant or massive mounts. 
The discrete shaft model excludes the effect of gyroscopic moment, but includes 
torsional, flexural and longitudinal motions. In example cases I and II, only flexural 
and torsional motions are predicted. The effects of mount stiffness, casing mass and 
bearing preload on the overall dynamic behavior have been investigated through 3 
example cases of single-stage rotor and geared rotor systems excited by rotating mass 
unbalance at shaft frequency and kinematic transmission error excitation at mesh 
frequency respectively. The results indicate that our proposed model is clearly superior 
to the conventional simple models as given in the literature. For instance, we are able to 
predict bearing and mount moment transmissibilities and improve casing flexural 
vibration prediction significantly using our theory; measurements made on case III 
validate our formulations partially. Conversely, the simple models are unable to account 
for any bearing moment transmissibility, and consequently predict zero flexural 
response on the casing. 

Natural frequencies of the example case III geared rotor system increase with 
higher mount stiffness and bearing preload, and lower casing mass. Transition mass 
and stiffness points which separate regimes of a natural frequency diagram are 
predicted. Through extensive parametric studies, we find that bearing transmissibilities 
may be reduced by using a higher bearing preload. However, this may result in higher 
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dynamic transmission error for the spur gear pair, mostly due to the coupled flexural- 
torsional motion. Mount transmissibility is affected by changes in casing mass, mount 
stiffness and mean bearing axial displacement, but bearing transmissibility is most 
sensitive to the mean bearing axial displacement. Similar conclusions can be drawn 
regarding the other two example cases. 

Our theory, though restricted to the linear and time-invariant dynamic system, is 
comprehensive. It can be used for analysis as well as design studies of other rotating 
mechanical systems with multiple shafts and gear pairs or multi-staged rotors. We are 
currently extending this formulation to examine the vibratory energy transfer over a 
broad band basis. 


CHAPTER V 


STATISTICAL ENERGY ANALYSIS 
5.1 INTRODUCTION 

It has been shown in Chapters H-IV, using classical lumped parameter and 
dynamic finite element techniques, that the proposed bearing model is clearly superior to 
the existing simple models for predicting vibration transmission through bearings in a 
geared rotor system. Although the proposed models have been shown to be reliable up 
to a moderately high frequency, it is conceivable that these models are inadequate at very 
high frequencies where the modal density is high. Classical vibration models do not 
predict modes accurately in this frequency regime, and even if it is possible to do so by 
employing closely spaced nodal points, such models require a significantly large 
computational effort Moreover, the vast amount of predicted response spectra at many 
spatial points would be difficult to interpret. Accordingly, asymptotic or statistical 
methods must be adopted; typical techniques include the statistical energy analysis [56- 
60], asymptotic modal analysis [61-63] and asymptotic analysis using infinite system 
impedances [64]. 

This study concentrates on the development of a broad band vibratory energy 
transfer model for a geared rotor system with the proposed bearing model using the 
statistical energy analysis (SEA) method. This method has been applied successfully to 
a wide variety of structural dynamic and acoustic systems with large number of modes 


155 


156 


[56-60,65-69]. However several unresolved research issues still exist [57,61,62,70, 
71]. The specific objectives of this study are to: (i) conduct modal analysis of a geared 
rotor system with flexible casing and mounts, (ii) investigate the feasibility of applying 
SEA to this problem, (iii) analyze the following 4 example cases using SEA: I. a plate- 
cantilevered rectangular beam, II . case I with circular shaft-bearing system replacing the 
cantilevered beam, IQ. a circular shaft-bearing-plate-mount system, and IV. a simple 
geared rotor system, and (iv) perform parametric studies to examine the characteristics 
of vibratory energy transfer through bearings, and mean square vibroacoustic response 
of the casing. The first and second example cases are revised and extended versions of 
a study performed by Lyon and Eichler [56,58]. The later two example cases are taken 
from Chapters HI and IV where these systems have been studied at low frequencies. 
Experimental validations are also included. 

5.2 MODAL ANALYSIS OF GEAR CASING AND MOUNTS 

Analytical and experimental modal analyses of a real gear casing and mounts have 
been conducted to investigate the feasibility of using SEA. Natural frequencies coj and 
modes <j>j are calculated using a commercial finite element method (FEM) program [39], 
and predictions for a rigidly mounted, stiffened gearbox are verified by the experimental 
modal analysis (EMA). 

5.2.1 Finite Element Model 

The rectangular gearbox as shown in Figure 5.1a is approximately 0.254m x 
0.279m x 0.330m (lO.Oin x ll.Oin x 13.0in), and all of its plates are 0.006m (1/4 in) 
thick made of 1020 steel except the regions near the bearings which are 0.025m (l.Oin) 
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(b) 


Figure 5. 1 Gearbox with rigid and flexible mounts. 
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thick. There are four circular holes for the bearings, two at each side plate supporting 
the shafts. Figure 5.1b illustrates the 0.254m (lO.Oin) tall flexible mount frame which 
is constructed from eight 0.006m (1/4 in) thick, 1020 steel angle beams with three 
different lengths. Fuselage sheet, as shown in Figure 5.1b, is attached horizontally to 
the flexible mount structure. The casing is supported at each comer of the base plate for 
all mounting conditions, and the mounts are attached to a rigid foundation. 

Two FEM models of the rectangular gearbox without its spur gears set, shafts and 
bearings are developed for the rigidly and flexibly mounted casings. These FEM 
models consist of four-noded quadrilateral plate elements with bending and membrane 
capabilities for the housing and attached fuselage, and two-noded shaft element with 
shear deformation and rotary inertia capabilities for the flexible mount skeleton and 
housing plate stiffeners. The boundary conditions are: (i) zero displacements and 
rotations at each comer of the base plate for the rigid mount, and (ii) similar conditions 
at each foot of the flexible mount. The interfaces between adjacent housing plates are 
assumed to be continuous. About 100 dynamic degrees of freedom are specified to 
reduce computational effort while still maintaining sufficient accuracy. Natural 
frequencies are computed up'to at least 4 kHz to cover the gear mesh frequency regime. 

5.2.2 Experiments and Model Validation 

Modal experiments have been performed on a NASA high precision gearbox with 
the spur gear set, shafts, and ball bearings installed. An approximate configuration of 
the NASA gearbox is shown in Figure 5.2. The nominal dimensions of the gear 
housing have been given in the previous section. The variable center distance gear-shaft 
pair is supported by four ball bearings. Four side plates and a base plate are welded 
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Figure 5.2 NASA gear test facility [55]. 
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together while the top plate is bolted to the side plates. The housing plates are stiffened 
internally, and the gear housing system is mounted rigidly to a massive foundation. 
Dynamic transfer functions are obtained only on the exterior of the gear casing structure 
using the GENRAD 2515 system [72]. For these experiments, 154 measurement points 
have been selected in the direction transverse to the plane of the casing plates with the 
reference point being approximately near the center of the top plate to avoid nodal points 
of interest. Natural frequencies and modes are then estimated using the MODAL PLUS 
program [73]. Here, the exponential method has been used to extract modal parameters 
and generate analytical functions for the transfer functions, while the circle fit method 
has been used to construct the modal vectors. 

Figure 5.3 compares predicted and measured C 0 j, and FEM is found to be in good 
agreement with EMA. For each mode <j>j, two simplified illustrations are shown in 
Figure 5.3: (i) mode shape of the three visible plates ,and (ii) mode shape of the three 
nonvisible plates in an approximate isometric view. The higher modes, not shown here, 
are also given by similar combinations of plate flexural motions. Comparison between 
theory and experiment for the higher modes is made on the basis of number of modes 
within each one-third octave bands over 500-4000 Hz in Figure 5.4a, because of the 
high number of participating modes observed. The results again indicate that FEM is in 
good agreement with EMA. 

5.2.3 Parametric Studies 

The rigidly mounted gear housing is observed to possess only the elastic modes of 
the casing plates. On the other hand, FEM model of the flexibly mounted gear casing 
indicates that the first six modes are translational and rotational rigid body modes of the 
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(o 1 = 501 Hz 
(497 Hz) 


co 2 = 598 Hz 
(584 Hz) 



w 3 = 627 Hz © 4 = 752 Hz 

(639 Hz) (781 Hz) 



®5 = 889 Hz 
(877 Hz) 


Figure 5.3 Mode shapes of the rigidly mounted, stiffened NASA gearbox as predicted 
by FEM (EMA results are given in parenthesis). Bold sign implies larger 
amplitude. 
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the casing as shown in Figure 5.5; for example, <»i = 54 Hz corresponds to the casing 
vibrating in the Y-direction as shown. In addition, coj are considerably lower, by 
approximately one order of magnitude as compared to those of the rigidly mounted gear 
casing. These rigid body vibration modes result from the complex elastic deformations 
of the flexible mount skeleton and fuselage sheet. The casing plate natural frequencies 
are also lowered, especially the first few, when the box is mounted flexibly. Figure 
5.4b compares the number of modes of the flexibly mounted gear casing to that of the 
rigidly mounted one. High modal counts above 1 kHz are seen. 

The introduction of gear casing plate stiffeners as shown in Figure 5.2 does not 
change the nature of the mode predictions; Also, the natural frequencies for this case 
only vary slightly; lower modes are affected more by the stiffeners than the higher ones. 
Note that the numbers of modes in each 1/3 octave band over 400-4000 Hz range remain 
nearly the same as evident from Table 5.1. 


Table 5.1 Number of modes for the stiffened and unstiffened gearboxes mounted 
rigidly. 


1/3 Octave Band Center 












Frequency (kHz) 

0.4 

0.5 

0.63 

0.8 

1.0 

1.25 

1.6 

2.0 

2.5 

3.15 

4.0 

Stiffened 

0 

1 

2 

3 

1 

5 

7 

6 

10 

11 

11 

Unstiffened 

1 

3 

2 

1 

2 

6 

8 

7 

11 

10 

11 
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5.3 JUSTIFICATION FOR USING SEA 
5.3.1 Modal Densities 

Analytical and experimental modal analyses of the NASA gear casing clearly show 
that modal density tends to be high at higher frequencies. Although the analysis is 
restricted to a gear casing system, it is reasonable to assume that similar results are valid 
for other rotating mechanical system housings. Accordingly, SEA can be used 
justifiably. 

Next, consider the modal densities of a shaft or beam n s and rectangular plate n c , 
given for bending motion with simply- supported boundary conditions [56] 



A c / 3p c (1 - p. 2 ) 
n c" h c V E 


(5.1) 


where p is the material density, A c is the plate surface area, L s is the shaft length, E is 
the modulus of elasticity, I s is the area moment of inertia of the shaft, co is the 
bandwidth center frequency, h is the plate or shaft thickness, jj. is the Poisson’s ratio, 
and the subscripts s and c denote shaft and plate respectively. For typical numerical 
values given later in Section 5.4 we find that: (i) n c = 0.227, and (ii) n s = 0.012 at 100 
Hz, n s = 0.004 at 1000 Hz and n s = 0.002 at 4000 Hz. Note that n c » n s as expected. 
Although the shaft has a fairly low number of modes in the frequency range of interest, 
SEA is still valid since the plate modal density n c is very high. 
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5.3.2 Literature Review 

A majority of publications are on the application of the SEA procedure to 
dynamic systems with high n such as structural-acoustical interactions in a fuselage 
[56,65,74], sound transmission through panels [60], and vibratory energy transmission 
in mechanical equipment [56]. Of interest here are analytical or experimental estimation 
of SEA parameters for simple structural systems described by a flat plate, a cylinder 
and/or a shaft [56,58,59,66]. In these studies, structural connections are often assumed 
rigid such as in the ideally welded case. 

Lyon and Eichler [56,58] in 1964 and Lyon and Scharton [59] in 1965 developed 
analytical expressions for the coupling loss factor t| in several connected structures, 
such as a plate bonded to a cantilevered beam [56,58]. Here r| was derived assuming a 
semi-infinite beam attached to an infinite plate and by further assuming that only a 
dynamic moment coupling at the joint can describe the motion/force transmission 
phenomenon. This problem is re-examined in this study and is then extended to a 
circular shaft-bearing-plate-mount system. In addition, Lyon and Eichler [56,58] also 
developed SEA models of two structures inter-connected through a single (scalar) 
stiffness element. Typical examples include two longitudinal rods connected via a linear 
stiffness element, two discrete masses coupled by a linear spring, and a plate attached to 
a single degree of freedom resonator [56]. A similar analysis on the longitudinal 
vibration of linearly coupled rods was performed by Keane and Price [69] in 1987. 
Loss factors of typical line or point connected structures such as a plate welded to a 
cylinder, cross beams and two perpendicular plates bolted or welded together, have been 
calculated assuming ideal rigid joints [56,58,59,75,76]. But a compliant bearing system 
problem is yet to be analyzed. 
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5.4 EXAMPLE CASE I: COUPLING LOSS FACTOR OF PLATE- 
CANTILEVERED BEAM SYSTEM 

First, we attempt to rework the plate-cantilevered beam problem of Lyon and 
Eichler [56,58] as shown in Figure 5.6. Only flexural motions of the plate and beam are 
considered in this case. Accordingly, Lyon and Eichler [56,58] developed an 
expression for the coupling loss factor which describes the vibratory energy transfer 
between the beam(s) and the plate(c) due to a moment coupling at the joint 


Tlsc(®) 


(2p s A s K s c s) 

com „ 


Z c 

Z c+ Z s 



(5.2) 


where k=^/ l s /A s is the radius of gyration, c=Ve/P s is the wave speed, m is the mass 
and Re( ) implies the real part of a complex variable. The driving point moment 
impedances for the plate 2^ and beam are [58,77] 


Z c (co) = 


16P c h c Kg c? 


CO 


O-lfmCW) ’ 


Z s (co) = 


Ps ^s^s c s^s 

CO 


( 1 — i ) 


(5.3) 


where k is the wavenumber and h is the plate or beam thickness. Here, note that Z s of 
equation (5.3) differs from Lyon and Eichler's expression [56,58] in the sign of the 
imaginary part which is probably a typographical error. They assumed that IZ C I»IZ S I, 
k s =k c , c s =c c and p s =p c , and simplified equation (5.2) to yield a frequency invariant 
expression for as 


n S c = W s /(4L s ) 


(5.4) 
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Figure 5.6 Schematic of example case I: a plate-cantilevered beam, (a) Experimental 
finite System, (b) Theoretical infinite system. 
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where W s and L s are the beam width and length respectively as shown in Figure 5.6. 
However, our calculations indicate that IZ C I>IZ S I but not IZ C I»IZ S I for the parameters 
used by Lyon and Eichler [58], as shown in Figure 5.7a. Using equation (5.2), r] sc is 
recomputed and compared with the approximate model given by equation (5.3) in Figure 
5.7b. Experimental results given by Lyon and Eichler are also plotted in Figure 5.7b. It 
can be seen that our calculation is better than earlier prediction. At low frequency, a 
discrepancy is observed between theory and experiment which is primarily due to the 
low modal density in this regime. The presence of a low natural frequency may be due 
to the compliant epoxy bond between the beam and plate. However, above the 
threshold frequency where many modes participate, shown as a vertical line in Figure 
5.7b, the slope of the least square straight line fit on the experimental data is nearly the 
same as the predicted , n sc . 

5.5 EXAMPLE CASE II: COUPLING LOSS FACTOR OF CIRCULAR 
SHAFT-BEARING-PLATE SYSTEM 

Next, we modify Figure 5.6 by inserting a ball bearing between the circular shaft 
(which replaces the beam in Figure 5.6) and the rectangular plate. Again, a semi-infinite 
shaft and an infinite plate are assumed. For SEA, we reduce the system to a plate 
subsystem and a shaft-bearing subsystem. The coupling loss factor T] sc is still given by 
equation (5.2), but Zg must be modified to account for the compliant bearing. 

Consider a shaft with boundary conditions shown in Figure 5.8. The bearing end 
is subjected to zero transverse velocity v ysa (t,z=0')=0 and a sinusoidally varying 
moment M xsa (t,z=0')=M xsa e 1(ot . Here, z=0' refers to the junction point between bearing 
and plate, and z=0 + is the junction point between the bearing and shaft. In the 
frequency domain, using the definition of bearing force vector 



Coupling Loss Factor, 
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(b) 


Figure 5.7 Comparison between Lyon and Eichler's approximation [56,58] and our 
proposed formulation for example case I. (a) Comparison of IZ C I/IZ C +Z S I. 
(b) Comparison of predicted r| sc with experimental results given by Lyon 
and Eichler [56,58]. 
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M xsa( £ °. z =0") = M xsa 

Vysa(o >,z=(T)=0 



z=0 


Figure 5.8 Boundary conditions for example case II: semi-infinite shaft-bearing 
system. 
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{ f(co) } b=[K] b m { v(C 0 ) } b /ito in terms of the velocity vector {v(co)} b and the proposed 
bearing stiffness matrix [K] bm , v ysa (co,z=0')=0 is shown to be equivalent to 
Fy Sa (o),O + )=(k{jy9 x /k| ) 0 x 0 x )M XS a+ v y S a(®>O + ) { (^by0x/^b0 x 0 x )^b9 x 0x — ^byy }/(ico). 

The governing equations for M xsa (co,z) and F ysa (co,z) of the shaft in terms of Vy^CD.z) 
are 


M EI s dv ysa 

*** iC0 dz 2 


F = - 


EI C d v 


ysa 


ysa ico ^3 


(5.5) 


In general, the bounded solution for v ysa (CO,z) is 


v ysa (co,z) = { b L e i k * z + b2 e k * z } 


(5.6) 


Using equations (5.5) and (5.6), we enforce the boundary conditions for the shaft at 
z=0 + to obtain the following closed form solution for the coefficient vector 
{b)={b|,b 2 ) T . 
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Hence, the bearing-shaft impedance Zg is then given by 


Z s (co) = 


M xsa ( £0 » 0) 


‘b0 x e x 


dv 


ysa 

~dz 


dv 


(5.8) 


(co, 0") 


ICO + 


k. a v (to, 0 ) -i-k^. _ , 

by0 x ysa ' b0 x 0 x dz 


ySa (CO, 0 + ) 


M 


xsa 


Although M xsa appears in equation (5.8), Z s is independent of the magnitude of M xsa 

dV ygg 

since v ysa (co,0 + ) and —— — (to, 0 ) are linearly proportional to in equation (5.5). 

Using equations (5.2), (5.7) and (5.8), T[ sc is computed for several bearings 
whose stiffness coefficients are tabulated in Table 5.2. Figure 5.9 compares these 
where the material and geometrical properties of the shaft and plate are equivalent to 
those used by Lyon and Eichler [56,58] in example case I. For a very soft bearing, Z s 
is dominated by the bearing parameters, and hence l2^ c l»IZ s l is valid. This implies that 
rjgc reduces to the frequency invariant expression given by equation (5.4) as shown by 
set A in Figure 5.9. On the other hand, rj sc for a rigid bearing (say set B or C) is 
typically smaller than rjsc for a soft bearing given by A. The extreme values of t^, as 
k b jj or bearing preload becomes very large, depend on the relative magnitudes of the 
bearing stiffness coefficients as shown in Figure 5.9 for sets D and E. Such deviation 
for a stiff bearing from set B is primarily due to the presence of off-diagonal stiffness 
coefficient such as k by Q x . However it is found that T) sc for most stiff bearings will 
approach set E in Figure 5.9 which is obviously identical to the prediction for the 
example case I using equations (5.2) and (5.3). 
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Table 5.2 Typical bearing stiffness coefficients of Figure 5.9 for example case II. 


Set 

*- 
c r 
v- 
v; 

! 

^by0 x (N) 

k b 0 x e x (Nm) 

A 

Very Compliant Bearings (k^y < 1E3) 

B 

1E5 

5E4 

2E4 

C 

5E5 

8E4 

5E4 

D 

1E8 

1E7 

1E6 

E 

1E8 

3E7 

1E7 


5.6 EXAMPLE CASE III: A CIRCULAR SHAFT-BEARING-PLATE 

SYSTEM 

5.6.1 Theory 

Consider a circular shaft-bearing-plate-mount system similar to example case II but 
with a shaft of finite length as shown in Figure 5.10a. Note that this system has been 
analyzed earlier using the deterministic vibration modeling technique for low frequency 
response. Recall the unconstrained end of the non-rotating circular shaft is subjected to 
a harmonically varying force F ys (t) = F ysa e icot + F ysm , where F ysm is the mean force. 
From Chapter III, we know that the longitudinal and torsional motions of the shaft, and 
the in-plane vibration of the plate can be neglected. Hence, the shaft bending vibration 
and plate flexural motion are of interest here. Also, the previous deterministic vibration 
models indicate that the coupling between these two motions is mainly due to the 
dynamic moment at the bearing provided the longitudinal shaft vibration is not excited. 
Accordingly, two subsystems which can be easily identified using SEA are the 
transverse modes of the shaft-bearing system and the flexural modes of the plate-mount 
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Flexible 
Shaft (s) 




Applied Force 



Input Power 
n s I 




Shaft-bearing system 
(bending modes) 

E s (co), n s (co) 

Tls 

► Energy Dissipation 

Hsc 

Energy Transfer 

f 

Plate-mount system 
(flexural modes) 

E c (co),n c (co) 

► Energy Dissipation 


(b) 

Figure 5.10 Example case III: (a) schematic of the circular shaft-bearing-plate-mount 
system and (b) A SEA model of the shaft-bearing-plate-mount system. 
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system as shown in Figure 5.10b. We now proceed to derive T| sc governing the 
vibratory energy flow between these two subsystems using the driving point junction 
impedance method suggested by Lyon [56]. 

Since the rectangular plate dimension is much larger than the bearing dimension, 
the plate is again assumed infinite. Consider the vibrational power flow n sc from the 
shaft-bearing subsystem (subscript s) to the plate (subscript c), due to the uniform 
external Gaussian random force over a frequency bandwidth Aco with center frequency 
co, 


11^(0)) =cqti sc (co) n s (co) 


E s (o» 

n s (co) 


E c (to) ' 
n c (<0) , 


(5.9) 


where nj and Ej (j = s, c) refer to subsystem modal density and total vibratory energy 
respectively. Since the plate is assumed to be reasonably well damped and geometrically 
large, equation (5.9) is approximated assuming n c »ns or to yield 


n K (fi)) 

= CD E s (co) 


(5.10) 


For the shaft, E s = m s <V s > where m s is the shaft mass and <V S > is mean square shaft 

transverse velocity averaged over Aco and shaft length L s . Using the expression for 
given by Lyon and Eichler [58] as discussed in the previous section, in terms of Z s 
and Zq is identical to equation (5.2) which is rewritten here for the circular shaft case as 


Tlsc^) 


4 El s 
coL s 


Re 


f ll 

z c 

o | 

Z c + Z s 


2 


(5.11) 


178 


For an infinite steel plate of thickness h c , the point moment impedance Z c is given by 
equation (5.3); it is rewritten as 

4 Eh ^ — ] 

Z C (G)) = { * ~ “tT^ 1 ^20’*M s )}‘ * Ms « 1 (5.12) 

The shaft-bearing impedance Z s in equation (5.1 1) is derived next by solving the 
boundary -value problem for the system shown in Figure 5.11a. The free end of the 
shaft has a vanishing dynamic shear force F ysa (o>,z=L s )=0 and bending moment 
M xsa (co,z=L s )=0. The bearing end is similar to Figure 5.8. Following the same 

argument used previously for example case II will lead to the two boundary conditions 
described by M xsa (co,z=0 + )=M xsa and F xsa (co,z=0 + ) = (k bx0y /k b0y0y )M ysa + 

v y Sa (co,0 + ) { (k bx0y /k b0y0y ) 2 k b0y0 y -k bxx } I (ico). Governing equations are still given 

by equation (5.5), but the solution for v ysa (o),z) is assumed to be of the following form 
[77] 

v ysa( t > z ) = { bi e -JM + 1>2 e i k s z + b 3 e -k s z + b 4 e k s z } (5.13) 

Using equations (5.5) and (5.13), the boundary conditions yield the following algebraic 
problem 

[B] { b } = M ysa { k b y 0x /k b 0 x 0 x , ico/(EIk s ), 0, o} 


(5.14) 
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^xsa^* z ® ) ^xsa 

v ysa(CO,z=0 _ )=° 
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-► z 


z=0 


M xsa(®’ z— ^«) ~ 0 


Flexible Shaft 


(a) 


xsa 

F ysa (co, z = L s ) =0 


v V sa(° ) , z =- L s)=0 


d v 


ysa 

ysa 


dz 


(co,z=-L s ) =0 

0 


0 

Bearing 


z=-L s 


M X sa(° ) » z = 0 ) =0 


F ysa (ci),z = 0) = F 


xsa 


Flexible Shaft 


z=0 


(b) 


Figure 5. 1 1 Boundary conditions for example case HI: finite shaft-bearing system. 

(a) Moment applied at the bearing end. (b) Force applied at the free end. 
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The nonzero elements of coefficient matrix [B] of dimension 4 are 


B n = t 


B 13 CO 


..2 

lk. a 

Elk l - ik + -r- — 
v byy k be x e xy 


Elk o - k L + 


by 0, 


s byy k 


b0 x 0xy 


1 B 12 = ffl 


; B 14 = W 


..2 

lk, Q 

3 by0 x 

- Elk . - ik . + 


s “ byy k 




b0 x 0 x j 


3 by0 3 

- Elk * — k, +■ 


s byy k 


b 0 x 0 x j 


B 21 -B 22 _ ~ B 23~ -B 24 -1 » B 31 _-iB 41 -ie '* ksLs * B 32 -iB 42 - i^ ksLs *» 

B 33 = " B 43 =_e_ksLs ; B 34 =B 44 =e k s L s (5.15) 


Both [B] and {b} = { b x , b 2 , b 3 , b 4 } T can be easily obtained numerically. The bearing- 

shaft impedance Zg is still given by equation (5.8). 

The same procedure may be applied to obtain the driving point force impedance for 

a harmonically varying transverse force Fy Sa (t,z=0). Note that the origin is redefined at 

the forcing point as shown in Figure 5.11b for convenience. The boundary conditions 

dv vsa 

are Fy Sa (co,z = 0)=Fy Sa , M xsa (co,z— 0)=0, Vy Sa (co,z=- L^)=0, and ^ (co, L s ) = 0. 
These conditions at z=-L^ can be rewritten for z=-L* like the previous case to obtain 

^xsa(®> z=— Lg) 55 - (kjj 0 x Q x — (co, — L s ) +kby 0 x v ysa (co,-L s )) j (i co) , 

dv ^ " 1 

Fy S a(co,z=— Ls) =— (kby 0 x — (co, — L s ) + ^byy^ysa^' - ■ (5.16) 

These prescribed boundary conditions again yield a set of algebraic problem similar to 
Equation (5.14). The nonzero elements of the coefficient matrix [B] of dimension 4 are 
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Bn - -B 2 i - B 12 - B 2 2 - ~Bi 3 ~ ^23 - _ Bi4 “ *®24 “ *0 5 
B 3i=<- j fe s k h« e +k b,e ; 


L s b0„0. 


by0, 


B 32 = (j k S Xe I + k b^, + EIk s)'' ilt,L ' 

B 33 = ^ -k s k b0 6 + k by0 “ Elk s ) e ; 


B„ = (k s k k0 „ +k k „ - Elk j )e -klL| ; 


34 '~ s ~b$ x B x ' “byO, 


B 4i = < ik s k b,0 ;i - k lw - iEIk ^ it,L ‘ ; 

B 42 = ( — ik s k bv0 - k bw + iHk s > e ~ ' l ' L ' 


s by0 byy 


3 x 


B 43 = ( k 4 k by0 l - k IW + EIk S ) e 
B 44 = <- k s k bye - k by y - Hk s>«' , ‘ ,L ' 


(5.17) 


The right hand side vector (b) of the algebraic problem is {0,F ysa co/ (Elk^),0,0 } T . 
Force impedance at the driving point is then given by Z s (©,z=0)=F ysa /v ysa (co,0). 
Accordingly, the input power is n s = (1/2) F ysa Re{(l/Z s )*} where Re{ } is the real 
part of the complex variable and ( )* implies the complex conjugation. 

We can now compute the vibratory energy transfer n sc through the bearing and 
steady-state subsystem energy levels E s and E c by applying the energy balances to both 
subsystems shown in Figure 5.10b; here T| cs = rj sc 


T ls +T l S c(®) 

“Ties (CO) 

[e s (0>)1 _j 

n s /a> 

-Tlac(fi>) 

Tlc+llcsC®) . 

1E C (0))J 1 

l o J 


(5.18) 
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E,(£D) 

E c (co) 


Hii^c + y 

©(^c + ^sc + ^cs) 
E s *1 sc 

©OlsHc + ^sc + ^cs) 


(5.19) 


Since E = m <V^>, the following velocity levels may be obtained at any center 
frequency to from either equation (5.18) or (5.19) 


<V C > ^ m s( *1 sc (°>) > 

<V S 2 > m chc +1 lcs(®)J 


(5.20) 


5.6.2 Validation and Parametric Studies 

In order to validate our SEA formulation, we compare the mean square mobility 

level of the plate with experimental data reported earlier in Chapter III. Note that 

although all nonzero bearing stiffness coefficients k b y are computed and given in Table 

3.6, only kjj yy =3.69E8, k by g x =3.52E5 and k b Q x e x =4.19E4 are used as they appear to 

be the most significant ones according to the proposed theory. Using equation (5.19), 
2 

<V > is computed and compared with experimental results m Figure 5.12. Theoretical 

predictions for three values of dissipation loss factor T| s =T| c ^r| s (co) are given since the 
choice of structural damping is critical to the SEA analysis. It can be seen from Figure 
5.13 that the experimental data are approximately bounded by r| s =0.0003 and r| s =0.03. 
Here T| s is also assumed to be frequency invariant which may not be valid in the 
experiment. Accordingly, comparison between theory and experiment is deemed to be 
excellent. 

Further comparison between theory and experiment can be made for the case of a 
semi-infinite shaft considered in the example case II. Using equations (5.2), (5.7) and 
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1/3 Octave Band Center Frequency (Hz) 


Figure 5.12 Comparison between theory and experiment [17] for example case in with 
a finite shaft. Here plate mobility level in dB is 101og[o<V^> re 
<vj>= 1 .0(m/s) 2 for F ysa = 1 .ON. 
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(5.8) for tj sc , the mobility levels are computed and are found to be given by straight 
lines as shown in Figure 5.13. These lines represent the asymptotic behavior of the 
system when the shaft is very long i.e. L s — > <». Also, note that Figure 5.13 is 
consistent with the trends reported by Lyon and Eichler [58] for the plate-cantilevered 
beam problem. Again, most experimental data are bounded within the range given by 
T] s = 0.0003 to 0.03. 

Consider the finite shaft length L s =1.32 m of high modal density n s . In Figure 
5.14 t] sc is compared for a finite shaft and semi-infinite shaft. It can be seen that the 
result for the semi-infinite shaft follows the average values of the finite shaft. 

Next, the effect of bearing preload or stiffness coefficients on the mean-square 
velocity response of the plate-mount system is evaluated. Equation (5.20) is used to 
predict the curves shown in Figure 5.15 while keeping other system parameters the 
same. The bearing coefficients used in this analysis are from Tables 3.3 and 3.4. Three 
bearing preloads used here are 115N, 190N and 285N which are referred to as low, 
medium and high preloads. We observe minor changes with preload except in the 
vicinity of 4 kHz. 



Coupling Loss Factor, 





<v?> 


8 ’ 



1/3 Octave Band Center Frequency (Hz) 


Figure 5.15 Effect of bearing preload on the normalized casing plate response (example 
case III). 



188 


5.7 EXAMPLE CASE IV: A GEARED ROTOR SYSTEM 
5.7.1 Assumptions 

As the final case, we investigate a generic geared rotor system with 4 bearings as 
shown in Figure 5.16a. This system is not only a natural extension to example HI, it 
can also be treated as an approximate model of the NASA gear test facility in Figure 5.2; 
for gearbox details, refer to Chapter IV. The intent of the SEA analysis is to predict 
asymptotic casing plate response and radiated sound pressure from the gearbox, on an 
order of magnitude basis. Accordingly, the following assumptions are made: 

1 . Consider the case of 4 identical rolling element bearings and 2 identical shafts 
carrying gears. 

2 . The vibratory source associated with the static transmission error excitation is at the 
gear mesh. The net exchange of energy through the gear mesh coupling is 
neglected and the vibratory energy is assumed to flow outward from the gear mesh 
source to the casing through shafts and bearings. The source is assumed to be far 
away from the bearings. 

3. Only shaft bending vibration is coupled to the casing plate flexural motion. 
Torsional modes are not included since they are relatively low in number compared 
to the shaft bending modes. In steady-state, both shafts have equal amount of 
vibratory energy (E s ) associated with their bending motions which are uncoupled 
from the external motor and load. 

4. The portion of torsional vibratory energy which is transmitted to the load is 
incorporated in the dissipation loss factor (T| s ) of the shaft-bearing subsystem. 
Here r| s is assumed to be given by Tj s =yn c ^rj s (co) where T| c ^rj c (co) is the loss 
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Figure 5.16 Example case IV: (a) a generic geared rotor system with casing and mounts 
and (b) boundary conditions for a finite shaft-bearing system. 
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factor of the casing; y is expected to be higher than one due to additional losses at 
the bearings and load. 

5 . Bearing dimensions are assumed to be much smaller than the wavelength on the 
casing plate. Also, bearings are not expected to be exchanging vibratory energy. 

6. Modal densities of shaft n s and casing plate are given by equation (5.1). 

7. A diffuse vibration field is assumed for the external casing-mount system. The 
energy dissipation mechanism for this system includes acoustic energy radiated 
from the plate in addition to the structural damping. 

Application of the SEA principle will result in 2 subsystems, one internal 
consisting of 2 shafts and 4 bearings and one external consisting of 2 casing plates and 
mounts similar to Figure 5.10b. The internal subsystem in this case can store the total 
vibratory energy of 2 shafts. Vibratory energy transfer as viewed through Figure 5.10b 
represents the algebraic sum of the energy transfers through 4 bearings. 


5.7.2 Coupling Loss Factor 

The coupling loss factor Tj sc between the internal (shafts-bearings) and external 
(casing) systems is derived, based on the formulation given by equations (5.11) and 
(5.12). This implies that vibratory energy transfer is associated with only the dynamic 
moment at the bearing, and the casing plate is infinite with respect to the bearing 
dimensions. The driving point impedance of the shaft Z s is derived next given the 
boundary conditions shown in Figure 5.16b. At the left bearing end (z=0‘), 
v y Sa (oj,z=0')=0, and M xsa (co,z=0')=M xsa which are equivalent to M xsa (co,z=0 + )=M xsa , 
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and Fy Sa (co,O + )=(k by 0 x /k b 0 x 0 x )M xsa +v ysa (co,O + ) { (kb y e x /k b e x 0 x ) 2 k bex 0 x -k byy } j (ico). 
The other bearing end at z=L* with zero translational and angular velocities is described 
at z=L‘ as 

dv 

M xsa(«.z=L s )=-(k b 0 x 0 x (co, L~ ) +k by 0 x v ysa (co,L*)) / (ico) ; 

dv 

Fysa( co > z= L s ) =_ (k by 0 x - ^ (co, L“ ) + k b yyV ysa (co,L s ))/(ico) (5.21) 


Using equations (5.5) and (5.13), the above boundary conditions are evaluated to 
formulate an algebraic problem of the type given by equation (5.14). The nonzero 
elements of the coefficient matrix [B] are 


B ll co 


B 13 = OO 


..2 

l k, 

EIk s 3 — ik +■■■ y * 

byy Ve x , 


; B 12 = ¥ 


.,2 

CT1 3 -i, b y0, 

- Elk „ - ik , + 


s “ byy k 




3 by0 x 

Elk c - k , + -. x 

s byy • k 


V 


be x 0xy 


; B 14 = i 


-mcl-t. + by9 ‘ 


s byy k 




B 21“ B 22 -_B 23~ _B 24 =1 » 

B 31 = *“-i k s k b« J 0 il + k bye il + EIk s ) e ; 

B 32 = ^ ksk b0 X 0 X + k by9 ^ + EBc s ^ e * b,L * ; 

B 33 = * _ k s k be j 0 ll + k bye }l _ EDt s ) s _k ‘ L ' I 

B 34 = < k S k b0 il 0, +k by0 I - EIk s)' k ' L ' ; 

B 41 =< ik s k b ,0 - k byy- iEIk s)'' ik ‘ L ‘ ! 
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B 42 = ( - * s k bye x - k byy + 

B 4 3 = < k s i 1>y e x - k b yy + Hk s>'' k,L ' i 

B 44 = ^“ k s k by0. - k byy “ ^ s > ek * L> 


(5.22) 


The right hand side vector { b } of equation (5.14) is M ysa { k by 0 x /k b Q x 0 x , 
2 l T 

ico/(EIk^),0,0 j , and the shaft impedance Z s is evaluated numerically using equation 
(5.8). 


5.7.3 Vibroacoustic Response 

Consider the total energy dissipation by the casing as the sum of the energy terms 
associated with structural damping dissipation and energy escape via sound radiation. 
Hence the total dissipation loss factor t| cT for the casing of area A c and mass rr^ is 


r\ cT ((o) = J\ c + 


tom c 


(5.23) 


where z 0 is the characteristic impedance of the surrounding medium, and a c is the 
radiation efficiency of the casing. In this analysis, two radiation efficiency models are 
used: (i) <r cl = 1.0 for an ideal radiator as for many gearboxes the measured radiation 
efficiency has been found to be close to unity [78], and (ii) a c2 for a simply-supported 
rectangular plate [79]. Several investigators have used o c2 successfully in SEA 
applications [60,66,67]. 
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Since this SEA model is very similar to the one developed for example case III, 

solutions given by equations (5.19) and (5.20) are still valid here provided the existence 

of multiple paths, as opposed to only one path in example case III, is recognized. The 
mean- square velocity levels of the shafts-bearings system <V S > and casing-mount 

system <V^> are 


< V, > (“) = 
<V^>(U) = 


n «(V + 4 , 'c«) _ 

m s w ( 'H s 1 ^ c t + ^*1 cT 1 ^ sc + cs) 
411 s TJ ^ 

m c CT + 4T1 cT^sc + 4t 1 s 7 ! cs) 


(5.24) 


Power injected into the system is developed by examining the internal static transmission 
error excitation e(t)=eh(co)e icot . For the spur gear pair, the gear mesh frequency 
co(rad/sec)=Ngfi zM 7t/30 where N g is the number of gear teeth and Q zM (rpm) is the 
driven shaft speed. By definition, the input power by the gear mesh elastic force 
F h (t)=k h e(t)=e h (co)e 1£0t is given by n s (co)= 1/2 (k h e h ) 2 Re { (l/Z s (co))* } . At a very high 
frequency, it is reasonable to assume that the boundaries of the shafts do not affect the 
gear mesh source regime. Hence the infinite shaft assumption should hold. 
Accordingly, Z s can be obtained analytically using the driving point shaft impedance 
given in Reference [77]; however, the gear mass m g must be included. 

Z s (co) =2p s A s c s (1 + i) + icom g (5.25) 


Sound power W radiated from the casing is computed using the following equation 
where \ is the total casing surface area 

-9 



194 


W(co) = z 0 A c <V^(co)>a c (co) (5.26) 

Assuming a source directivity Q((o) associated with geared system mounting condition, 
the sound pressure level Lp (dB re 20|iPa) in the far field at distance r from the casing is 

Q(co) 4 

L p (co) = L w (0» + 10 log 10 + R^) (5:27) 

where R(co)=aS/(l-a) is the room constant, a(co) is the average absorption 
coefficient, S is the room surface area and Lyy (dB re l.OpW) is the sound power level. 

5.7.4 Experimental Validation and Parametric Studies 

The SEA formulation for this example case is verified by comparing results with 
vibroacoustic responses measured on the NASA gearbox [55]. A detail description of 
the experimental setup and its system parameters have been given in Chapter IV. 
Additionally unweighted Lp at r « 0.38 m, directly above the surface of the top cover 
casing plate as shown in Figure 5.2, has been measured with a B&K type 2230 sound 
level meter. Predicted and measured Lp are compared in Figure 5.17a for y=10, 
a c =a c 2 and bearing reload 8 zm =0.04mm. The second term in equation (5.27) is 
dominated by Q/(4ni 2 ) since 4/R « Q/^Ttr 2 ) in this case due to a ~ 1.0 and room 
surface area S being very large. Predictions are found to be within ±10dB of the 
measured values for typical structural dissipation loss factor 0.004<r| c <0.04. Figure 

5.17b shows the comparison of spatially averaged mean-square casing acceleration level 
(or<V^>) between theory and experiment. Experimental curve represents the averaged 

value of the measurements made at 3 casing plate locations (top plate, side plate with 
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i(r i(r 10 


1/3 Octave Band Center Frequency (Hz) 
(a) 



1/3 Octave Band Center Freqeuncy (Hz) 


(b) 

Figure 5. 17 Comparison between theory and NASA experiment [55] for example case 
IV (y=10, Q=2, g c =g c 2 , 5 zm =0.04mm). (a) Sound pressure level L p . 
(b) Spatially averaged casing acceleration. 
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bearings and side plate without bearings). Again, predicted acceleration spectra, which 
are similar to that found for the sound level, are found to be in reasonable agreement 
with the measured spectrum given the appropriate values of r| c and y. In general, we 
observe that the response level decreases with increasing T| c . Also, comparisons 
suggest that r) c in this system may be frequency dependent. 

Next we vaiy y but keep r) c =0.02 constant in Figure 5.18. Comparison between 
theory and measurements also indicates that y=10 is the best fit for the experimental data 
especially at the higher frequencies. Here, radiation efficiency o c i and bearing axial 
preload ^zm — 0 .04mm have been specified. 

Now we investigate the effects of casing plate radiation efficiency O c and bearing 
preload or mean axial displacement 5 zm on Lp. Figure 5.19 compares a cl and a c2 . 
Based on the comparison with experiment it seems that the radiation efficiency of the 
NASA gearbox is better modeled with O cl although the measured Lp curve is mostly 
between O cl and a c2 curves. Differences between a cl and o c2 are significant at lower 
frequencies, but the variation never exceeds lOdB over the entire frequency range of 
interest. It may be noted that since the acoustic energy radiated W(co) is significantly 
smaller than the energy dissipated by the system, virtually no change is found in the 
predicted casing acceleration spectra by varying o c . The effect of bearing preload on Lp 
is shown in Figures 5.20a and 5.20b for a c2 and a c i respectively. Figure 5.20 
indicates that Lp is lower when the bearing preload is increased. A similar effect is seen 
for the casing acceleration level which is consistent with the deterministic model 
prediction of Chapter IV. 
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Figure 5.20 Effect of bearing axial preload 5 zm on L p for example case IV (n c =0.02, 
Q=2, 7=10). 
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5.8 CONCLUDING REMARKS 

The vibration transmission through bearing has been analyzed using the SEA 
technique. A new procedure has been developed to compute the coupling loss factor 
which relies on the solution of the boundary value problem at the bearing-casing 
interface. This scheme incorporates the bearing stiffness matrix developed earlier as a 
part of the deterministic s vibration models in Chapters H-IV. Experimental validations 
verify the proposed theory even though a very simple model for the geared rotor system 
has been presented. A more detailed SEA model including energy sharing between 
subsystems is required to analyze this system. This is left for future research. 


CHAPTER VI 


CONCLUSIONS AND RECOMMENDATIONS 
5.1 SUMMARY 

A new precision rolling element bearing model is proposed for the analyses of the 
vibration transmission through bearings and overall system dynamics in rotating 
mechanical equipment including geared rotor systems. Current bearing models, based 
on either ideal boundary conditions or purely translational stiffness element description, 
can not explain how the vibratory motion may be transmitted from the rotating shaft to 
the casing and other connecting structures in rotating mechanical equipment. For 
example, a vibration model of a rotating system based upon the existing bearing models 
can only predict purely in-plane type motion on the flexible casing plate given only the 
bending motion on the shaft. However, experimental results have shown that the casing 
plate motion is primarily flexural or out-of-plane type. Chapter II clarifies this issue 
qualitatively and quantitatively by developing a new mathematical model for the 
precision rolling element bearings from basic principles. A comprehensive bearing 
stiffness matrix [K]t> m of dimension 6 is proposed which clearly demonstrates a 
coupling between the shaft bending motion and the flexural motion on the casing plate. 
A numerical scheme which involves a solution to nonlinear algebraic equations is 
proposed for the estimation of the stiffness coefficients given the mean bearing load 
vector. And, a second method which requires the direct evaluation of these stiffness 
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coefficients given the mean bearing displacement vector is also discussed. Some of the 
translational stiffness coefficients of the proposed bearing matrix have been verified 
using available analytical and experimental data. Further validation of [K] bm is not 
possible as coupling coefficients are never measured. Also, parametric studies on the 
effect of unloaded contact angle, preload, or bearing type are included. These results 
lead to a complete characterization of the bearing stiffness matrix. 

Chapter III extends the proposed bearing stiffness formulation and demonstrates 
its superiority over existing models in vibration transmission analyses for a generic 
single shaft-bearing-plate-mount system. The bearing stiffness matrix [K]j, m is 
incorporated in discrete system models using lumped parameter and finite element 
modeling techniques. Shaft, plate and mount flexibilities are also included in such 
models. The stability issue associated with the proposed bearing model is addressed 
analytically using Liapunov's stability method and the system is found to be dynamically 
stable provided the preloads are sufficiently high. Eigensolution and forced harmonic 
response to the following rolling element bearing system example cases are obtained 
using our formulation and results are compared with the predictions yielded by the 
current vibration models: I. rigid shaft and plate system freely suspended, n. rigid shaft 
and plate supported on flexible mounts, and III. an experimental setup consisting of a 
flexible shaft, two ball bearings, a rectangular plate and the supporting structure. 
Analytical results indicate that our proposed model is indeed capable of predicting plate 
rigid body angular motion or plate flexural motion as excited by shaft motion. Such 
predictions are not observed in simple vibration models. Also, lower degrees of 
freedom models, developed by several previous investigators tend to underestimate the 
resonant frequencies and force or moment transmissibilities as compared with our multi- 



203 


degree of freedom models. Comparisons between our model and the available 
experiments have been found to be reasonably good. 

The overall dynamics of a geared rotor system which includes a spur gear pair, 
shafts, rolling element bearings, prime mover and load (attached to the geared rotor 
system through flexible torsional couplings), rigid or flexible casing, and compliant or 
massive mounts is studied in Chapter IV. Linear time-invariant, discrete dynamic 
models of a generic geared rotor system with proportional viscous damping are 
developed using lumped parameter and dynamic finite element techniques which are then 
used to predict the vibration transmissibility through bearings and mounts, casing 
vibration motion, and dynamic response of the internal rotating system. The proposed 
bearing formulation is also included. Each rotating shaft is modeled as an Euler beam in 
the lumped parameter model and as a Timoshenko beam in the dynamic finite element 
model, but the gyroscopic moment is not included. Eigensolution and forced harmonic 
response studies due to rotating mass unbalance or kinematic transmission error 
excitation for the following example cases are obtained using our formulation and 
compared with simple models currently available in the literature and/or experiment: I. a 
single-stage rotor system with flexibly mounted rigid casing consisting of two bearings 
as a special case of the geared rotor system, n. a spur gear pair drive supported by four 
bearings installed in a flexibly mounted rigid casing, and III. an experimental setup 
consisting of high precision gear and pinion, and four identical rolling element bearings 
contained in a flexible casing mounted rigidly on a massive foundation. In example 
cases II and III of Chapter IV, the gear mesh stiffness is assumed to be linear and time- 
invariant. Analytical predictions show that our theory is indeed capable of predicting 
bearing and mount moment transmissibilities in addition to the force transmissibilities. 
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Also, flexural vibrations of the casing plate are predicted well as our theory is in good 
agreement with measurements made on case in of Chapter IV; such predictions are not 
seen in simple models due to inadequate bearing formulation. 

In Chapters n-IV, the deterministic vibration models of geared rotor system with 
proposed bearing stiffness matrix incorporated have been shown to be reliable up to 
moderately high frequencies. But these models are inadequate at very high frequencies 
when the modal density is high. Classical vibration models usually do not predict 
modes accurately in this frequency regime, require large computational effort and 
produce spectra at many spatial points which are difficult to interpret. Accordingly, we 
have used the statistical energy analysis (SEA) method to predict the mean-square 
vibratory response of internal and external subsystems in a geared rotor system. The 
feasibility of applying SEA to this system is also investigated by performing modal 
analysis of a gearbox. Four example cases are analyzed using SEA: I. a plate- 
cantilevered rectangular beam, n. case I with circular shaft-bearing system replacing the 
cantilevered rectangular beam, HI. a circular shaft-bearing-plate-mount system, and IV a 
simple geared rotor system. In the first two example cases, we have revised and 
extended Lyon and Eichler's plate-cantilevered rectangular beam problem [56,58] to 
improve the coupling loss factor prediction and to formulate the vibratory energy 
transfer problem through rolling element bearings. The third and final example cases 
compute the system response spectra and compare them with measurements. Good 
agreement is found between theory and experiment provided proper values of the 
dissipation loss factors and bearing preloads are used. The NASA gearbox radiation 
efficiency is found to be nearly unity at higher frequencies. 
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5.2 CONTRIBUTIONS 

In this dissertation, several contributions related to the bearing and geared rotor 
systems vibration studies are evident: 

1 . A comprehensive bearing stiffness matrix of dimension 6 have been proposed and 
developed for the analyses of vibration transmission through precision rolling 
element bearings and overall geared rotor system dynamics. Analytical expressions 
for the bearing stiffness coefficients have been given for several types of ball and 
roller bearings. Also, two methods have been proposed for the estimation of the 
bearing stiffness coefficients which involves solution to nonlinear algebraic 
equations and direct evaluation of these coefficients, respectively. Using 
Liapunov's stability method, the proposed bearing system models have been shown 
to be stable provided preloads are sufficiently high. 

2. The proposed bearing model clearly demonstrates a coupling between the shaft 
bending motion and the flexural motion on the casing plate through three example 
cases including comparison with available experiments. Such predictions are not 
observed in existing bearing vibration models. An overall geared rotor system 
vibration model with the proposed bearing formulation incorporated have been 
developed analytically and verified experimentally. This system includes spur gear 
pair, flexible shafts, flexible or rigid casing, motor, load, and compliant or massive 
mounts. Effects of casing and mounts on the internal geared system vibration and 
on the overall system behavior have been evaluated. Bearing moment transfer 
spectra in addition to force transfer spectra have been computed when the proposed 


206 


bearing formulation is included in tire system models. Simple bearing model can 
only predict bearing force transmissibility spectra in such cases. 

3 . Analytical expressions for the coupling loss factor, a statistical energy analysis 
(SEA) parameter, of a shaft-compliant bearing-plate problem have been developed 
which include the newly formulated bearing stiffness matrix. The coupling loss 
factor formulation for a plate-cantilevered rectangular beam problem developed 
earlier earlier by Lyon and Eichler [56,58] has been revised. The proposed 
prediction agrees better with the experimental prediction given by Lyon and Eichler 
[58] than the previous approximate formulation. The mean-square vibroacoustic 
response of a geared rotor system predicted using SEA at high frequencies has been 
validated by comparison with experiment. 

5.2 FUTURE RESEARCH 

Several areas of potential research problems based on the present study of the 

vibration transmission through bearings, are identified as follows: 

1. Analyze vibration transmission through hydrodynamic bearings, using the 
proposed work on rolling element bearings. Develop an experimental methodology 
to estimate the bearing stiffness matrix and transfer properties. Also, develop 
bearing diagnostic techniques using vibration transmission theory proposed in this 
dissertation. 
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2. Extend the proposed overall geared rotor system vibration model with spur gears to 
helical, bevel and worm gear drives. Develop vibration models using the proposed 
bearing model to analyze multiple transmission paths. Improve gear diagnostic 
procedures using analytical bearing transfer properties. 

3. Refine the current statistical energy analysis formulation of a generic geared rotor 
system to investigate the vibratory energy transfer among smaller subsystems and 
to include the effects of rotating shaft torsional modes. Generate analytical and 
experimental schemes to predict coupling loss factors for typical shaft-bearing-plate . 
system. Finally, establish gearbox design methodology for reduced airborne and 
structure-borne noise. 
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